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Foreword 


Determining  the  scattered  field  from  objects  in  a  waveguide  has 
been  difficult  and  is  of  great  importance  to  the  U.S.  Navy. 
Since  future  detection  methods  will  exploit  for  both  long-  and 
short-range  bistatic  detections  of  submerged  objects,  determining 
whether  various  system/detection  scenarios  are  feasible  will 
depend  on  both  experiments  and  numerical  model  studies.  Numerical 
studies  are  particularly  appealing,  since  they  are  considerably  less 
expensive  than  experiments  and,  in  terms  of  a  time  frame,  many 
numerical  studies  can  be  done  for  any  given  experiment.  More¬ 
over,  numerical  studies  can  also  be  used  to  suggest  experiments. 

This  report  proposes  a  numerical  scheme  that  will  adequately 
describe  scattering  from  realistic  objects  (submarines)  in  an  ocean 
waveguide.  The  type  of  object  used  in  this  report  was  a  spheroid, 
and  the  object  was  assumed  to  have  Dirichlet  boundary  conditions. 
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Executive  Summary 


A  method  was  developed,  tp  describe  acoustical  scattering  from  an 
object  in  a  waveguide  by  using  normal  mode  theory  to  describe  the  incident 
Held.  Each  mode  is  decomposed  at  the  object  into  a  pair  of  upgoing  and 
downgoing  plane  waves.  A  transition  matrix  was  used  (developed  via 
the  extended  boundary  condition  method  of  Waterman)  to  determine  the 
resulting  near-field  scattered  field.  The  far-field  scattered  field  was 
determined  by  invoking  Huygens’  principle.  This  far-field  solution  satisfies 
all  boundary  conditions  and  preserves  continuity  of  the  solution 
throughout  all  space. 

The  examples  show  that  the  object  is  correctly  coupled  to  the  waveguide. 
This  was  done  by  showing  that  the  object’s  scattered  field  acts  as  a 
secondary  source  and  that  this  scattered  field  obeyed  the  same  boundary 
conditions  as  the  point  source  field. 

This  method  of  determining  the  scattered  field  from  a  three-dimensional 
object  allows  one  to  determine  not  only  the  correct  target  strength  (intensity), 
but  also  to  properly  determine  the  phase.  In  addition  this  method  allows 
the  investigation  of  the  interaction  between  the  incident  field  and  the 
waveguide,  the  incident  field  with  the  object,  and  the  object’s  scattered 
field  with  the  waveguide.  In  this  manner  a  better  understanding  of  the 
ongoing  physical  processes  can  be  obtained. 
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Theory  and  Application  of  Scattering  from 
an  Object  in  an  Ocean  Waveguide 


Synopsis 

The  treatment  of  scattering  from  submerged 
objects  in  an  unbounded  environment  is  of 
considerable  interest  to  both  the  academic  and 
technological  communities.  Several  approaches  have 
yielded  results  for  different  classes  of  problems  and 
have  proven  manageable  for  the  free  environment 
case.  The  solution  of  scattering  from  objects  in  a 
waveguide  has  been  a  difficult  and  almost  intractable 
problem  to  do  in  an  exact  framework,  due  to  the 
coupling  of  effects  from  the  object’s  scattered  field 
with  that  of  the  boundaries.  It  is,  however,  a  problem 
that  is  of  great  importance  to  the  U.S.  Navy 
community,  since  future  detection  methods  are  going 
to  exploit  both  long-  and  short-range  bistatic 
detections  of  threat  submarines.  The  ability  to 
determine  whether  various  system/detection  scenarios 
are  feasible  will  depend  on  both  experiments  and 
numerical  model  studies.  Numerical  studies  are 
particularly  appealing,  since  they  are  considerably 
less  expensive  than  experiments  and  in  terms  of 
time,  many  numerical  studies  can  be  performed  for 
any  given  experiment.  Moreover,  numerical  studies 
can  be  used  to  suggest  experiments.  The  objective 
of  this  work  is  to  develop  a  model  suitable  for  such 
studies  in  an  approximate  but  accurate  framework. 


Summary 

The  study  described  in  the  appendix  presents  a 
coherent,  self-consistent  method  to  describe  acous¬ 
tical  scattering  and  propagation  from  an  object  in  a 
waveguide.  The  method  is  based  on  Huygens’  prin¬ 
ciple.  The  object  was  a  spheroid  with  Dirichlet 
boundary  conditions.  The  study  began  by  developing 
the  extended  boundary  condition  (EBC)  method*'^ 
to  describe  acoustical  scattering  from  three- 
dimensional  objects  in  free  space.  This  method  was 
used  to  describe  scattering  from  objects  with  both 


Neuman,  and  Dirichlet  boundary  conditions  and  has 
been  extended  for  elastic  objects.  Next,  the  theory 
of  acoustic  propagation  in  a  fluid  was  developed 
from  first  principles  and  both  normal  mode  and  ray 
theory  were  developed  to  describe  propagation  in 
the  simplest  type  of  acoustic  waveguide. 

The  mathematical  expression  that  describes  the 
object’s  far-field  scattered  field  is  developed  here 
through  the  use  of  Huygens’  principle.  The  near¬ 
field  scattered  field,  developed  by  allowing  a 
transition  matrix  (developed  by  the  EBC  method) 
to  map  the  incident  field  onto  the  scattered  field, 
was  coupled  to  the  Green’s  function,  which  satisfies 
the  waveguide.  This  application  of  Huygens’ 
principle  leads  to  a  manageable  direct  solution  of 
the  problem.  This  method  also  satisfies  all 
appropriate  boundary  conditions  and  yields  a 
continuous  solution  throughout  space.  This 
expression  is  then  applied  to  a  Pekeris  waveguide 
and  finally  to  a  more  general  range-independent 
waveguide. 

The  examples  presented  in  the  last  two  chapters 
of  the  appendix  were  designed  to  demonstrate  that 
the  objea’s  scattered  field  is  correctly  coupled  to  the 
waveguide.  This  was  accomplished  by  showing  that 
the  object  acted  as  a  secondary  source  in  the 
waveguide.  Its  field  obeyed  (satisfied)  the  same 
boundary  conditions  for  the  waveguide  as  the  point 
source  field. 


Recommendations 

The  results  obtained  are  for  a  waveguide  that  is 
range  independent;  that  is,  the  sound  speed  profile 
in  the  water  column  and  the  sediment  vary  only  in 
depth  and  not  in  range,  and  the  bathymetry  does 
not  change  with  respect  to  range.  In  addition,  the 
object  is  expected  to  be  in  the  far  field  of  the  source. 
The  model  could  simulate  range  dependency  by  cither 


applying  the  adiabatic  approximation  or  by  coupling 
the  discrete  (propagating)  modes.  Additionally,  if  the 
continuous  spectrum  of  the  incident  field  is  included 
and  not  just  the  discrete  spectrum,  then  the  object 
could  be  located  near  the  source  where  the  continuous 
spectrum  could  play  an  important  role.  The  addition 
of  range  dependency  would  be  the  most  valuable 
addition  to  the  model. 

The  incident  field  was  allowed  to  interact  only 
once  with  the  object;  that  is,  the  primary  interaction 
of  the  incident  field  with  the  scatterer  was  retained, 
and  the  subsequent  interactions  of  the  scattered  field 
between  the  boundaries  of  the  waveguide  and  the 
scatterer  were  ignored.  This  is  not  a  limitation,  since 


the  resulting  field  would  consist  of  high-angle  energy 
and  would  ^ssipate  quickly  in  an  ocean  environment. 
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CHAPTER  1 


INTRODUCTION 

Goals  of  Research 

This  research  project  will  answer  the  following  question.  What  is  the 
acoustic  pressure  at  a  field  point  in  an  underwater  acoustic  waveguide  when  a 
three  dimensional  object  is  present?  This  research  project  concerns  classical 
wave  propagation,  diffraction  and  scattering  effects.  An  inherent  assumption 
is  that  the  object's  dimension's  (length  and  width)  are  small  compared  to  the 
waveguides  dimensions.  We  therefore  arc  assuming  that  the  object  will  make  a 
small  perturbation  to  the  acoustic  field  produced  by  a  monochromatic  point 
source.  We  assume  then  that  the  wave  equation  will  properly  describe  the 
acoustic  wave  propagating  in  the  waveguide.  We  assume  also  that  the  field  can 
be  expressed  as  a  near  field  scattered  field  and  a  far  field  scattered  field.  The 
near  field  scattered  field  will  only  be  valid  near  the  object,  before  it  interacts 
with  the  waveguide  and  the  far  field  scattered  field  will  be  appropriate  in  the 
far  field  of  the  object,  that  is  when  the  scattered  field  has  interacted  with  the 
waveguide  boundaries.  With  these  assumptions,  this  project  can  be  subdivided 
into  four  distinct  sub  areas.  First  the  incident  field  must  be  determined,  that  is 
the  acoustic  field  in  the  absence  of  the  object  or  scatterer.  Second,  the 
incident  field  which  insonifies.  the  object  is  mathematically  acted  on  by  a 
transition  matrix  which  transforms  the  incident  field  into  the  near  field 
scattered  field.  Next,  the  near  field  scattered  field  is  coupled  to  the  waveguide 
via  invoking  Huygens'  principle  (see  chapter  4).  Finally,  we  obtain  the  total 
field  by  coherently  adding  the  incident  and  scattered  field. 
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Because  of  ihc  way  the  problem  divides  mathematically,  we  have  divided 
this  discussion  into  two  parts.  The  first  is  concerned  mainly  with  theory  and 
results  pertaining  to  free  space  scattering  (chapter  2)  and  acoustic 
propagation  in  a  waveguide  (chapter  3).  This  is  done  in  order  to  first  define 
the  physics  that  defines  scattering  from  an  object  (via  a  matrix  mapping)  and 
the  propagation  of  this  and  the  incident  field  in  a  waveguide.  In  part  2 
(chapters  3-6)  of  this  dissertation  we  present  the  various  methods  that  were 
used  prior  to  the  present  work.  We  then  present  formally  this  new  method 
based  on  Huygens'  principle  (chapter  4)  and  in  Chapters  5-6  apply  the  method 
to  two  acoustic  waveguides.  The  first  is  a  waveguide  composed  of  an  isovelocity 
water-layer  over  an  isovelocity  fluid  half-space  and  the  second  is  a 
multilayered  waveguide  consisting  of  multiple  isovelocity  sound  and  density 
layers  over  an  isovelocity  half  space. 

Appendix  A  validates  the  Extended  Boundary  Condition  (see  chapter  2) 
code  by  showing  a  comparison  of  the  scattered  field  generated  by  this  method 
with  the  analytical  solution  for  the  case  when  a  plane  wave  is  incident  upon  a 
sphere.  Appendix  B  shows  a  derivation  of  a  result  needed  in  the  development 
of  the  acoustic  field  in  a  waveguide.  Appendix  C  is  a  validation  of  the 
propagation  models  used  in  this  dissertation.  Appendix  D  is  a  brief  description 
of  the  attenuation  term  used  in  the  two  waveguides.  Appendix  E  shows  the 
comparison  of  the  scattered  field  using  the  proposed  method  based  on 
Huygens'  principle  along  with  a  different  method  described  in  Chapter  4. 

These  results  were  presented  by  the  author  at  the  Second  Joint  Meeting  of  the 
Acoustical  Society  of  America  and  the  Acoustical  Society  of  Japan  14-18,  Nov. 


1988  in  Honolulu. 


CHAl^ER  2 


FREE  SPACE  SCATTERING 

In  this  chapter  we  give  a  brief  overview,  pf  the  usual  methods  used  to 
describe  acoustical  scattering  prior  to  the  introduction  of  the  Extended 
Boundary  Condition  (EBC)  method.  We  then  develop  the  EBC  method,  which  is 
used  in  this  dissertation  to  describe  the  near  scattered  field.  We  complete  this 
chapter  by  showing  some  representable  free  space  scattering  results  from  this 
method  for  both  far  and  near  field.  The  references  listed  represent  a  fraction 
of  the  work  being  done  by  the  scattering  community  using  this  method  and  is 
not  meant  to  be  all  inclusive. 


Early  Free  Space  Scattering  Techniques 

Historically  there  are  three  predominate  methods  employed  to  solve 
scattering  and/or  diffraction  problems.  They  are  separation  of  variables, 
variational  techniques,  and  the  direct  numerical  solution  of  integral 
equations. 

The  separation  of  variables  technique  has  been  used  for  objects  bounded 
by  quadric  (separable)  surfaces.  The  technique  can  be  used  as  long  as  the 
scalar  wave  equation  governing  the  problem  can  be  solved  by  separation  of 
variables.  This  method  solves  the  wave  equation,  (a  linear  partial  differential 
equation)  by  separating  the  original  equation  into  a  set  of  ordinary 
differential  equations,  each  involving  only  one  variable.  The  general  solution 
is  a  product  of  the  solutions  to  the  ordinary  differential  equations.  For 
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example,  consider  a  spheroidal  object.  One  can  use  the  spheroidal  coordinate 
system  to  re*  sent  the  incident  and  scattered  wavefunctions  by  a  series  of 
spheroidal  wavefunctions.  This  requires  that  the  wave  equation  be  separated 
in  the  spheroidal  coordinate  system.  This  approach  has  been  carried  out  by 
Burke,  who  has  solved  the  scattering  problem  for  a  plane  wave  incident  upon 
an  impenetrable  spheroid  having  Dirichlet  boundary  conditions  (1)  (which 
fixes  the  value  of  the  wavefunction  on  the  surface)  as  well  as  for  a  rigid 
spheroid  having  Neuman  boundary  conditions  (2)  (which  fixes  the  value  of 
the  normal  derivative  of  the  wavefunction  on  the  surface).  He  has  also  solved 
the  problem  of  a  plane  wave  incident  upon  a  penetrable  spheroid.  (3)  Yeh  has 
solved  the  problem  for  a  plane  wave  incident  upon  a  penetrable  liquid  prolate 
spheroid. (4)  The  liquid  prolate  spheroid  is  the  limiting  case  of  an  elastic 
spheroid  with  zero  shear  modulus.  The  main  reason  this  approximation  w'as 
adopted  was  because  the  vector  wave  equations  are  not  separable  in  spheroidal 
coordinates.  The  secondary  reason  being  that  due  to  the  complexity  of 
spheroidal  coordinates,  laborious  computations  had  to  be  carried  out  in  order 
to  obtain  numerical  values.  There  are  other  coordinate  systems,  one  could  use 
to  solve  various  scattering  problems,  spherical,  cylindrical,  etc.,  depending  on 
the  type  of  object  or  objects  one  is  concerned  with.  Considering  only  confocal 
quadric  surfaces  and  more  specifically  ellipsoids,  there  are  nine  coordinate 
systems  which  can  be  obtained  front  ellipsoidal  coordinates.  (5)  To  solve  the 
scattering  problem,  a  large  part  of  the  computational  effort  goes  into  the 
evaluation  of  the  wavefunctions  themselves  except  for  the  sphere  and  the 
circular  cylinder,  for  which  efficient  recursion  relations  arc  available.  (5) 
Levine  and  Schwinger  (6)  described  the  variational  method  for  the 
diffraction  of  a  scalar  plane  wave  by  an  apeiture  in  an  infinite  plane  screen. 
In  describing  diffraction  from  general  bodies  the  main  effort  goes  into 
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evaluating  matrix  elements,  which  consist  of  repeated  surface  integrals 
requiring  fourfold  numerical  quadrature  or  volume  integrals  requiring 
sixfold  numerical  quadrature.  In  each  case  the  integrals  have  singular 
kernels. 

The  integral  equation  method  has  the  advantage  of  generality,  for 
usually  the  integral  is  invariant  under  coordinate  transformation.  Once  the 
Greens  function  is  found  for  the  appropriate  geometry,  the  solution  of  the 
scalar  wave  equation  may  be  found.  The  problem  with  this  method  is  that  for 
many  cases  the  integral  cannot  be  integrated  in  closed  form  and  numerical 
values  are  then  extremely  difficalt  to  obtain.  However,  with  the  advent  of  the 
digital  computer,  this  method  has  been  used  by  various  individuals.  (7-12) 

More  recently  the  integral  equations  have  been  classified  into  two  types, 
depending  upon  the  observation  points.  (13)  Surface  integral  equations  have 

the  observation  points  lying  on  the  surface  of  the  scatterer  and  the  Extended 
integral  equations  have  the  observation  points  inside  the  scatterer. 

Waterman  (14)  described  in  his  1969  paper  a  then  new  matrix 
formulation  of  acoustic  scattering.  He  had  previously  developed  the  same 
method  for  electromagnetic  scattering.  (15)  The  equations  generated  most 
nearly  resembled  those  obtained  using  the  variational  method.  The  advantage 
was  that  for  both  surface  and  volume  type  scattering,  elements  of  the  matrix  to 
be  inverted  are  described  by  a  single  surface  integral  with  no  singularities  in 
the  integrand.  Waterman  used  the  interior  boundary  equation  as  a  constraint 
in  solving  the  exterior  boundary  equation.  As  he  stated  in  the  paper  this 
approach  was  not  necessarily  new  as  it  appears  that  Smythe  applied  it  in 
electrostatics  in  1956.  (16)  The  next  section  presents  the  derivation  of  the 
scattering  matrix  using  this  method  derived  by  Waterman.  The  method  of 
obtaining  the  scattering  matrix  has  been  known  by  various  names,  T-Matrix 


6 


method.  Extended  Boundary  Condition  method  and  the  Null  Field  Method.  We 
believe  that  the  Extended  Boundary  Condition  method  best  describes  this 
method.  The  method  has  since  evolved  and  is  used  in  many  areas  of  classical 
scattering  determinations.  (15,17-22) 

The  Extended  Boundary  Condition  (EBC)  Method  to  Determine  the  Scattered 

Eidi 


The  theoretical  basis  for  the  Extended  Boundary  Condition  method  is 
contained  in  the  mathematical  formulation  of  Huygens'  principle.  (23)  This 
description  follows  not  only  Watermans’  but  also  Werby  and  Chin-Bing.  (24) 

We  start  with  the  Helmholtz-Poincare'  integral  representation  of  the  total  field 
Ut  exterior  to  the  bounded  object,  which  has  the  following  form,  (23) 


(2-1) 


Uj(r)  =  U.(r)  + 


3Gtr  aU  ,  (r') 


ds 


where  U+(r')  is  the  scalar  wavcfield  on  the  object  surface,  G(r,r’)  is  the 
outgoing  Green's  function,  r'  is  taken  on  the  surface  of  the  object,  and  n  is  a 
unit  vector  normal  to  the  surface  of  the  object.  The  surface  G  is  taken  to  be 
the  surface  of  the  bounded  object.  In  order  to  obtain  a  unique  solution,  we  will 
need  an  additional  expression.  We  will  use  the  field  interior  to  the  object.  This 
Helmholtz-Poincare’  expression  is  as  follows,  (23) 


(.:-2) 


0  =  U.(r'' 
1 


4. 


U,(r-)-^ 


,r')  aU^(r') 

-G(r%r')-gV- 


ds 
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where  r"  is  a  point  taken  in  the  interior  of  the  object.  This  equation  was 
recognized  by  Waterman  as  providing  a  constraint  to  eliminate  the  unknown 
surface  quantities  U+,  These  two  equations  yield  the  extended  boundary 
condition  equations.  (25)  These  equations  will  now  be  put  into  a  form 
amenable  to  numerical  computation.  For  this  example  we  allow  the  object  to  be 
impenetrable,  which  requires  that 

(2-3)  U+=0 . 


Equation  (2-1)  becomes 


(2-4) 


U  fr)=U.(r)- 
t  1 


G(r,r') 


au+(r') 
in  . 


ds 


and  Eq.(2-2)  becomes 


(2-5) 


0  =  U.(r'')- 
1 


G(r%r') 


3U+(r') 

Ti . . 


ds 


In  order  to  solve  these  expressions,  it  is  convenient  to  express  Ui(r),  U+(r')  and 
G(r,r')  as  some  suitable  series  expansion,  which  upon  truncation  leads  to 
matrix  equations  that  can  be  solved  using  digital  computers.  The  Green's 
functions,  G(r,r')  is  a  normal  operator  and  can  be  expressed  by  the  following 
biorthogonal  series  (26) 


G(r,r')  =ik  SR e9.(r  ^)(p.(r  >) 


(2-6) 


where  r<  and  r>  is  the  lesser  or  greater  of  the  two  points  r  and  r’  relative  to 
the  origin  of  the  object.  The  incident  wavefield  Ui  is  known  and  can  be 
expressed  as,  (27) 


(2-7) 


where  Re  is  the  regular  part  of  Cpnfr)-  The  fact  that  this  expansion  Eq.  (2-7) 
can  be  obtained  follows  from  the  Hilbert-Schmidt  theorem. (27)  Using  this  and 

Eq.(2-6),  we  see  that  Eq.(2-5)  becomes. 


(2-8) 


au.(r') 

Sa^Re(Pn(r)  =  iklJ  Re(p^^(r)(Pn(rO— ^-5  -  ds 
n  no 


(2-9) 


=  ik 


au^.(r') 

5n 


ds 


U+(r')  is  now  written  in  some  complete  set  of  suitably  chosen  basis  function 
satisfying  completeness  on  the  object  surface.  The  basis  functions  used  were 
the  same  as  those  chosen  by  Waterman,  namely  Re  (pn(r).  where  Re  is  the 
regular  part  of  <pn(r)*  We  have  not  found  any  problem  arising  from  irregular 
values,  nor  have  we  had  problems  due  to  poorly-convergent  solutions. 
Waterman  has  shown  that  Eq.  (2-10)  satisfies  closure  on  the  surface  of  a  rigid 
object: 


Using  Eq.  (2-10)  in  Eq.  (2-9)  results  in  the  following  for  the  incident  field 
expansion  coefficients  an 


(2-11) 


a  =  i  k  Xb  _ 
n  ni 


m  •'a 


aRecp  (r') 
n-— 


(2-12) 


=  ik 

m 


where 


■■mn 


3  Recp  (r') 

<P  (r') - - ds  .  In  matrix  notation,  Eq.  (2-11)  can  be 

^  3n 


written 


(2-13) 


a  =  ik  Qb 


and  the  expression  for  Us(r')  becomes. 


(2-14) 


U,(r')=a„P„(r') 


Recognizing  that  Us  =  Ut  -  Uj,  the  expression  for  the  scattered  field  expansion 
coefficients  fp,  becomes. 


(2-15) 

(2-16) 


=  ~ik  Xb_ 

*"  n  '  '  tr  m 
n  n 


3Req)jj^(r') 


Re(p„(r)T^(r') - - ds 


3Rc<p  (r') 
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where 


ReQ, 


mn 


-J, 


Re(Pn(r) 


aRe9„(r') 

dn 


In  matrix  notation,  fn  becomes, 


(2-18) 


f  =-ikb»eQ 


Equation  (2-13)  can  be  solved  for  b, 

(2-19)  b  =  ;:i^ 

k 


Inserting  Eq.(2-19)  in  Eq.(2-18)  one  obtains, 

(2-20)  f=-(ReQ)Q“S 

o  r 

(2-21)  f  =  Ta 

where  T  =  “(ReO  Q“^  .  Notice  that  the  matrix  T  is  only  a  function  of  the 
boundary  conditions  and  the  shape  of  the  object.  Once  the  matrix  T  is  known, 
the  scattered  field  can  be  determined  from  any  chosen  incident  field  a.  The 
Extended  Boundary  Condition  method  yields  a  unique  solution  to  the  exterior 
acoustic  problem,  and  is  efficient  and  numerically  stable  for  a  large  variety  of 
bounded  object  shapes.  Tne  Extended  Boundary  Condition  method  has  evolved 
since  its  inception  in  1965.  In  the  area  of  acoustic  scattering,  Werby,  (28) 
Werby,  Tango  and  Green,  (29)  and  Werby  and  Chin-Bing  (24)  have  made 
significant  improvements  to  the  general  method.  Not  only  is  this  method 
applicable  to  solid  surface  of  revolutions  but  also  to  shells  (30-36)  and  other 
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clastic  objects.  (37*41)  Resonance  effects  are  observed  and  can  be  predicted 
using  this  method.  (33,36,39,42,43) 

Far  Field  Scattering  Examples  Using  the  Emnded  Boundary  Condition  M^J hod 

In  this  section  we  will  show  some  examples  of  the  scattered  field  in  the 
form  of  scattering  amplitude  vs.  angle,  normally  called  the  bistatic  angular 
distribution.  (39)  Bistatic  angular  distributions  are  defined  by  the  case  of  non¬ 
coincident  source  and  receiver,  where  the  receivers  are  omni  directional. 
Bistatic  angular  distributions  are  dependent  on  object  geometry  and  can  be 
useful  in  determining  such  features  of  object  shape  as  symmetry  or 
elongation.  The  far  field  bistatic  angular  distribution  is  defined  as  the 
following,  (38) 


(2-22) 


+  i  skp 
i)  e 


p™(cot(e))f;^.A“(0.) 


J 


where  A™^0.^is  the  incident  field  expansion  coefficient,  which  is  dependent 

upon  the  angle  at  which  the  source  insonifies  the  object,  a  is  the  semi  major 

m 

axis  of  the  object  and  f^^,  is  the  free  space  scattering  field  as  determined  by 

the  EBC  method  (see  Eq.  (2-21)).  Thcicforc  the  far  field  bistatic  angular 
distribution  is  dependent  upon  two  angles  (incident  and  observed),  and  not 
just  one. 

An  example  of  the  characteristics  of  a  typical  scattering  object  are 
shown  in  Table  1.  The  incident  field  was  a  monochromatic  point  source.  It  is 
assumed  that  the  object  is  sufficiently  far  from  the  point  source  so  that  the 
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incident  wave  front  has  infinite  radius  of  curvature  that  is,  it  is  a  plane  wave. 
The  three  frequencies  chosen  are  given  in  Table  2,  along  with  other 
information. 


TABLE  1 

OBJECT  CHARACTERISTICS 


Object  Type 

Aspect  Ratio 

Boundary 

Length(L) 

Width(W) 

(L/W) 

Condition 

(m) 

(m) 

Solid  Spheroid 

5 

Dirichlet 

50. 

10. 

TABLE 2 

EXTENDED  BOUNDARY  CONDITION  INFORMAllON 


Frequency 

(Hz) 

KL/2 

Sound  Speed 
(m/s) 

% 

(m) 

Highest 

Matrix  Order 

#  Matrix 

95.5 

10. 

1500.0 

15.7 

17 

4 

191.0 

20. 

1500.0 

7.9 

32 

7 

382.0 

40. 

1500.0 

3.9 

64 

13 

Figure  1  compares  the  angular  distribution  with  the  incident  direction 
along  the  major  axis  of  symmetry.  The  object  is  oriented  along  the  O-ISO 
degree  axis  in  the  figure  and  the  incident  field  is  coming  from  the  direction  of 
the  arrow.  Angular  distribution  plot  (a)  is  for  a  frequency  of  95. .5  (Hz)  or  a 
KL/2  of  10.  Where  K  is  tlic  wavenumber  and  L/2  is  half  the  object's  length. 

The  corresponding  wavelength  is  1.*^  ‘  (m).  The  object  is  approximately  3 
wavelengths  long.  This  would  represent  a  low  frequency  case.  Note  the 
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broadness  of  the  field.  Even  at  this  low  frequency  the  field  is  predominantly 
in  the  forward  direction.  As  we  increase  in  frequency  (and  KL/2)  we  note  that 
backscattering  (toward  the  source)  is  reduced,  while  it  is  being  focused  in  the 
forward  direction.  The  amplitude  is  also  increasing  with  increasing 
frequency. 

Figure  2  illustrates  the  case  when  the  object  is  insonified  at  an  angle  of 
45  degrees  relative  to  the  major  axis  of  symmetry.  The  top  plot  is  for  a 
frequency  of  95.5  (Hz)  and  the  bottom  is  for  382.  (Hz).  Note  that  for  the  low 
frequency  case,  the  field  is  oscillatory  in  appearance.  We  also  see  two  peaks, 
one  in  the  forward  direction  (225  degrees)  and  the  other  in  the  specular 
direction  (135  degrees).  The  peaks  become  more  focused  in  the  forward 
direction  as  we  increase  the  frequency.  The  forward  diffracted  wave  becomes 
more  highly  focused  with  increased  frequency.  Also  note  that  the  wave-like 
appearance  of  the  field  at  low  frequency  virtually  disappears  in  the  high 
frequency  case.  Again  we  see  that  the  field  tends  to  increase  in  amplitude  as 
frequency  increases. 

Figure  3  illustrates  the  case  when  the  object  is  being  insonified  at  an 
angle  of  90  degrees  or  broadside  incidence.  Note  for  the  low  frequency  case, 
the  forward  and  backscattered  amplitude  are  about  equal.  As  we  increase  in 
frequency  we  see  that  the  forward  diffracted  wave  becomes  more  focused  and 
its  amplitude  is  larger  than  the  backscattered  (90  degrees)  wave.  In  fact  we 
see  that  the  backscattered  amplitude  remains  approximately  constant.  This  is 
because  the  backscattered  amplitude  (at  90  degrees  incident  and  90  degrees 
observed)  is  dependent  upon  geometrical  effects  and  is  independent  of 
frequency.  Now  we  turn  our  attention  to  some  near  field  examples. 
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270 


Fig.  1.  Angular  Distribution  of  Free  Field  Fcm  Function  for  an  Objwt  Being 
insonified  Along  the  Principle  Axis  of  Symmetry  for  KL/2  Values  of  10, 20.  and  40. 
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R-equency  •  95^  (Hz) 
KL/Z-10.0 


Frequency  •  191.0  (Hz) 
10/2-20.0 


90 


Frequency  •  382.0  (Hz) 
KL/2-40.0 


Fig.  2.  Angular  Distribution  of  Free  Field  Forai  Function  for  an  Object  Being  Insonlfied 
45  Degrees  Relative  to  the  Principle  Axis  of  Symmcuy  for  KU2  Vtdues  of  10. 20,  and  40. 
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Frequency  •  95^  (Hz) 
KU2.10.0 


Frequency  •  191.0  (Hz) 
KU2-20.0 


180  h-40-20  g;?;^20-40H  0 


Frequency  •  382.0  (Hz) 
KL/2.40.0 


Fig,  3.  Angular  Distribution  of  Free  Held  Form  Function  for  an  Object  Being  Insonified 
90  Degrees  Relative  to  the  Principle  Axis  of  Synunetiy  for  KL/2  Vdues  of  10, 20,  and  40. 
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Near  Field  Scattering  Examples  Using  the  Extended  Boundary  Condition  Method 


In  this  section  we  will  show  representable  examples  of  the  near  field 
form  function.  The  parameters  of  Tables  1  and  2  from  the  last  section  will  be 
used  in  this  section.  The  near  field  form  function  takes  the  following  form. 


(2-23) 


The  mathematical  surface  on  which  the  form  function  will  be  evaluated  is  a 
spheroid  with  an  aspect  ratio  of  5:1  a  semi-major  axis  of  125(m)  and  a  semi¬ 
minor  axis  of  25(m).  This  surface  is  not  the  surface  of  the  object  but  a 
mathematical  surface  enclosing  the  object. 

Figure  4  illustrates  the  field  for  3  different  frequencies  as  it  is  being 
insonified  along  the  axis  of  symmetry.  Note  for  the  low  frequency  case  (top 
figure)  the  field  is  oscillatory  in  appearance  while  in  the  intermediate  (middle 
figure)  and  high  (bottom  figure)  frequency  cases  the  field  is  more  strongly 
focused  in  the  forward  direction.  The  amplitude  of  the  field  increases  with 
increasing  frequency. 

Figure  5  illustrates  the  case  when  the  object  is  being  insonified  at  an 
angle  of  45  degrees  relative  to  the  principle  axis  of  symmetry.  We  see  the  same 
traits  as  in  the  previous  figures,  with  the  low  frequency  case  appearing  wave¬ 
like  while  the  intermediate  and  high  frequency  cases  are  void  of  this  feature. 
In  addition  we  see  the  reflected  field  is  in  the  specular  direction.  Finally,  we 
see  that  the  amplitude  of  the  forward  diffracted  field  increases  with 
increasing  frequency. 


He.  4.  Angular  Disiribuiien  of  Near  Free  Field  Fam  Funciion  for  an  Obj^  Being 
Insonificd  Along  the  Principle  Axis  of  Symmetry  for  KU2  Values  of  10, 20,  and  40. 


90 


Frequency  •  95.5  (Hz) 
KL/2-10.0 


Fig.  5.  Angular  Distribution  of  Near  Free  Field  Form  Function  for  an  Object  Being 
insonified  45  Degrees  Relative  to  the  Principle  Axis  of  Symmetry  for  KU2  Values 
of  10, 20,  and  40, 
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Figure  6  illustrates  the  case  when  the  object  is  being  insonified  at  an 
angle  of  90  degrees  from  the  principal  axis  of  symmetry.  We  note  immediately 
that  the  field  is  not  being  focused  as  we  increase  in  frequency  but  the 
amplitude  is  increasing  in  the  forward  direction.  We  also  note  that  the 
backscattered  field  remains  approximately  constant  for  the  three  cases  shown. 

To  observe  the  effect  that  range  has  on  the  field,  the  surface  on  which 
the  field  was  solved  was  enlarged  from  a  semi-major  axis  length  of  I25(m)  to 
250(m)  and  a  semi-minor  axis  length  of  25(m)  to  50(m).  The  intermediate 
frequency  of  191  (Hz)  was  selected  for  this  evaluation.  Figure  7  illustrates  the 
results  for  the  near  field.  Three  incident  angles  were  investigated  0,  45,  and  90 
degrees  relative  to  the  principle  axis  of  symmetry.  By  the  time  the  surface  is 
at  a  range  of  10  times  the  object's  dimensions  (which  is  for  this  case  25x5)  the 
field  is  approximating  the  far  field  result  rather  well  for  the  0  and  45  degree 
incident  case.  Although  this  is  not  the  case  at  90  degrees,  the  trend  is  correct; 
that  is  the  near  field  at  the  outer  surface  is  approaching  the  far  field  result, 
but  not  as  rapidly  as  for  the  cases  when  the  object  is  insonified  at  an  angle  of  0 
and  45  degrees. 

We  have  shown  the  angular  distribution  plots  because  they  will  be  used 
in  later  chapters.  There  are  other  ways  of  showing  free  space  scattering 
results  such  as  showing  the  forward  or  backscattered  form  function  vs  KL/2. 

In  this  way  one  can  sec  resonance  locations  (in  the  case  of  elastic  spheroids)  or 
Franz  waves  (38)  which  occur  when  the  specular  and  circumferential  waves 
interact  with  one  another  at  the  field  point.  To  include  these  would  be  outside 
the  scope  of  this  dissertation  since  they  have  no  bearing  on  the  problem.  We 
will  now  turn  our  attention  to  acoustic  propagation  in  a  waveguide. 


90 


Frequency  •  95^  (Hz) 
KL/2-10.0 


Frequency  - 191.0  (Hz) 
KL/2*20.0 


0 


180 


0 


270 


Frequency  •  382.0  (Hz) 
KL/2>40.0 


2  6.  Angular DismbuUonofNearFreeHeldFonnFunciionforanOb^lBeing 

unified  ^  Degrees  Rcladvc  to  the  Principle  Axis  of  Symmeuy  for  KU2  Values 
10, 20,  and  40. 


90 


Frequency  - 191.0  (Hz) 
KL/2-20.0 


0 


90 


Frequency  - 191.0  (Hz) 
KL/2«20.0 


0 


Fig.  7.  Angular  Distribution  of  Near  Free  Field  Form  Function  for  an  Object  Being 
Insonified  0, 45,  and  90  Degrees  Relative  to  the  Principle  Axis  of  Symmeuy  for  a 
KU2  Value  of  20. 
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Chapter  3 

ACOUSTIC  PROPAGATION  IN  A  WAVEGUIDE 


In  this  chapter  we  will  describe  how  sound  propagates  through  a  shallow 
water  waveguide.  We  assume  that  the  sound  is  a  small  amplitude  wavelike 
disturbance  on  the  hydrodynamic  background.  We  first  describe  the 
properties  of  the  fluid  equations.  We  will  derive  these  equations  from  first 
principles,  and  then  obtain  the  equations  governing  the  sound  field  which  we 
discuss  in  more  detail  from  two  points  of  view,  using  Normal  Mode  theory  and 
Ray  theory.  We  will  finish  by  showing  the  connection  between  the  two.  This 
section  follows  the  procedure  outlined  by  Ahlowalia  and  Keller.  (1) 


!  a. 


Lets  begin  by  considering  a  large  fluid  mass,  with  mass  density  of  p. 
Within  the  fluid  mass  lets  define  an  infinitesimal  cube  whose  sides  are  of 
length  dx,  dy,  dz.  The  center  of  the  cube  is  labeled  c.  The  pressure  at  c  is 
denoted  by  P.  Figure  8  illustrates  the  cube. 


X 
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FORCE 
FACE  b 


Fig.  8.  Geometry  for  the  fluid  cube. 


The  force  at  c  due  to  hydrostatic  pressure  on  face  a  of  the  cube  is, 

(3-1)  facea=[p-|^^]dydz 

and  for  face  b 

(3-2)  faceb=^P  +  ^Y)^y 

the  net  force  in  the  x  direction  due  to  the  pressure  on  these  two  faces  is, 

(3-3)  +  dz  =  --^dxdydz 

\  dx  2  J  \  dx  2  J  dx 

By  analogy  the  same  procedure  can  be  done  for  the  net  force  in  the  y 
and  z  direction. 


net  force  in  y  direction  =“ 


iE 

ay 


dxdy dz 


(3-4) 
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(3-5) 


3P 

net  force  in  z  direction  =— ^dxdydz 

dz 


The  resulting  net  force  is  then, 


(3-6) 


=  -VPdxdydz 


now  taking  the  limit  as  dx,  dy  and  dz  approach  zero,  gives  for  the  net  force 


(3-7) 


Equation  (3-7)  represents  the  net  force  per  unit  volume  at  the  point  c.  This  is 
just  the  internal  stress  at  point  c  due  to  hydrostatic  pressure. 

The  above  description  was  for  a  fluid  without  an  acoustic  source.  At  this 
time  wc  inaoduce  an  acoustic  source  in  the  fluid  medium.  This  is  equivalent  to 
having  an  external  force  present.  We  denote  the  external  force  by  efa  where  e 
is  a  measure  of  the  source  strength.  We  also  introduce  the  force  per  unit 
volume  due  to  gravity,  pg.  Now  the  total  force  per  unit  volume  at  point  c  is 


equal  to, 
(3-8) 


Fn=  -VP+pg+efa 


We  define  u  to  be  the  velocity  at  point  c.  Then  the  acceleration  is. 


(3-9) 


du  _  StL  ^  8u  dx  ^  ^  8u  dz 

dt  dt  dx  dt  dy  dt  dz  dt 

=  !“+(„,  V)u 

ot 


26 


Using  Newton’s  2nd  law,  ma=m^=SF.  .  we  can  write  down  an  expression 
that  equates  the  sum  of  the  forces  to  the  mass  density  times  du/dt. 

pj^|L+(u  «V)uj=  _VP+pg+  efa 

(3-10)  [If +  (“  g  +  -^ 

This  equation  is  called  the  equation  of  momentum. 

Derivation  of  the  Continuity  Equation 

To  derive  the  continuity  equation,  we  need  to  define  a  volume  of  liquid 
(V),  surrounded  by  a  surface  O.  Gauss’s  Theorem  (2)  states  that  the  rate  of  flow 

of  the  fluid  out  of  the  volume  equals  the  flow  of  the  volume  through  the 
surface  O. 

(3-11)  J^(V*A)dv  =Jj,A -do 

where  A  is  a  vector  field  representing  the  flow  of  the  fluid.  Now  the  net  flow 
of  the  fluid  out  of  the  volume  through  the  surface  results  in  a  reduction  of  the 
density  inside  the  volume. 


I^pu  •  d©  = 


(3-12) 


Now  substituting  Eq.  (3-11)  in  Eq,  (3-12), 
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(3-13) 


J^V.(pu)dv  =  -| 


Eq.  (3-13)  is  the  continuity  of  mass  equation. 


Derivation  of  the  Continuity  Equation  for  Entropy 


The  procedure  is  the  same  as  the  last  derivation.  We  define  S  to  be  the  entropy 


density. 


•^^:^dv  =  -|  (Spu)  •  do 


(3-14) 


=  -  jV  •  (Spu)d  V 


which  through  the  limiting  process  becomes, 


the  second  term  is  zero  and  one  is  left  with, 


(3-15) 


^  +  (VS) .  u  =  0 


which  is  the  continuity  of  entropy  equation. 


Equation  of  State 
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Pressure  can  be  considered  a  function  of  the  fundamental  quantities 
density  and  entropy.  ( 1 ) 

(3-16)  P  =  P(P.S) 

which  is  just  the  equation  of  state. 

Bounded  Fluid 

Now  the  fluid  will  be  bounded  above  by  a  free  surface  defined  by 
z=Ti(x.y,t)  and  below  by  a  rigid  surface  by  z=-h(x,y).  The  pressure  above  the 
free  surface  is  constant  and  is  equal  to  Pq.  There  are  two  boundary  conditions 
on  the  top  surface,  the  continuity  of  the  pressure  across  the  surface  and  the 
continuity  of  the  normal  component  of  velocity  across  the  surface. 


Pq=  constant 


Fig.  9.  Waveguide  Geometry. 


Using  the  first  boundary  condition  we  can  equate  the  pressure  on  the  opposite 
side  of  the  boundary  z=r\. 

(3-17)  P(x.y,ri(x,y.t).t)  =  Pq  when  z=Ti(x.y.t) 

To  find  the  normal  component  of  velocity,  we  need  to  construct  the  unit 
normal  vector.  Let  O  be  a  scalar  function  defining  the  surface. 


(3-18) 


<J)(x,y,z,t)=z-Ti(x,y,t) 


The  unit  vector  normal  to  <I>  is  defined  as 


Vit> 

|V<1)|  D  3x  D  3y  D 


V(t>~  -  X 


in 

dx 


+ 


z 


and 


in 

ay 


(3-19)  where 


+  1  =D 


Evaluating  the  unit  vectors  x,  y,  and  z.  we  obtain 
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•  X 


c  0  s(a)  =  -2l2!L  •  X  •  X  T  -I"  •  X 
D  9x  D  3x  D 


1  ^ 

Dax 


•y  =COS(p)  =_|-|l.y 
D  dx 


Liii 

Ddx 


y  + 


D 


•  y 


lin 

Ddy 


(3-20) 


Vd) 


•  z  s  c  0  ^  -y) 


Ddx 


«  z 


xii 

Dax 


.Z+-.Z 


D 


Now  equating  the  normal  component  of  velocity  across  the  upper  surface, 


(3-21) 


1  dxf  3^1  V  1  dyf  inV  Idz 

Ddtl  axj  Ddt  l,"ay  J  Ddt 


Since  the  derivatives  of  x  and  y  with  respect  to  t  arc  zero,  Eq.  (3-20)  can  be 
rearranged  into  the  following. 


(3-22) 


dz 

dt 


w  when  z  =  T|(x,y,t) 


This  is  the  requirement  on  the  velocity  components  at  the  top  surface.  For  the 
bottom  surface  let 


(3-23) 


d)(x,y,z,t)=z+h(x,y) 
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For  a  rigid  bottom,  the  normal  component  of  the  velocity  vanishes.  If  we 
apply  the  same  procedures  as  were  used  for  the  top  surface,  the  results  for  the 
bottom  surface  are. 


(3-24) 


,  V  1  ah 

D-sr 


cos(r)  =  5 


where 


Ds 


(3-25) 


u-^  +v-^+  w  =  0 

ax  ay 


when  z  =  -h(x,y) 


The  last  expression  is  the  requirement  on  the  velocity  components  at  the 
bottom  surface.  When  e,  the  measure  of  the  source  strength,  is  zero,  the  Eqs. 
(3-10),  (3-13),  (3-15),  (3-16),  (3-17),  (3-22),  (3-25)  constitute  what  can  be  called 
the  basic  state  equations.  When  c  =  0  this  corresponds  to  a  fluid  containing  no 
source.  Because  of  this,  the  velocity  u  =  0.  The  equation  of  momentum  Eq.  (3- 
10)  becomes 

(3-26)  VP  =  -pg. 


Then  the  continuity  equation  for  density  becomes. 


(3-27) 


i.c.  p  is  a  function  of  coordinate  only,  p=p(z). 


The  continuity  equation  for  entropy  becomes, 
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(3-28) 


i 


as 

at 


=  0, 


and  the  .ntropy  is  a  function  of  the  z  coordinate  only,  S  =  S(z).  The  pressure  is 
still  a  function  of  density  and  entropy,  i.e.  P  =  P(r,S).  The  pressure  at  z  =  0  is 
equal  to  the  pressure  on  the  other  side  of  the  top  surface,  Pq. 


Acoustic  Quantities  anj. Equations 


i 


i 


We  will  now  derive  what  Ahlowalia  and  Keller  (1)  call  the  acoustic 
quantities.  The  approach  will  be  to  consider  the  equations  which  depend  on 
the  source  term  e.  By  differentiating  these  equations  with  respect  to  e  and 

evaluating  t.  ;  results  at  e  =  0,  we  obtain  the  acoustic  equations.  The  acoustic 
quantities  are  the  individual  functions,  i.e.  u,  P  etc.  differentiated  with  respect 
to  e  and  evaluated  at  e  =  0.  Once  the  acoustic  equations  are  obtained,  the  wave 

equation  for  the  acoustic  quantity  P  will  be  obtained  along  with  the  boundary 
conditions. 


i 


t 


! 


Thus 


(3-29) 


(3-30) 


(3-31) 


(3-32) 


(3-33) 


dt 


e  =0 


f>  =  iP| 

dt 


e  =  0 


P= 


Bt 


£=  0 


<:-  9S 
Bt 


e  =  0 


dt 


E  =  0 
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Starting  with  Eq.  (3-10),  the  continuity  of  momentum  equation,  and  evaluating 
it  at  e  =  0  yields 


(3-34) 
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We  stated  earlier  that  the  velocity  u  =  0  when  e  =  0.  However  that  is  not 
necessarily  the  case  for  ; 


(3-35) 


nor  is  it  for 


ah. 

at  ■ 


(3-36) 


0 


2 

au  a  u 


at  aeat  I 

E 


*0 

0 


Therefore  Eq.  (3-34)  becomes, 


(3-37) 


Now  consider  the  continuity 


=  JLvf>  +  i-vp  +  . 

at  P  ^  p2  ^  P  • 


of  mass  equation,  Eq.  (3-13). 


(3-38) 


=  0 


the  first  term  is  zero  because  is  equal  to  zero  when  evaluated  at  e  =  0. 

ae 

Therefore  Eq.  (3-38)  becomes. 


(3-39) 


A  ae  yJ  at 


Now  consider  the  continuity  equation  for  entropy. 
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(3-40) 


V—  •  u 

+  [vs .  -^1 

L  ae  J 

L  aeJ 

the  second  term  is  zero  because  —  is  equal  to  zero  when  evaluated  ax  e 

ae 

0.  Therefore  Eq.  (3-40)  becomes, 


(3-41) 


Now  consider  Eq.  (3-16),  the  equation  of  state. 


P(ftS)  = 


9P(ftS) 


o 

II 

1 — 1 — 1 

3  cj 

• 

CO 

+ 

1  of  state. 

=  iPiP 

+  iPJS 

3p  9e 

*  A  • 

as  de 

(3-42) 


ap  ^  as 


Now  consider  Eq.  (3-17), 


(3-43) 


P(x,y.v)-f|^  -0 

3x  d£  ^  dy  de  3n  de  dl  de 

£=0  E=0  £=0  e=0 


when  z  =  1]  =  0,  we  can  define  f*  at  z=0  as, 


(3-44) 


aP(x,y,q,t) 

- 

^iLiL 

ap  at 

de 

dx  de 

dy  de 

at  de 

now  Eq.  ('3-43)  becomes. 


(3-45)  F  +  -^fi  =  0  when  z  =  11=0. 


Now  consider  Eq.  (3-21) 


(3-46) 


^  H 

+  u-^ 

^  dv  ^ 

+  v-^ 

- 

at 

de  dx 

dx 

A 

-  «  ae  ay 

ay 

de 

e  =  0 


£  =  0 


E=  0 


E=0 


E  =  0 


note  that  u  =  v  =  w  =  0  when  e  =  0  since  u  =  0.  And  since  q=  e  fa(x,y,t). 


(3-47) 


in 

de 


=  0  md  ^ 

= 

dx 

E  =  0 

3y 

£  =  0 

=  0 


E  =  0 


therefore  Eq.  (3-46) 


becomes, 


Now  for  the  bottom  condition,  Eq.  (3-25) 
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(3-49) 


when  z  =  -h(x,y).  The  second  and  fourth  terms  are  zero,  and  Eq. 
(3-49)  becomes. 


(3-50) 


w  +  u— +  v4^  =  0,  when  z  =  -h(x,y). 
dx  dy 


Equations  (3-37),  (3-39),  (3-41),  (3-42),  (3-45),  (3-48),  and  (3-50)  are  called  the 
acoustic  equations.  For  instance,  F  acts  like  a  change  in  pressure  due  to  source 
emissions.  These  equation  must  now  be  combined  in  order  to  obtain  a  single 
equation  for  the  pressure  and  two  boundary  condition  equations.  We  start  by 
taking  the  partial  derivative  of  Eq.  (3-39)  with  respect  to  time;  this  leads  to 


(3-51) 


The  second  term  is  zero  since  the  density  is  not  a  function  of  time. 
Equation  (3-51)  becomes. 
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From  Eq.  (3-37)  we  have  an  expression  for  Using  this  expression  in 

dt 

Eq.  (3-52)  gives  after  some  simplifying, 


(3-53) 


It  is  important  to  remember  that  these  quantities  are  a  function  of  z  only. 
Therefore  the  third  term  in  Eq.  (3-53)  becomes. 


(3-54) 


^p(2)_3P^ 
,P(z)  Bz) 


,P(z) 


z(-pg) 


which,  when  substituted  into  Eq.  (3-53)  gives. 


(3-55) 


•> 

a  p 


at 


^2  ap 
-V  = 

2  ^  dz 


-V.f 


Now  taking  the  partial  derivative  with  respect  to  time  of  Eq.  (3-41)  gives. 


at 


•  vs +  6 


=  0. 


(3-56) 
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Since  the  entropy  depends  only  on  z  in  the  basic  state,  the  partial  derivative 
with  respect  to  time  of  the  entropy  is  zero.  Equation  (3-56)  becomes, 


(3-57) 


3.' 


Now  substituting  Eq.  (3-37)  for  ~  in  Eq.  (3-57)  gives, 

at 


(3-58) 


2  f 

i|.+  .vs  =  o, 

at  L  p 


Remembering  that  S  =  S(z)  and  P  =  P(z),  we  obtain  for  VS  and  for  VP , 


(3-59) 


T;c  ap 

VS  =  2-:t-  and  VP  =  Z-:^ 
dz  dz 


Equation  (3-58)  becomes  after  some  simplifying. 


(3-60) 


as  ifai>  .  las  ^aas 

.  2  pUz  p  az 

ot 


Now  taking  the  partial  derivative  of  Eq.  (3-42)  with  respect  to  time  twice  gives. 


d^P  ap  ap  a^S 

,  2  “  ap  -  2  as  2 
at  ^  at  at 


(3-61) 
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Substituting  Eq.  (3-55)  for 


— ^  ar.s  Eq.  (3-60)  for  — j  in  Eq.  (3-61)  yields, 

at  at 


(3-62) 


a>  apfT,2f,..3p  ^  .  1  apasrat^  . . 


at 


p  as  az 


Rearranging  Eq.  (3-62)  gives. 


(3-63) 


at 


2  ap 


p  as  az 


^  „ap3p 

+  g- 


ap  az 


We  can  neglect  the  last  two  terms  of  Eq.  (3-63)  since 


0,  is 

az 


0  and  .^  =  0 . 
dz 


This  results  in  the  wave  equation. 


(3-64) 


VP- 


.It 

c"  a/ 


=  where 


ap 


-1 


The  quantity  c  is  just  the  speed  of  sound  in  the  fluid.  We  will  now  obtain  the 
equation  governing  the  condition  for  the  top  boundary  condition.  Starting 
with  Eq.  (3-45) 


P 


=  0 


(3-65) 


Equation  (3-65)  can  be  rewritten  as 
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<3-66) 


f’  +  Tl(-pg)  =  0 


but  fi  <  <  f*  when  z  =  0.  We  can  therefore  neglect  this  term  and  Eq.  (3-66) 
becomes, 


(3-67) 


f*  =  0  when  z  =  0. 


Consider  now  the  bottom  condition.  We  start  by  taking  the  partial  derivative  of 
Eq.  (3-50)  with  respect  to  time.  This  yields, 


(3-68) 


2  .  2 

3w  3u  3h  .  3  h  3v  3h  .  3  h  « 

3t  3t  3x  3t3x  3t  3x  3t3y  "  ’ 


since  h  is  a  function  of  x  and  y  only,  the  partial  derivative  of  h  with  respect  to 
time  is  zero.  Equation  (3-68)  becomes. 


(3-69) 


3w  ^  ^  3h,  ^  3v  3h  _  q 
3t  3t  3x  3t  3x 


We  will  use  expressions  for 


4^  and  from  Eqs.  (3-37). 
ot  ot  dt 


Equation  (3-69) 


becomes  after  the  substitution. 
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Since  the  bottom  is  defined  as  being  rigid,  no  force  is  transmitted  through  the 
bottom,  which  says  that  fl  =  fi  =  f3  =0  when  z  =  -h(x.y).  Finally  since p»0  we 

can  write  Eq.  (3-70)  as. 


(3-71) 


^iJL+ 9£.iiL  =  o 

3z  3x  3x  3y  3y 


when  z  =  -h(x,y). 


This  is  the  bottom  boundary  condition.  We  have  now  defined  a  boundary  value 
problem  for  the  determination  of  the  acoustic  pressure  using  the  wave 
equation,  Eq.  (3-64)  which  governs  the  problem,  subject  to  the  boundary 
condition  at  the  surface,  Eq.  (3-67)  and  at  the  bottom,  Eq.  (3-71).  In  order  to 
solve  the  problem,  we  need  to  know  how  the  sound  speed  varies  with  depth,  the 
depth  of  the  waveguide  and  the  source  distribution.  In  underwater  acoustics  it 
is  common  to  consider  the  source  to  be  a  point  source  producing  a 
monochromatic  time  dependent  wave.  In  this  case,  we  now  write  the  pressure 
as  a  function  of  distance  (r)  and  time  (t), 

(3-72)  f>(r,t)  =  e“‘“‘p(r). 


Substituting  Eq,  (3-72)  into  the  wave  equation,  Eq.  (3-64)  yields. 
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Defining  k  e  ^  and  n  (z)  e  where  Cq  is  a  reference  sound  speed  and  n(z) 
is  the  index  of  refraction  we  have. 


(3-74) 


and  Eq.  (3-73)  becomes, 


(3-75) 


V^P(r )  +  k ^n  ^(z)P(r )]=  V  .  f , 


Suppressing  the  time  dependence  of  Eq.  (3-75)  and  defining  q(r)  as  that 
quantity  in  parenthesis,  wc  obtain 


(3-76) 


q(r )  s  V^P(r)  +  k^n  ^(z)P(r ) 


Equation  (3-75)  can  be  expressed  as. 


(3-77) 


q(r)  =  V.f, 


Equation  (3-77)  is  sometimes  referred  to  as  the  reduced  wave  equation  or  as  the 
Helmholtz  equation.  (2)  The  boundary  conditions  arc  as  follows, 


(3-78) 


f*(r,  t)  =0  =>  P  =0  when  z  =  0  and 


f  =  who.  z  =  .h(x,y). 

dz  dx  dx  dy  dy  dz  dx  dy 


(3-79) 
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The  pressure  due  to  Eq.  (3*77)  subject  to  the  boundary  conditions  Eq,  (3-78)  and 
Eq.  (3-79)  do  not  determine  P  uniquely,  since  the  homogeneous  case  (q(r)=0) 
allows  solutions  which  represent  waves  coming  in  from  infinity.  We  will  need 
to  impose  an  additional  condition  in  order  to  uniquely  determine  P.  The 
additional  condition  is  that  the  wave  number,  k  will  become  complex, 

(3-80)  k-»-^+ia!’ 

This  condition  makes  the  outgoing  wave  decay  to  zero  at  infinity  while  the 
incoming  wave  will  be  infinite  at  infinity.  The  solution  will  be  bounded  at 
infinity  and  thus  eliminate  the  incoming  wave.  The  pressure  P(r)  now  takes 

the  following  form, 

(3-81)  P(r)s  lim  P(r  ,a). 

a  -+0 

We  stated  earlier  that  the  source  was  assumed  to  be  a  point  source,  A  point 

source  can  be  represented  by  a  delta  function.  The  source  distribution  term, 

V  •  f  g  can  be  represented  by, 

(3-82)  = 

where  Tq  is  the  source  position.  The  point  source  will  be  normalized  to 
produce  unit  strength.  Assuming  cylindrical  coordinates, 

r-ro)dr  =  l. 


(3-83) 


Evaluating  the  integral  in  Eq.  (3-83), 


45 


(3-84) 


J  8^r  -  r^^dr  =  J  f(r)5(r)^z  -z^^r  dOdrdz 
=  2kt 


therefore 


Y-ur- 


8(r)S(x-Zj) 


2jrr 


..and  {(r)=^ 


The  Helmholtz  equation  becomes. 


(3-85) 


V^P(r)  +  k^n^(z)P(r)  = 


2nr 


Recall  that  P(r,t)=  P(r)  e"‘®^ 

We  will  now  consider  a  simple  waveguide,  namely  the  homogeneous 
waveguide  of  constant  depth.  We  will  solve  for  the  pressure  using  normal 
mode  theory  and  ray  theory.  We  will  show  the  connection  between  the  two 
methods.  The  normal  mode  solution  is  most  useful  at  distances  where  only  a 
finite  number  of  propagating  modes  arc  present  and  the  ray  representation  is 
most  useful  near  the  source  where  only  the  incident  field  and  a  few  reflected 
waves  need  to  be  considered  since  spherical  spreading  diminishes  the  rest. 
Since  both  representations  should  give  the  same  answer  for  the  pressure  (at 
least  mathematically)  at  any  range,  we  will  convert  the  ray  representation 
into  the  normal  mode  representation  in  order  to  shov  the  connection. 


i 


Normal  Mode  Theory 
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Normal  Mode  theory  has  been  used  to  solve  a  variety  of  classical 
phenomena.  Problems  involving  mechanical  vibration  reduce  to  that  of  the 
harmonic  oscillator  at  small  amplitudes  of  vibration  as  long  as  they  remain 
within  the  clastic  limits  of  the  material.  The  concept  is  used  in  the  context  of 
finite,  discrete  systems  of  masses  and  springs,  stretched  strings  and 
membranes,  etc.  The  motion  of  these  systems  results  in  a  number  of  so-called 
normal  modes  of  vibration,  each  mode  behaving  in  many  ways  like  an 
independent  harmonic  oscillator.  (3,4)  Normal  modes  have  also  been  used  in 
bounded  continua,  such  as  the  acoustic  modes  of  a  room.  (5)  The  above 
examples  have  infinitely  many  but  discrete  cigenfrequencics.  The  concept  of 

normal  modes  can  also  be  applied  to  unbounded  continua.  An  example  is  the 
treatment  of  electromagnetic  modes  of  infinite  space  used  in  the  quantum 
theory  of  fields.  (6)  Biot  and  Tolstoy  (7)  have  generalized  the  procedure  to 
conservative,  unbounded,  mechanical  medium  of  any  type.  In  principle, 
normal  mode  theory  can  be  used  to  provide  a  unified  point  of  view  of  all  types 
of  mechanical,  electromagnetic,  and  electromechanical  waves. 

A  water  waveguide  can  be  visualized  by  starting  with  a  rectangular  box 
of  dimension  a,  b  and  h  corresponding  to  x,  y  and  z  coordinates.  Now  allow  a 
and  b  to  approach  infinity.  If  we  chose  appropriate  boundary  conditions  at  z=0 
and  at  z=h,  one  can  generate  what  is  commonly  called  the  ideal  or  perfect 
waveguide.  (8)  The  water  that  occupies  the  space  between  the  horizontal 
interfaces  at  z=0  and  z=h  has  homogeneous  properties.  The  ideal  waveguide, 
sometimes  called  the  slab  waveguide,  does  not  represent  a  'real'  world 
environment.  A  waveguide  that  better  represents  a  real  ocean  waveguide 
would  allow  the  speed  of  sound  in  the  water  to  vary  in  depth  and  range,  i.e. 


< 


t 


i 


c(2,r),  where  c  is  the  local  speed  of  sound  at  the  depth  position  z  and  the  range 
position  r.  In  addition,  the  bathymetry  would  vary  with  range  and  there 
would  be  multiple  layers  of  sediment  beneath  the  water  layer,  each  with  its 
own  sound  speed,  density,  sound  attenuation,  shear  speed  and  shear 
attenuation.  We  will  consider  two  variations  of  the  ideal  waveguide  in 
Chapters  5  and  6.  We  will,  however,  at  this  point  proceed  with  a  discussion  of 
the  ideal  waveguide. 

For  a  waveguide  of  constai.*  depth  and  sound  speed,  we  have  c(z)  =  Cq  = 
constant,  and  n(z)  =  Cq/c(z)  =  1.  The  governing  equation  is  the  wave  equation, 

Eq.  (3-85),  here  written  in  it's  homogeneous  fonn, 

(3-86)  V^P(r)  +  k^P(r)  =  0 

The  boundary  condition  at  the  surface  is  that  the  pressure  at  the 
surface  is  zero,  i.e.  it  is  a  pressure  relcaci  surface. 

(3-87)  P=0  at  z=0. 

The  bottom  boundary  condition  is  that  the  normal  derivative  of  the  pressure  is 
zero,  or  that  the  bottom  is  rigid. 


The  Laplacian  in  cylindrical  coordinates  is, (9) 
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(3-89)  v^()=i^r^]+-i 

^  drK  dr  J  2 


1  d\)  .  3^) 
r^  30^  dz^ 


For  this  waveguide  we  assume  cylindrical  symmetry  so  the  theta  term  vanishes 


=  0.  The  Laplacian  then  becomes, 


(3-90) 


Since  the  pressure  is  a  function  of  both  z  and  r  we  can  write. 


(3-91) 


P  ~  9(z)\|;(r) 


and  Eq.  (3-86)  becomes. 


(3-92) 


Equation  (3-92)  can  be  separated  into  two  parts  (possible  only 

when  c=constant  or  c=c(z)),  one  pan.  containing  terms  dependent  in  r  and  the 

second  part  containing  terms  dependent  in  z. 


(3-93) 


2  2 

1  1  dw  d  w  1  d  (D  2 

— i-  —  -f  r  -  4-  4.  ic  =0 

Vr  dr  .2  q>  ,  2 
L  dr  J  dz 
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I 


Denoting 


•—  by  Vr.  ^  Eq.  (3-93)  now  becomes. 


I 


» 


(3-94) 


^ r 


•*-Vrr 


k  a 


2 


where  on  the  RHS  is  the  separation  constant.  The  resulting  Ordinary 
Differential  Equations  are 


(3-95) 


.p^^  +  epk  =k 


2  2 
a  (p 


•  (3-96) 


I  .2  2 

+  r^r  =-•«  a  \|/  . 


I 


J 
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Solving  the  depth  dependent  equation  Eq.  (3-95)  first,  we  know  that  the 
general  solution  will  be  of  the  form, 


(3-97) 


(p(z)  =  A  s  i  n(k  z'J  1-  a^  )+  Be  os  (k  fjl  -a^  ). 


Due  to  the  surface  boundary  condition,  Eq.  (3-87),  the  constant  B  in  Eq.  (3-97) 
must  be  equal  to  zero  if  Eq.  (3-97)  is  to  be  a  solution.  Applying  the  bottom 
boundary  condition  to  Eq.  (3-97)  we  have. 


(3-98) 


^^=kVl-a^[A  co{kzVr^)]=0 

oz 


-h 


=  coJ^kVl-a^(-h))=0 
=  c  0  s(k  hVl  -  a^)  =  0 


50 

In  order  for  this  to  be  true  the  following  condition  must  be  satisfied, 


(3-99) 


Solving  for  a,  we  obtain 


(3-100) 


2  2 
k  h 


we  see  that  there  is  a  distinct,  discrete  solution  an,  for  each  n.  Equation  (3-100) 
can  be  written  as. 


(3-101) 


n  =0, 1,2,3,-  •• 


The  final  solution  for  the  depth  equation  is  of  the  form. 


(3-102) 


(Pn(z)  =  A„ 


n  =  0, 1,  2, 3,  •  •  • 


Now  consider  the  solution  to  Eq.  (3-96),  which  is  Bessels'  Equation  of  the  first 
kind;  solutions  are. 


(3-103) 


VW  =  CiJo(kar)  -f  C2Yo(kar). 


# 
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We  have  the  requirement  that  for  any  arbitrary  z  value,  P  is  finite  at  r=0.  This 
H*nuires  that  C2  be  equal  to  zero.  Rewriting  Jo(kar)  in  terms  of  a  sum  of 

i.%cl  Functions  (10)  of  the  first  and  second  kind  we  obtain, 

(1) 

H„  =J„+iY„ 

(3-104)  000 

(2) 

=J„-iY„ 

(3-105)  000 

(3-106)  \l/(r)  =  CHjJ\kar)  +DH^^^\kar) 

To  solve  for  the  constants  in  Eq.  (3-106),  we  note  that  for  the 
radiation  condition,  (as  r  ->  eo  ),  p  behaves  as  an  outgoing  wave.  Therefore,  in 
the  limit  as  r-^  «»,  \|r(r)  needs  to  behave  as  an  outgoing  wave.  Therefore, 

choose  the  Hankel  function  which  represents  an  outgoing  wave.  We  will  use 

the  large  argument  approximation  to  the  Hankel  function  in  order  to 
determine  which  Hankel  function  to  keep.  The  asymptotic  form  of  the  Hankel 
functions  are  (11) 


(3-107) 


(3-108) 
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Equation  (3*107)  is  an  outr  'ing  wave  while  Eq,  (3-108)  is  an  ingoing  wave. 
Since  we  want  an  outgoing  wave,  we  therefore  choose  D=0.  And  we  arc  left 
with 


(3-109) 


V(r)  =  CH 


(1) 

0 


(kar) 


Our  expression  for  the  pressure  now  takes  the  following  form, 

(3-110)  P(r)  =  (p(z)\|/(r)  =  s  i  -  a^  j  a^^  r) 

Now  the  ajj  are  positive  and  either  real  or  imaginary.  If  it  is  real  then 
we  have  propagating  modes  (trapped  modes).  If  an  is  imaginary  we  have  non¬ 
propagating  modes  (evanescent  modes).  (12,13)  That  is  to  say  that  if. 


and  if  an  =  i(l) 

&) 

Thus,  if  an  is  imaginary,  P(r)  decays 
determine  the  number  of  propagating 
Eq.  (3-100)  the  condition  that 


exponentially  as  r  increases.  We  can  also 
modes.  Since  an  is  positive  we  have  by 


(3-112) 


'"(■■■  4)  (kr) 


we  know  from  Eq.  (3-101)  that  an  is, 


(3-113) 


n  =0,1, 2,3,- --.M 


We  can  now  set  the  limit  on  the  number  of  propagating  modes,  by  the 
following  condition. 


(3-114) 

(3-115) 

It  2 

Our  expression  for  the  pressure  now  has  the  following  form 


(3-116)  P(r,  z)  =  pn  5  ‘  n 

In  order  to  solve  for  An  we  evaluate  V^P  using  Eq.  (3-90), 
Eq.  (3-117) 
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Wc  have  from  Eq,  (3*85)  the  inhomogeneous  wave  equation  that 


(3-118) 


5(z  z^)s(r)_ 

2jrr 


=  IA„sin(kz^l-a;, 
n 


L-  lA 

2  r  9r 


3r 


(1) 


H, 


With  help  from  Appendix  B  we  have. 


(3-119) 


+ 


lA 

r  ar 


4iS(r) 
2jtr  ’ 


So  that  Eq.  (3-118)  becomes 


(3-120) 


Utilizing  the  orthogonality  condition  for  the  depth  eigenfunctions 
(sec  Eq.  (3-102)),  (14)  wc  find  that  An  equals 


(3-121) 


The  expression  for  the  pressure  now  becomes, 


P(r,  z)  =  -^  Ss  i  n(kZp7^-aj]s  i  r|kz  a„r) 


(3-122) 


Rav  Theory 
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Wc  begin  again  with  the  time  independent  wave  equation  (Helmholtz) 
Eq,  (3-85). 


(3-123) 


V^P(r)  +  k^n^(z)P(r)  = 


•  2nT 


the  boundary  conditions  arc  the  same  as  for  the  normal  mode  case  namely  that 
the  pressure  is  zero  ai  the  surface  and  the  normal  derivative  of  the  pressure  is 
equal  to  zero  when  evaluated  at  the  bottom.  While  we  assumed  cylindrical 
symmetry  for  the  normal  mode  case,  we  will  assume  spherical  symmetry  for 
the  ray  representation.  The  range  variable  R  is  defined  as, 


(3-124) 


The  Laplacian  for  spherical  coordinates  is,  (15) 


(3-125) 


Due  to  the  spherical  symmetry,  the  partial  derivatives  with  respect  to  6  and  V 
are  zero.  Therefore  Eq.  (3-125)  becomes, 


(3-126) 


and 


(3-127) 


We  will  use  the  homogeneous  wave  equation  in  order  to  obtain  a  general 
solution  and  then  use  the  inhomogeneous  wave  equation  to  obtain  the 
particular  solution.  Starting  with  the  homogeneous  Helmholtz  equation, 
which  is  Eq.  (3-85)  with  n(z)  =  1  and  the  source  distribution  term  equal  to  zero 

(3-128)  V^P  +  k^P=0. 

2 

Substituting  in  Eq.  (3-127)  for  ^  P  we  obtain. 


The  general  solution  to  Eq,  (3-129)  is,  (16) 


±ilcR 

P=  “ 


(3-130) 


R 

+  ikR  -  ikR 


We  will  use  the  radiation  condition  for  spherical  waves  in  order  to  decide 
which  term  in  Eq.  (3-130)  to  keep.  The  radiation  condition  is,  (17) 


(3-131) 


lim  r[p'  (R)-ikP.(R)l=0 

R-»co  lO  0  J  , 


+  ikR  ikR  ikR 

When  P„(R)  =  A-2 - ,  P'  (R)  =  ikA-^ - A- -  and  Eq.  (3-130)  becomes, 

0  R  0  R  2 

R 


(3-132) 


1  i  m 
R  ^  oo 


Ae 


ikR 


R 


-ikR  -ikR  -ikR 

When  P  (R)  =  B^-— ,  P' (R)  =  -ikB-^-- - B^-—  and  Eq.  (3-130)  becomes, 

U  R  U  R  2 

R 


i 


therefore  B  =  0  and  Eq.  (3*130)  becomes. 


(3-134) 


ikR 

P„(R)=A^ 


In  order  to  determine  the  constant  A,  we  will  use  the  inhomogeneous  wave 
Equation,  Eq.  (3-85).  We  integrate  Eq.  (3-85)  over  a  volume  v  and  then  take  the 
limit  as  R  approaches  zero: 


(3-135) 


liml  V«(VP)dv  +  k  LPdv=-  6/z- z„'\^^dv 

R_»o[J^  Jv  ^  O^'^nr 


Evaluating  the  first  term  with  help  from  Gauss's  theorem,  (18) 


(3-136)  /X?.(VP)dv  =  4VP*d0 


where  VP  = 


becomes, 


...  ikR  .  ikR  2 

- -  ---  and  da  =  RR  sin(e)deJ9.  Equation  (3-136) 

^  R  . 


(3-137) 


V  •  (VP)dv  =  /A[ikRc‘*^^-  e**^^]s  i  n(0)d6d(p 


=  4r,Ac‘*^^[ikR-  l]. 


i 
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Now  evaluating  the  second  term  in  Eq.  (3-135), 


k^LPdv  =k^A 


^  .  n  2ji 

J  Re’^^dR  Is  i  n(0)d8  Jdcp 
0  0  0 


(3-138) 


=  -4  7rA[c ’^^(ikR-l)+  l] 


and  finally  evaluating  the  last  term, 


(3-139) 


V  0  ® 

=  -l 


Equation  (3-135)  becomes. 


1  i  m  {47tAe’^^(ikR  - 1)  -  4TCA[e'^^(ikR-  l)+  l]  =  - 1}  , 
R->0 


(3-140)  lim{4jrA=l}, 

R-»0 


therefore  A  =  ant'  Eq.  (3-130)  becomes. 


P 


0 


(R)  = 


ikR 
e _ 

4nR  ■ 


(3-141) 
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So  far  wc  have  evaluated  only  the  direct  ray.  We  will  now  evaluate  the  ray 
which  suffers  a  single  interaction  with  the  surface  and  the  bottom.  Any  ray 
which  reflects  from  the  surface  or  bottom  obeys  the  Law  of  Reflection,  (19) 
which  says  that  the  angle  of  incidence  is  equal  to  the  reflected  angle. 

Following  the  same  procedure  as  for  the  direct  ray,  the  phase  and  amplitude 
can  be  found  for  each  ray.  The  ray  that  reflects  once  from  the  surface 
appears  to  come  from  a  source  located  at  r=0,  z=-zo  (see  Fig.  (10)).  The  pressure 
due  to  this  reflected  ray  is  equal  to 


(3-142) 


P 


Sref 


(r»z) 


ikR' 

-  where 

47tR' 


R  = 


and  the  pressure  due  to  the  direct  ray  is  given  by  Eq.  (3-141).  The  total 
pressure  then  at  the  point  (r,z)  due  to  the  direct  and  surface  reflected  ray  is 
just  the  sum  of  the  individual  pressures. 


(3-143) 


P(r,z)  = 


ikR  ikR' 

£ _ .  T>  £ _ 

47tR 


where  Rg  is  the  surface  reflection  coefficient.  The  boundary  condition  for  the 
top  surface  states  tliat  the  pressure  is  zero  at  z=0.  Using  this  condition  wc  can 


determine  the  surface  reflection  coefficient  Rg. 
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(3-144) 


P(r,7.)|  =0=-^ 

Z  =0  A. 


i  k +  z , 


(3-145) 


whence  Rs  =  -1 


Equation  (3-143)  now  becomes, 


(3-146) 


P(r,z)  = 


,  ikR  ikR' 
1  e  e 


44  R  R' 


To  determine  the  pressure  at  a  point  P(r,z)  due  to  the  direct  ray  and  ^  ray 
which  interacts  with  the  bottom  once,  we  proceed  as  before.  The  ray  which 
interacts  once  with  the  bottom  appears  to  come  from  the  point  r=0,  z=(2h+zo) 
(see  Fig.  (10)).  The  pressure  due  to  this  bottom  reflected  ray  is. 


(3-147) 


where  R*=  /r^ +^z-2h  +  z^^ 


The  total  pressure  at  the  point  (r,z)  due  to  the  direct  and  the  bottom  reflected 
ray  is. 


ikR  ikR 


(3-148) 
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where  Rb  is  the  bottom  surface  reflection  coefficient.  Utilizing  the  bottom 
boundary  condition  we  can  determine  the  bottom  surface  reflection 
coefficient  (Rb)- 


it  follows  that  Rb  =  1-  Equation  (3-148)  becomes 


(3-150) 


P(r,z)  = 


J_ 

4n 


ikR 

e 

R 


#  This  procedure  can  be  generalized  for  rays  that  incur  multiple  reflections. 

Each  of  these  multiply  reflected  rays  appears  to  come  from  an  image  point. 
These  image  points  are  located  at,  z  =  ±z^  -H2n  h,  n  =  0,  ±  1,  ±  2,  •  •  • 

#  The  pressure  then  is  a  summation  of  all  these  interactions, 
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There  is  not  a  one  to  one  correspondence  between  modes  and  rays.  In  the 
homogeneous  waveguide,  all  images  are  required  to  form  one  mode  and  vice 
versa.  However,  since  both  theories  give  the  same  answer  in  the  waveguide, 
they  must  be  equal.  We  should  therefore  be  able  to  convert  from  one 
representation  to  the  other.  With  this  in  mind  we  will  show  that  the  Normal 
Mode  solution  and  the  Ray  or  Multiple  Reflection  method  from  the  last  section 
are  mathematically  connected.  We  begin  with  Eq.  (3-151),  the  final  form  for 
the  pressure  in  the  Ray  representation. 


which  can  be  written  as, 

(3-153)  P(r,  z).-L^  Zj- 1)"  {P(r.  z-  -P(r.  z  J 


writing  (-1)"  =  we  write. 


(3-154) 


where 


and 


using  the 


(3-155) 


where 


setting  a 


(3-156) 


where 


P(r,z)=P(r,  Y-)-P(r,y^.) 


n  ^ 

+  «.  ikv*’  +(Y--2nh)  -inn 


J  r^  +  (y_- 


2n  h) 


+  00  ‘H/r*  +(Y+-2nh)^  -inn 


2n  h) 


Poisson  Sum  Formula  (20) 


n  = 


m  =-oo 


>  / 


+  *o 


F(q)  =  — -=5=-J  e  ****** f(to)d to 


=  2jt,  we  obtain 


^If(2jtn)=-^  XF(m). 
n  =-oo  V  Itc  m  =  -» 


F(m)=  ‘™‘“f((o)d(o. 


Wc  can  now  write  f(27tn)  as. 
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If  wc  set  t=y-C0h/7t,  and  dco=-7!:/h  dt  and  restrict  ourselves  to  P(r,Y.)  we  find, 


(3-160) 


+  oo 


—  eo 


If  wc  introduce  0  by  i=r  sinh(0),  and  let  q=(n+.5)7t/kh  then  the  integral  of 
Eq.  (3-160)  becomes, 


+  ao 


+  00 


f  iJ7T7  +  i(n  +  i)^ 


(3-161) 


J 

—  oo 


7 


-dt  = 


ikTjl+  sinh  (©)  +  ikqt 


j 

—  oo 


r^/l  +s  i 


+  s  i  n  h  (©) 


+  oo 

I  ikrVl-^  coshCsO  , 
=  J  e  d0 


where  ©'  =  0+ 1  a  n  h  \n ) 


Now  set  sinh(0')  =s/r,  so  that  cosh(0')d0'  =  ds/r;  then 


(3-162)  cosH(eO  =  :^/l^  d©'= 


Vr^7’ 


and  Eq.(3-161)  becomes, 


(3-163) 


,  +“ 

+  00  _  ^ 

I  ikryi  -  q*  cosl/sO  ,  a 

J e  d©  =  — 


|*_ikyi  Jt^  +s^ 

TTTTy 


-ds 
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r  cosh(©)d© 


from  Magnus  et  ai.,  (21) 


(3-164) 


(1)  i 

J  Vx  +t 


Eq.(3-163)  then  becomes. 


(3-165) 


ivjl-q^  +s^ 

(r^s^) 


2  2 


i  k  r  /l  -  (n  +  f-]  -f 


Finally,  Eq,(3-160)  becomes. 


(3-166) 


From  Eq.(3-101),  —  =fn  +  — where  we  have  defined  Yn  by 
k  V  2  Ah 


(3-167) 


and  now  Eq.(3-166)  becomes. 


(3-168) 


with 

n 


P(r,y-)  = 


le 

n 


(1  ) 

»0  M 


and  Y-=z-Zo 


Finally  Eq.(3-168)  and  the  equivalent  equation  for  z+zq  can  be  written 


Subtracting  Eq.  (3-169)  from  Eq.  (3-170)  (see  Eq.  3-153),  one  obtains  an 
expression  for  P(r,z), 


(3-171)  P(r 


,z)=— He  -e  JH  (Kj^r) 


8h 


1 

8h 


y  -'Vf  ‘Vo  -■‘Vo’lO*  , 

&  U  -e  )H  (K„r; 


8h 


Ic 

n 


(1) 


(2isiii(Y„Zo))H  (K^r) 


(1) 


Since  the  sum  is  from  n=-  "<>  to  n=+  «»,  the  first  term  in  the  RHS  vanishes,  and 
the  second  term  equals  4  s  i  n^Yj^z^s  i  n^7  ^^z^) .  The  result  is 


(1) 


(3-172)  >’(^^)-2r?‘”(l’n^)'‘"(l’n^o)”o  C'n') 


which  is  Eq.  (3-122),  the  normal  mode  representation  for  the  acoustic 


pressure. 
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CHAPTER  4 

METHODS  DESCRTOING  SCATTERING  FROM  OBJECTS  IN  A  WAVEGUIDE 


In  this  chapter  we  begin  by  reviewing  a  few  of  the  previous  methods 
used  to  describe  acoustical  scattering  from  objects  in  a  waveguide.  While  the 
methods  given  are  not  exhaustive,  they  represent  the  major  work  in  this  area. 
We  then  formulate  the  present  method  based  on  Huygens'  principle.  We 
conclude  this  chapter  by  obtaining  a  general  expression  for  the  scattered  field 
far  from  the  object. 


Ptiviouslv  Used  Methods 


The  problem  of  describing  the  scattering  of  acoustic,  electromagnetic, 
etc.  wave  energy  from  an  object  in  free  space  is  difficult.  One  has  to  contend 
with  numerical  instabilities,  ill  posed  problems  etc.  When  one  allows  the 

object  to  be  in  a  waveguide,  the  problem  greatly  increases  in  difficulty.  One  of 
the  earliest  methods  which  allowed  for  numerical  answers  was  the  sonar 
equation  approach  to  scattering  from  an  object  in  a  waveguide.  This  method 
simply  takes  the  transmission  loss  (-201og(P/Po))  from  the  source  to  the  object, 
adds  the  target  strength  due  to  a  single  plane  wave  incident  on  the  object  and 
then  adds  the  transmission  loss  from  the  object  to  the  field  point  (receiver). 
Mathematically  it  takes  the  following  form. 


U  =TL^„  +  TS  ,,.  +  TL„, 
s  so  obj 


(4-1) 
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The  problem  is  that  there  is  only  a  single  number  that  one  obtains  from  this 
process.  There  are  no  insights  as  to  what  physical  processes  arc  going  on  with 
the  scattered  field.  This  method  however  gave  surprisingly  good  results  when 
the  field  point  (receiver)  was  in  the  forward  direction,  and  when  the 
frequency  was  high  enough  so  that  the  scattered  field  was  highly  focused  in 
the  forward  direction. 

Another  approach  taken  to  solve  this  problem  was  to  use  ray  theory  and 
follow  individual  rays  as  they  interacted  with  the  boundaries  and  the  object. 

In  theory  this  method  seems  very  plausible;  however  it  becomes  intractable  in 
practice.  This  method  has  been  developed  by  Evans  (1)  and  Werby  and  Evans 
(2)  for  an  object  in  a  half  space  and  the  theory  has  been  worked  out  for  an 
object  in  a  waveguide,  A  similar  method  has  been  presented  by  Hackman  et  al 
(3-5)  for  elastic  objects  including  shells. 

A  more  straightforward  method  was  developed  by  Evans  (6)  which 
decomposed  the  incident  field  (composed  of  a  sum  of  normal  modes)  into  an 
equivalent  plane  wave  representation.  Each  of  these  plane  waves  are  allowed 
to  insonify  the  object  in  order  to  generate  the  resulting  scattered  field.  This 
field  is  then  projected  onto  a  cylinder  which  encloses  the  object  and  spans  the 
depth  of  the  waveguide.  This  method  is  known  as  the  projection  method  and  its 
results  are  compared  with  the  present  method  in  Appendix  E. 

A  method  based  on  matched  asymptotics  was  devised  by  Collins  (7,8)  to 
derive  an  approximate  expression  for  the  field  scattered  from  a  spheroid.  The 
matched  asymptotic  solution  consists  of  an  inner  solution  valid  near  the 
scatterer  and  an  outer  solution  valid  away  from  the  scattcrcr.  The  two  parts  of 
the  solution  are  matched  in  the  region  in  which  they  arc  both  valid.  The 
resulting  field  is  then  propagated  in  the  waveguide.  Sen  (9)  has  reported  the 
use  of  matched  asymptotics  to  describe  scattering  of  acoustic  waves  in  a 
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waveguide.  He  has  reported  that  the  method  can  be  used  to  determine  globally 


valid  pressure  field  junction  conditions  near  a  boundary  discontinuity. 

In  addition,  Klcshchcv  has  written  several  papers  dealing  with 
scattering  from  objects  in  waveguides  and/or  sound  channels.  (10-12) 


In  1690,  C.  Huygens  published  Traite'  de  la  Lumiere.  In  it  he  discussed  a 
new  principle  pertaining  to  the  propagation  of  light.  Huygens  proposed  that 
at  any  instant  t=to.  a  point  source  of  light  generates  a  disturbance  which  is 
propagated  into  the  surrounding  medium  as  an  isolated  spherical  wave, 
expanding  at  the  velocity  of  light.  The  disturbance  at  time  t=ti  is  due  to  a 
succession  of  disturbances  at  intermediate  times  between  t=to  and  t=ti.  The 
actual  effect  then  at  time  t=ti  is  the  result  of  all  the  secondary  disturbances. 

In  order  to  determine  the  effect  at  time  t=ti  caused  by  an  initial  disturbance  at 
t=to  we  then  calculate  the  state  at  some  intermediate  time  t=t’,  then  at  t=t",  etc., 
until  we  are  at  time  t=ti. 

There  are  some  limitations  imposed  by  this  theory  namely  that  the 
rectilinear  propagation  of  light  can  be  accounted  for  only  by  assuming  that  a 
secondary  wave  has  no  effect  except  at  the  point  where  it  touches  its  envelope, 
and  in  addition  it  is  assumed  that  since  a  disturbance  on  the  envelope  of  the 
secondary  waves  propagates  in  both  the  forward  and  backward  direction,  only 
the  forward  propagated  wave  was  to  be  considered. 

The  present  method  uses  the  basic  concept  of  Huygens'  principle  namely 
that  the  field  propagating  away  from  the  object  and  is  not  contributing  to  the 
field  at  the  receiver  can  be  neglected.  Huygens  reasoning  for  making  this 
basic  premise  was  based  on  empirical  observations.  We  make  this  assumption 


based  on  the  type  of  propagation  model  used.  Propagation  models  either  allow 
for  one  way  or  two  way  propagation.  If  the  model  performs  one  way 
propagation  then  the  field  always  propagates  away  from  the  source  and  no 
energy  via  backscattering  from  rough  surfaces,  rough  bottoms  or  objects  arc 
considered.  If  on  the  other  hand  the  propagation  model  allows  energy  to  be 
received  from  the  above  mentioned  mechanism  then  the  field  at  the  receiver 
can  be  dependent  on  what  environmental  or  physical  changes  happens  in  the 
waveguide  beyond  the  receiver.  For  example,  if  beyond  the  receiver  there 
happens  to  be  a  sea-mount  (an  underwater  mountain)  then  acoustical  energy 
could  be  scattered  back  toward  the  receiver  affecting  the  field  at  the  receiver. 
In  other  words,  for  a  one  way  model,  the  field  at  the  receiver  is  dependent  on 
what  happens  to  the  field  as  it  approaches  the  object.  For  a  two-way  model,  the 
field  at  the  receiver  is  dependent  on  what  happens  before  and  after  the  field 
reaches  the  receiver. 

We  will  now  formulate  the  scattering  problem.  The  geometry  used  to 
describe  the  scattering  is  shown  in  Fig.  1!.  There  are  two  regions,  region  1 
which  contains  the  scatterer,  and  region  2  where  the  field  will  be  determined, 
region  2  is  free  of  sources.  The  surface  S2  encloses  region  2.  The  surface  S2  is 
taken  to  be  at  a  large  distance  from  the  region  of  interest.  The  surface  SI  is  an 
arbitrary  mathematical  surface  enclosing  the  scatterer.  For  free  space  or 
unbounded  problems  the  surface  S2  could  be  extended  an  infinite  distance 
away  from  the  area  of  interest. 


REGION  2 


For  a  problem  dealing  with  scattering  from  an  object  in  a  waveguide,  the 
surface  S2  can  still  extend  to  infinity  but  this  necessitates  that  we  use  the 
Gre>m's  function  which  satisfies  the  waveguide  boundary  condition.  See  Fig. 
12. 


Fig.  12.  Specific  Scattering  Geometry  for  a  Waveguide. 
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To  establish  Huygens'  principle  for  the  scattering  problem  we  need  to  develop 
a  relation  between  the  known  value  of  the  field  on  the  surface  SI  and  the 
unknown  value  at  the  field  point  (receiver).  Assume  that  the  field  is  a  solution 
to  the  scalar  Helmholtz  wave  equation, 

(4-2)  (v^+k^)v(x)  =0 

and  the  Greens  function  for  the  Helmholtz  wave  equation  is  given 

by. 


(4-3) 


(v^+k^)G(x  -x')  =  8(x  -xO 


We  obtain  this  relationship  by  invoking  Green's  theorem  (13)  which  equates 
an  integral  of  a  vector  function  over  a  closed  surface  to  an  integral  of  a 
related  function  over  the  volume  enclosed  by  the  surface, 


(4-4) 


/  (v|fV^G  -  GV^v )d V'  =  J(  vVG -  GVy)  •  n '  ds' 


where  in  this  case  the  surface  S  is  SI  +  S2.  The  surface  integral  is  equal  to  zero 
when  the  observation  point  (receiver)  and  the  singularities  of  \jr  lie  outside 

the  surface  S.  (14)  If  the  field  point  is  inside  the  surface  (in  the  bounded 
volume)  then  the  scattered  field  becomes. 


(4-5)  Uj  =  y(x)  =  J(v(x')vtj(x  -  x')- G(x- x')vV(x'))  •  n'ds' 

s 


This  integral  over  S  is  divided  inio  two  parts,  over  the  surface  SI  bounding  the 
scatterer,  and  over  the  surface  at  infinity  S2.  Since  the  field  is  propagated 


from  the  surface  SI  to  S2,  it  is  composed  of  outgoing  waves.  We  are  then  left 
with  the  integral  over  SI.  The  scattered  field  then  in  region  2  becomes, 


(4-6)  Uj  =  J^(\j/(x')v'G(x-x')-G(x -x')vV(x'))*n'ds' 

We  replace  ijr  with  the  free  field  scattered  field  f  which  is  valid  at  the  surface 
SI.  We  then  have, 

(4-7)  Uj  =  J^j(f(x')vt3(x  -x')  -  G(x  -  x')v'f(x'))  •  n'ds' 

This  is  the  expression  for  the  scattered  field  far  from  the  object.  We  apply  this 
expression  to  two  different  waveguides  in  the  remaining  two  chapters. 
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CHAPTERS 

APPLICATION  OF  HUYGENS'  METHOD  TO  AN  ISO\^OCITV  WAVEGUIDE 


In  this  chapter  we  utilize  the  expression  developed  in  Chapter  4  to 
describe  scattering  from  mi  object  in  an  isovelocity  waveguide.  We  begin  by 
solving  for  the  acoustic  pressure  in  the  absence  of  an  object.  This  will  be  the 
acoustic  field  incident  upon  the  object.  We  then  show  the  development  of  the 
near  field  scattered  field  (the  field  that  is  valid  a  short  distance  from  the 
object).  We  then  develop  the  far  field  scattered  field.  We  end  this  chapter  with 
an  example  which  shows  that  the  object  acts  as  a  source  in  the  waveguide 
which  implies  that  the  object  is  now  fully  integrated  into  the  waveguide.  In 
Appendix  E,  there  is  a  comparison  of  this  method  with  the  projection  method 
described  in  Chapter  4. 

We  utilize  Normal  Mode  theory  to  determine  the  incident  field.  We 
consider  a  waveguide  consisting  of  an  isovelocity  layer  of  water  over  an 
isovelocity  half-space.  The  bathymetry  will  be  flat  and  the  half-space  will  not 
support  shear.  Sound  speed  in  the  half-space  will  be  greater  than  that  of  the 
water,  as  vill  its  density  and  attenuation  likewise  will  be  greater  than  the 
water's.  Such  a  waveguide  is  sometimes  called  the  'Pekeris'  waveguide,  (1) 
since  Pekeris  was  the  first  to  thoroughly  investigate  it.  As  mentioned  in  the 
introduction,  we  will  solve  the  acoustic  field  scattered  from  the  object  first 
using  the  Pekeris  waveguide  and  then  using  the  more  general  multilayered 
waveguide.  We  begin  by  deriving  the  incident  field  for  the  Pekeris 

waveguide.  The  following  development  is  based  on  Clay  and  Medwin.  (2) 


Jh.gpry 
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The  following  assumptions  will  be  made  when  deriving  the  incident 


field. 

1 )  The  object  is  in  the  farfield  of  the  source,  this  allows  the  use  of 
the  asymptotic  form  of  the  Hankel  Function. 

2)  There  is  a  pressure  release  air-water  interface,  this  would 
correspond  to  having  a  boundary. .  condition  such  that  the 
acoustic  pressure  evaluated  at  the  surface  is  equal  to  zero, 
i.e.P(0)=0. 

3)  The  pressure  is  continuous  across  the  water  layer  (subscript  0) 
and  the  semi-infinite  half-space  (subscript  1).  This  corresponds 
to  the  following  boundary  condition. 


P  =P  : 
0  1 


4) 


The  vertical  component  of  displacement  is  equal  across  the 
interface  of  layer  1  and  layer  2.  This  corresponds  to  the 
following  boundary  condition, 


(%)  =("i)  = 

z  z 


± 

dz 


fa.,] 

dz 

1  3t  J 

_  a 

_  a 

at 

1  az  j 

at 

;  az  > 

ra<Pi] 

1  az  J 

az  > 

5) 


The  source  produces  a  symmetric  field  about  the  z-axis. 
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We  wi  vrite  the  particle  displacement  d  and  particle  velocity  u  in  terms 


of  the  dispL  .lent  potential  O: 


(5-1) 


d=V<l>  andu=v4*. 

dt 


Newton's  Second  Law  becomes 


(5-2) 


Solving  for  the  pressure  P,  we  obtain 


(5-3) 


where  p  is  the  density  of  the  medium.  For  a  harmonic  point  source,  the 
pressure  is 


(5-4) 


P  =  (o'^ptp 


The  governing  equation  is  the  wave  equation, 


V^P  = 


a^p 


at 


(5-5) 
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where  c  (water  sound  speed)  is  time  independent.  Substituting  Eq,  (5-4)  into 
the  wave  equation,  Eq.  (5-5)  gives  the  time  independent  wave  equation  or 
Helmholtz  equation,  (3) 


2 

d  (p 


dt 


0  r 

(5-6)  ‘P=® 

where 

We  will  use  the  separation  of  variable  technique  (4)  to  solve  this 
equation.  The  Laplacian  in  cylindrical  coordinate  is  (5) 


(5-7) 


1  d 

2,2 
r  30 


Because  of  azimuthal  symmetry  the  Laplacian  becomes 


(5-8) 


3 


3r 


The  solution  to  Eq.  (5-6)  will  take  the  following  form, 


(5-9) 


(p=U(r)Z(z) 
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Using  Eq.  (5-8)  and  (5-9)  in  Eq.  (5-6)  yields  two  ordinary  differential 
equations,  in  range  and  depth.  The  range  solution  is. 


(5-10) 


ir(r)+|-U'(r)+A(r)  =  0 


where  K  is  the  separation  constant,  and  the  depth  solution  is. 
(5-11)  Z'(z)+(k^-K^)z=0. 


where  k  =  “y  is  the  characteristic  wave  number  and  Equation  (5-11)  can  be 
rewritten  as. 


(5-12) 


Z'(z)  +  Z=  0 


2^2 

where  Y  =k  -K  and  is  the  vertical  component  of  the  wave  number.  The 
wave  number  k  is  constant  in  the  waveguide,  while  the  horizontal  component 
K  and  the  vertical  component  y  will  vary  with  mode  number  (Ic=7  +  TJ). 
Consider  now  Eq.  (5-10), 


(5-10) 


U*(r)+fu'(r)+K^U(r)  =  0 


this  is  Bessels'  equation  of  zeroth  order,  whose  solution  is  the 
cylindrical  Bessel  function  of  the  first  kind  JQ(Kr)  (6)  thus 

U(r)  =  Jjj(Kr) 


(5-13) 


The  asymptotic  form  of  Jq  is  (7) 
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when  Kr  » 1 


(5-14) 


a/IjtkF 


The  first  exponential  represents  an  incoming  wave  with  infinite  amplitude  at 
infinity  and  the  second  exponential  represents  an  outgoing  wave  with  zero 
amplitude  at  infinity.  We  are  interested  only  in  the  outgoing  wave,  therefore 
the  solution  for  the  radial  term  takes  the  following  form. 


(5-15) 


U(r). 


VljtKr 


-(Kr-J) 


Now  consider  the  depth  solution,  Eq.  (5-12): 


(5-12) 


Z'(z)  +  Y^  Z=0 


the  general  solution  has  the  following  form,  (8) 


(5-16) 


Z(z)  =  As  i  n(Yz)  +B  cos(yz) 


4 


Implementing  the  surface  boundary  condition,  that  is  P=0  at  z=0: 


(5-17) 

(5-18) 

(5-19) 


P(r,z.t)=0 

P(r ,  z,  t )  =  (o^p<|>  e‘“  «  0)^pU(r  )Z(z)c'“* 
P(r,0,t)=<i)^p  U(r)Z(0)e‘“‘ 


which  implies  that  Z(0)=0,  hence  B=0  in  Eq.  (5-16).  The  depth  Eigenfunction 
then  takes  the  form. 


(5-20) 


Z(z) »  A  s  i  n(Yz) 


At  this  time  wc  will  derive  the  characteristic  equation  for  the  waveguide. 
To  accomplish  this  we  define  the  reflection  coefficient  to  be  (9) 


(5-21) 


The  reflection  coefficient  for  the  upper  layer  is  given  by 


(5-22) 


i2«l. 

Ru=-l=c 


which  implies  2<I>u  =  Jt  or<bu  =  j.  When  the  incident  angle  ©q  is  less  than  the 
critical  angle,  6^,  the  reflection  coefficient  r  the  lower  interface  is  given 
by  (10) 


PlCtC0^9,)-p,c^cos(e^) 

(5-23)  + 


4 
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-  1  I 

where  by  Snell’s  law  ©  =sin  — ^sin^0  '<  •  When  the  incident  angle  is 

1  \  0/ 

\  0  / 


larger  than  the  critical  angle  and  less  than  —  the  reflection  coefficient  for 
the  lower  interface  is  given  by 


(5-24) 


Rj  =  c  ^  when..e^<e^<f 


and  it  follows  that 


(5-25) 


/  1  2 

where  b  =  /  - —  sin  f©,")-!.  Utili^ing  the  boundary  condition 

>  JroJ 


z=h,  and  Eqs.  (5-22,  5-24)  the  characteristic  equation  becomes, 


(5-26) 


-UTh  +  , 

=e  =1 


rewriting  the  exponential  using  trigonometric  functions  gives. 


(5-27) 


+  «l>j  -  yh^^  +  i  s  i  +  <l)j  -yh  =  1 


I 


which  implies  2^0^^ +<I>j -yh^  =2jtn,n  =0,±  1,  ±2,  ••  • 


o  r 


Vnjh -<>„ =  (m  -  1)k  m=1.2.3. 


o  r 

(5-28)  Yroh  =  [m-l]«  +  0  . 
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I 


I 


We  can  rewrite  the  horizontal  and  vertical  components  of  the  wave  number  as 
the  following. 


''m  =  '‘sin(6j=isiii(ej 

and 

(5-29)  = 


Equation  (5-28)  can  be  solved  for  (O  with  the  help  of  Eq.  (5-29)  with  the  result 
that 


(5-30) 


%[(■” 

hcos(e„) 


which  can  be  rewritten  as 


I 


I 


I 


I 


I 


(5-31) 


(dh  CO 


K^ra) 


1 ''  - 1 
m  -  ^  }rt  =  <1>  -tan 

] 

V  2j 

i 


I 


This  can  be  put  in  the  following  form. 
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1 

( 

PQ 

P  0) 

% 

-^b 


(5-32) 


=  tan 


-  1 


Defining 


#  (5-33) 


-k 


2 

1 


where  k  =  .  We  now  rewrite  the  right  hand  side  of  Eq.  (5-32)  as  the 

following. 


(5-34) 


tan 


-1 


“  K 


^Pikcos(ej 


=  tan 


-1 


Pft 

u  m 


Pr^m 


so  that 


(5-35) 


b_^ 
U  m 

Pi  Ym 
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or 


(5-36) 


'‘"o.Sni' 


m  y 


The  left  hand  side  of  Eq.  (5r36)  can -  be  rewritten;  as 


/  /  \ 


(5-37) 


.  ^<=°<'‘7  J 

±siil(hTjj,)  ^m) 


After  using  Eq.  (5-37)  in  Eq.  (5-36),  the  following  equality  now  holds, 


or 

(5-38) 


-C0l(hTj 


PjY 


ta 


Equation  (5-38)  is  the  characteristic  equation  for  the  isovelocity  water  column 
overlying  an  isovelocity  half  space.  It  will  be  used  to  solve  for  the  cigenvalm 
Ym.  which  are  solutions  to  the  depth  eigenfunction  equation,  Eq.  (5*12) 
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Now  we  will  derive  the  depth  eigenfunction  for  the  second  layer.  The 
sound  pressure  is  exponentially  damped  because  the  incident  sound  Held  is 
totally  reflected.  Using  Snell's  law,  (11) 


(5-39) 


This  implies  that 


or 


,f.in 


(5-40) 


cos(ei)  = 


1- 


The  argument  in  the  radical  is  positive  if  6^  <  6^  that  is,  if  the  incident  angle 
is  less  than  the  critical  angle.  At  the  critical  angle. 


(5-41) 

and 


(5-42) 


cos^e^j  =11- 


fill 

1 

sln(eo)  =  0. 


When  6^  >  0^  .  the  argument  inside  the  radical  of  Eq.  (5-42)  is  negative: 


cos(ej)  = 


1- 


(5-43) 


Modes  will  propagate  only  in  the  upper  layer  only  for  0^  >  0^  . 
Therefore, 


k  j  CO  s^0j^  =i 


(5-44) 


Defining  by  the  following: 


(5-45) 

we  find 


P 


m 


-k 


2 

1 


(5-46) 


so  that 
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(5-47) 


=  ib 


m 


2  2  2 
where  b„=K„-k,. 
m  m  1 


The  depth  eigenfunction  for  the  second  layer  takes  the  following  form. 


->Pn)(z-h) 

Z„(z)=C„e  ^  +D„e  ® 
m  'in  in 


(5-48) 


=  C  e  “  +D„e 
m  m 


bm(z-h) 


The  flrst  term  on  the  right  damps  the  signal  as  the  wave  penetrates  into  the 
second  layer.  The  second  term  increases  the  strength  of  the  signal  as  it  gets 
farther  in  the  sediment  (and  farther  from  the  source).  We  therefore  choose 
DqisO.  Equation  (5-48)  now  becomes, 

-b  (i-h) 

(5-49)  Z  (z)  =  C  e  ®  z2!h 

m  in 


We  use  the  boundary  condition  at  the  water/fluid  bottom  interface  in  order  to 
solve  for  the  modal  coefficient  Cqj.  This  requires  that  the  pressure  across  the 

interface  be  continuous: 


(5-50) 


PoZojjj(h)*PjZi^(h)  (refer  to  page  77) 


# 
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or 


P,>;n(T 


m'')  =  Pfa' 


-bJH-h) 


(5-51) 

Therefore 


(5-52) 


Now  the  depth  eigenfunction  for  the  fluid  bottom  takes  the  following  form: 

Pa  -b_(z-h) 

(5-53) 


Before  we  can  solve  for  the  acoustic  pressure  Lt  the  waveguide,  we  need 
to  solve  the  modal  coefflcient  Am.  The  deptl;  eigenfunction  equation  for  the 
water  layer  Eq.  (5-20)  can  be  rewritten  as  Avllows: 


(5-54) 
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To  obtain  the  modal  coefficient  A01  we  use  the  orthogonality  of  depth 
eigenfunctions. 


m  ^  n 


(5-55) 


=j: 


m  =  n 


Since  we  expand  the  source  function  5(z-zq)  in  terms  of  these  orthogonal 
eigenfunctions  Zq: 


(5-56) 


6(z-zo)*  1  A„Z„(z) 

n  si 


so  that 
(5-57) 

0  r 


Po^m(*o)  * 
Po^m(*o) 

- vI - 


(5-58) 


The  sound  pressure  from  Eq.  (5*4)  is 
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(5-59)  P  =  e)2p,»  =  fl,2p 

^  01 

Now 

P(r,z,t)=a)2p  lA„Z„(z)U^(r)T(t) 

”  IH 

(5-60) 

Finally. 

(5-61) 

where  the  constant  C  has  absorbed  all  numerical  constants,  coefficients,  and 
the  time  dependent  term,  and  we  have  introduced  the  attenuation  term.  e*"'.  A 
brief  discussion  of  the  attenuation  term  is  presented  in  Appendix  D.  This  is  the 
field  in  the  absence  of  the  object.  To  obtain  the  total  Held  with  a  scattering 
object  present,  one  would  add  the  scattered  field  to  this  incident  field.  We  now 
proceed  to  insonify  the  object  using  this  incident  field. 

The  derivation  follows  Norton.  (12)  The  geometry  is  shown  in  Fig.  13. 

The  object  has  as  body  fixed  axis.  Xq.  Yq.  Zq  .  where  Zg  is  the  axis  of  symmetry 

and  is  parallel  to  the  waveguide  boundaries.  The  angle  that  the  symmetry  axis 
of  the  object  makes  with  the  vertical  plane  containing  the  source  and  the 
center  of  the  object  is  labeled  a. 


I 


I 


source 


Fig.  13.  Waveguide  Geometry  with  Object. 


The  connection  between  the  cylindrical  coordinates  of  the  source  and  the 
body-fixed  rectangular  coordinates  of  the  object  is  given  by 

(5-62)  r,  =  r„-co<a)Z„-$iii«4Z„, 

Since  we  required  that  the  object  be  in  the  far  field  of  the  source, 
rj  =  rgo  and  the  factor  .^k^^r  ^  s  >s  nearly  constant  in  range  at  the 


object.  Equation  (S-61)  can  be  rewritten  using  the  exponential  form  of  the 
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sin  (  YqZ)  term  and  Eq.  (5>62),  resulting  in  the  following  expression  for  the 
incident  field  in  the  rectangular  system. 


(5-63) 


U.(r,z)=CS 
‘  n 


sin(Yn^s^®  c 


so 


•C 


i(Yn»K|,  sii<a),KuCoa(a))#(Xo,Yo,Zo) 


-iY„z  ..  -iK„r.„ 

*  n  obj  n  so  ” « 


SO 


•  e 


i( -  Y„.  K„  s  I  Ko ) ,  K„  c  0 <a) ) .  (X^,  Y^,  Z„) 


We  have  specified  the  depth  of  the  field  point  to  be  at  the  depth  of  the 
object,  Zq^j  and  Zs  is  the  source  depth.  Note  that  there  are  two  sets  of  plane 

waves,  a  downgoing  and  upgoing  set.  They  are  expanded  in  a  spherical 
representation  so  as  to  be  compatible  with  the  T-matrix,  since  the  T-matrix  is 
in  a  spherical  representation. 

A  spherical  coordinate  system  is  established  at  the  center  of  the  object. 
I'be  angle  6  is  measured  from  the  Zq  axis,  and  <p  measures  the  rotation  in  the 
^0*  ^0  Plauc.  The  object  is  invariant  under  rotation  in  <p.  The  angles  are 
defined  as  follows, 

0=cos-^^),(p=tan-^^),  P=7^o+ 

0 


(5-64) 
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The  exponential  tenns  of  Eq.(5*62)  can  be  rewritten  as 

tR<5 

(5-65)  c 

Where  Randp  are  given  by  the  following: 

(5-66)  K  =  (±  Yjj  i  +Kjj  si  n(a)j  +  cos(a)  k) ,  p  =  (J^  i  +  j  +  ^  k) 

We  now  expand  Eq.  (5-65)  in  spherical  waves.  The  expansion  follows  Morse 
and  Feshbach.  (13) 


(5-67)  =  I  (21  +  l)iV(kr)  £  f}^^^p'”(cos(e))  p“(cos(e')) 


•  [  (c 0 <ni(^ c 0 s(in (p) )  4^  (s  i  n(m  s i n(m (p7 )  ] 


where  6  and  (p  are  relative  to  K  and  6'  and  (p'  are  relative  to  p  .  If  we  let 
l  =  e 

a=  (  ) ,  and 

U  =  0 


(5-68) 


^1  47c(l+m)! 


e„=l  6^=2  ,  m  ?t0 
0  ro 


then  Eq.  (5-67)  becomes. 


r'  .k  m^*  -V,,  [sin(mo)o  =  0l 

(5-69)  e  ^  =  4n  I  1(7“  P,  (cos©  '  \  j 

1  [cos(mq))  o=lJ 


1 

siii(m<p‘)o  =  0 

j(k  r)P^(cos(e'))^ 

► 

co^m(p')o  =  1 

Equation  (5-69)  can  be  simplified  by  making  the  following  definition: 
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(5-70) 


=  4  Jti*  Y™  (e,<(J 
i«  o 


sin(m<p)o  =  0 
co^m<p)a  =1 


r  m  fsin(m<p')o  =  0 

Re'P  (P)=fY™l  j,(kr)P,(cos(eO)- 

VW  ‘  corfm®')o=l 


(5-71) 


=j  (kr)Y  (0V) 
1  I9  c 


Where  Re  denotes  the  regular  part  of  4'.  Utilizing  Eqs.  (5-70)  and  (5-71),  Eq. 
(5-69)  takes  the  following  fonn. 


(5-72) 


e  '^=  I  a”  Re'P(p) 
l,m,o  1.0 


The  incident  field  (Eq.  (5-63))  in  terms  of  the  spherical  coordinates  is  now 
given  as 


(5-73) 


U.(!))-I  I  (P.  Re  (P) 
'  n  l,m,o 


where  n  is  the  index  for  the  number  of  modes.  The  partial  waves  coefficients 


( 


of  the  incident  field  are  then  given  by  the  following, 


(5-74) 

a™,  - - - — *(6  -c  )a 

n,l,o  2iVj^^Kjjr  *»  ® 


so 


where  Zg  is  the  depth  of  the  source,  ZqJjJ  is  the  depth  of  the  object  and  rgg  is  the 
range  from  the  source  to  the  object.  The  scattered  field  about  the  object  is 
written  as  (see  chapter  two) 


(5-75) 


f  (P)  =  I 

If  nil  G 


p”  1 

Hi  If  C  1 1  O 


where 


P 


m 

n,l,o 


o 


m 

n,r,o 


and 


(8,ip) 


The  scattered  field  given  by  Eq.  (5-75)  is  valid  only  near  the 

object  before  it  has  a  chance  to  interact  with  the  waveguide  boundaries. 

Notice  that  the  spherical  Bessel  function  used  in  Eq.(5-71)  to  describe  the 

regular  spherical  wave  (Re'P™  (p))  is  now  replaced  by  the  spherical  Hankel 

1,  o 

function  to  obtain  the  outgoing  spherical  wave  ( 'P™  (p) ).  Having  the 

1,  o 

scattered  field  about  the  object  we  can  generate  the  expression  for  the 
scattered  field  far  from  the  object. 
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The  scattered  near  field  does  not  satisfy  the  boundaiy  conditions  of  the 
waveguide.  The  solution  for  the  scattered  field  far  from  the  object  (from 
chapter  4)  is: 


(5-76) 


U3(r^,z)  = 


(f(r) 


3f(r) 


where  the  surface  a  is  a  closed  surface  enclosing  the  object.  Figure  14  shows 
the  geometry  of  the  waveguide  with  the  object  centered  at  the  origin  of  the  Xo, 
Yo,  Zo  coordinate  system. 


Fig.  14.  Geometry  of  Waveguide  from  Object  to  Receiver. 


We  now  need  to  determine  the  normal  derivative  of  the  scattered  field  and  the 
Green's  function  that  satisfies  the  waveguide  boundary  conditions,  and  its 
normal  gradient  on  some  arbitrary  surface  sufficiently  near  the  object. 


The  appropriate  Green’s  function  (GCr^.z)  )  is  given  by  Eq.(5-61), 
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(5-77) 


G(rj,z)  =  C2:- 
n 


iK  r  -  a„ r  . 
sin(Y„zpsln(Y„z)e  "  ‘e  "  ‘ 


'nVV 


where  zr  is  the  receiver  depth  and  rg  ,  to  within  first  order,  is  given  by 


(5-78) 


ts  =  r^jj-r  €08(0) 


where  6  is  the  angle  between  T  and  and  Tq^  is  the  distance  from  the 

object  to  the  receiver.  The  Green’s  function  after  substituting  Eq.  (5-78)  can 
be  rewritten  as 


(5-79) 


G(r  .z)  =  CI 


•  ,  •,  ~®n*^or 

sin(YnZr)c  e 


''nVVT7 


sin(YjjZ)c 


-iK  rco<0)  a  rcos(0) 


/t  r  C  0 
*n/  “^or 


Substituting  Eq.  (5-79)  into  Eq.  (5-76)  yields  for  the  scattered  field 


(5-80) 


sln(Y  z  )c‘’^’'''’''e"“"'‘’^ 
Us(rs.z)=CI - S . r.  ^  - B. 


n 


or 
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where 


(5-81) 


Bn  = 


ffo)  — 

-iK„r  cos(e)  a„r  co<e) 
sin(7jjZ)c  e 

/,  rcos(e) 

•'o 

V  ^or  J 

sin(YjjZ)e 


-  iK^r  co<0)  a  r  cos(e) 


7^ 


r  cos(OI 


.^[f(p)].dS 


or 


Finally  combining  C  and  Bu  into  a  single  term  Ajj,  gives  for  the 
scattered  field 


(5-82) 


U,(r,,.)=XA„ 

n 


-o  „r 

^  • _ f  \  n  or  n  or 

Sin(YnZr)c  c 


Vv^n^or 


This  expression  satisfies  all  boundary  conditions  at  the  waveguide 
interfaces  and  is  continuous  throughout  all  space.  Note  that  the  solution, 

Eq.  (5-82),  is  in  the  form  of  a  guided  wave  as  one  would  expect  for  the  field  far 
from  the  source.  Now  that  the  scattered  far  field  has  been  determined,  some 
numerical  examples  will  be  presented. 

Examples 

In  this  example  we  will  follow  the  acoustic  field  as  it  propagates  from  the 
source  to  the  object,  to  the  receiver.  We  will  determine  the  total  field  at  the 
receiver  by  using  the  principle  of  superposition.  The  characteristics  of  the 
sample  waveguide  arc  shown  in  Table  3. 


TABLES 


WAVEGUIDE  CHARACTERISTICS 


Water  Layer 


Depth 

Compression 

Wave 

Density 

Attenuation 

(m) 

Sound  Speed 

(m/s) 

(g/cm**3) 

(nepers/m) 

150 

1500 

.  1.0 

1.0e*7 

Half  Space 

Depth 

Compression 

Wave 

Density 

Attenuation 

(m) 

Sound  Speed 

(m/s) 

(g/cm**3) 

(dBA) 

infinite 

1600 

1.5 

.5 

The  source  will  have  a  frequency  of  100  (Hz).  Table  4  lists  the  source,  object 
and  receiver  locations.  For  this  environment  there  are  7  propagating  modes, 
as  determined  by  Eq.  (5-36). 


TABLE 4 

LOCATION  PARAMETERS 


Source  Location 

Object  Location 

Receiver 

Location 

r  (m) 

z  (m) 

r  (m) 

z  (m) 

r  (m) 

z  (m) 

0 

50 

5000 

75 

6500 

25 

Figures  IS  and  16  depict  the  waveguide  in  a  vertical  and  horizontal  view.  The 
horizontal  view  is  looking  at  the  waveguide  from  the  surface  down. 


RECEIVER 

(6500,25) 


104 


r- 

z 
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SOURCE 


OBJECT 

(5000,75) 


Fig.  15.  Vertical  (side)  View  of  the  Waveguide. 


OBJECT 

(5000,75) 


RECEIVER 

(6500,25) 


Fig.  16.  Horizontal  (top)  View  of  the  Waveguide. 


The  source,  object  and  receiver  can  be  at  different  depths  in  the 
waveguide;  a  requirement  on  the  object  is  that  it  must  be  in  the  far  field  of  the 
source.  This  is  because  we  arc  assuming  that  only  the  discrete  propagating 
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modes  contribute  to  the  pressure  at  the  receiver  and  because  we  use  the 
asymptotic  expression  for  the  Hankel  function,  (in  the  solution  for  the  range 
equation)  which  is  only  valid  at  long  ranges.  Another  requirement  is  that  the 
object  is  not  so  near  an  interface  that  the  surface  on  which  the  near  field 
scattered  field  is  determined  comes  in  contact  with  the  interface.  The  objects 
characteristics  are  shown  in  Table  5. 

TABLES 

OBJECT  CHARACTERIST1C:S 
Length  (m)  Width  (m)  Boundary  Condition 
50  10  Dirichlet 

We  will  now  follow  the  field  as  it  propagates  in  the  waveguide.  Figure  17 
illustrates  the  magnitude  of  the  normalized  pressure  field  verses  depth  at  a 
range  of  300  (m).  Notice  that  the  pressure  at  the  surface  goes  to  zero  as  the 
boundary  condition  requires.  There  is  an  interference  pattern  due  to  the 
summation  of  the  7  individual  modes.  Figures  18*24  illustrate  the  contribution 
from  each  mode.  Note  that  mode  3  (Fig.  20)  contributes  the  least  to  the  overall 
field  and  that  Mode  6  (Fig.  23)  and  mode  7  (Fig.  24)  also  contribute  very  little. 
Modes  2  (Fig.  19)  and  5  (Fig.  22)  have  the  largest  contributions  to  the  field. 
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DEPTH  (m) 


0.000  0.002  0.004  0.006  0.008 


NORMALIZED  PRESURE 

Fig.  21.  Point  Source  Pressure  vs.  Depth  at  a  Range  of  300  (m) 
from  the  Source.  Mode  4  Only. 
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We  now  allow  the  field  to  propagate  to  the  object.  Figure  25  depicts  the 
magnitude  of  the  normalized  pressure  field  at  the  object,  5000  (m)  from  the 
source.  Notice  that  the  interference  pattern  differs  from  that  shown  in  Fig. 

17,  near  the  source.  This  is  because  the  higher  modes  attenuate  more  quickly 
than  do  the  lower  modes.  Notice  that  the  magnitude  of  the  field  has  decreased 
to  about  1/5  its  value  at  the  source.  Figures  26*32  depicts  the  contributions 
from  the  individual  modes.  Notice  that  only  the  magnitude  of  each  of  the 
modal  fields  has  changed,  not  the  shape.  Note  also,  that  the  relative  strength 
of  mode  6  (Fig.  31)  and  7  (Fig.  32)  have  reversed  because  of  ittenuation.  due  to 
interaction  with  the  waveguide  interfaces.  Modes  1  (Fig.  26)  and  2  (Fig,  27) 
now  contribute  more  to  the  field  than  any  other. 

We  now  look  at  the  scattered  pressure  field  as  it  propagates  from  the 
object.  Since  the  scattered  field  is  not  cylindrically  symmetric,  we  will  look  at 
the  pressure  in  the  vertical  plane  at  90  degrees  relative  to  the  front  of  the 
object  (backscatter)  and  0  degrees  (along  major  axis  of  symmetry).  Figure  33 
illustrates  the  vertical  normalized  pressure  field  at  a  distance  of  300  (m)  from 
the  object  and  at  an  angle  of  90  degrees  (backscatter).  Note  that  the  strength 
of  the  field  is  approximately  2  orders  of  magnitude  less  than  the  incident  field. 
Figures  34*40  illustr  ite  the  contributions  from  the  individual  modes.  Notice 
also  that  while  the  structure  of  the  modes  are  similar  to  the  incident  field 
modes,  their  relative  contributions  are  not  the  same  as  for  the  incident  field. 
Modes  1,  (Fig.  34)  3,  (Fig.  36)  5,  (Fig.  38)  and  7  (Fig.  40)  contribute  much  more 
than  the  corresponding  incident  field  modes.  This  is  because  the  object,  acting 
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NORMALIZED  PRESSURE 


Source  Pressure  vs.  Depth  at  a  Range  of  5000  (m) 
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Fig.  33.  Objects'  Scattered  Pressure  vs.  Depth  at  a  Range  of  300  (m) 
from  the  Object.  Scattered  Angle  90  (deg.) 
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+0  2.0C-6  4.0C-6  6.0e-6  8.0e-6  l.Oc-5  1.2e-5 


NORMALIZED  PRESSURE 

i'  Scattered  Pressure  vs.  Depth  at  a  Range  of  300  (m) 
Scattered  Angle  90  (deg.).  Mode  3  Only. 


DEPTH  (m) 


Fig.  38.  Objects'  Scattered  Pressure  vs.  Depth  at  a  Range  of  300  (m) 
from  the  Object.  Scattered  Angle  90  (deg.).  Mode  5  Only. 


+0  2.0e-6  4.0e-6  6.0e-6  8.0e-6  l.Oe-5  1.2e- 


NORMALIZED  PRESSURE 


Scattered  Pressure  vs.  Depth  at  a  Range  of  300  (m) 
Scattered  Angle  90  (deg.).  Mode  6  Only. 
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Fig.  40.  Objects'  Scattered  Pressure  vs.  Depth  at  a  Range  of  300  (m) 
from  the  Object.  Scattered  Angle  90  (deg.).  Mode  7  Only. 


as  a  source,  is  introduced  into  the  waveguide  and  the  energy  which  radiates 
from  it  must  behave  as  any  other  source  would  behave,  propagating  according 
to  the  restrictions  imposed  on  it  by  the  boundary  conditions.  These 
restrictions  manifest  themselves  through  the  eigenfunctions.  If  we  look  at 
the  7  depth  eigenfunctions  for  this  waveguide  (Figs.  41-47)  we  see  that  modes 
2,  (Fig.  42)  4,  (Fig.  44)  and  6  (Fig.  46),  each  have  a  null  at  or  very  near  the  mid¬ 
point  of  the  waveguide.  Modes  1,  (Fig.  41)  3,  (Fig.  43)  5,  (Fig.  45)  ai  (Fig.  47) 
have  a  maximum,  (or  a  minimum)  at  or  near  the  mid-point  of  the  waveguide. 
When  any  source  is  placed  at  the  mid-point  of  the  waveguide  the  odd 
numbered  modes  will  be  highly  excited  while  the  even  modes  will  not. 
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DEPTH  (m) 


-0.12  -0.08  -0.04  0.00  0.04  0.08  0.12 

MODE  7 

Fig.  47  Mode  7  Depth  Eigenfunction. 

This  is  what  is  happening  with  the  scattered  field. 

Figure  48  depicts  the  vertical  normalized  pressure  field  at  a  ran^e  of  300 
(m)  and  an  angle  of  0  degrees  from  the  object.  Notice  that  the  field  is 
approximately  1  order  of  magnitude  smaller  than  the  backscattered  field  and  3 
orders  of  magnitude  less  than  the  incident  field.  It  is  clear  that  the  scattered 
field  is  not  symmetric  about  the  object. 

It  is  customary  to  show  acoustic  propagation  in  terms  of  transmission 
loss  versus  range,  where  transmission  loss  is  -20  log  P  where  P  is  the 
normalized  pressure.  The  pressure  is  normalized  relative  to  the  pressure  from 
the  point  source  at  a  range  of  l(m).  Figure  49  shows  a  comparison  between  the 
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Fig.  48.  Objects'  Scattered  Pressure  vs.  Depth  at  a  Range  of  300  (m) 
from  the  Object.  Scattered  Angle  0  (deg). 

point  source  transmission  loss  and  the  forward  and  backscattered  transmission 
loss  of  the  object.  The  object  is  at  a  range  of  5000  (m)  from  the  source.  This  was 
obtained  by  solving  for  the  scattered  field  1500  (m)  from  the  object,  in  both 
the  forward  and  backward  direction.  This  graphically  shows  the  difference 
between  the  two  fields.  The  total  field  would  be  a  coherent  sum  of  the  two. 
Expressed  logarithmically  (in  dB),  there  would  be  little  or  no  change  in  the 
incident  field,  since  the  scattered  field  is  so  small. 
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Fig.  49.  Comparison  of  the  Transmission  Loss  vs.  Range  of  the  Point  Source  and 
the  Objects'  Scattered  Field. 


Figure  49  gave  a  two  dimensional  view  of  the  incident  and  scattered  field. 
We  now  turn  our  attention  to  representing  the  incident,  scattered  and  total 
field  in  three  dimensions.  The  idea  is  illustrated  in  Fig.  50.  A  portion  of  the 
waveguide  is  defined  by  the  ihree  dimensional  rectangular  region;  the  source 
is  located  somewhere  in  this  volume.  Next  we  superimpose  a  solution  grio  at 
some  depth  in  the  waveguide.  In  general  the  point  source  and  the  solution 
grid  will  not  be  at  the  same  depth.  At  each  point  of  the  grid,  we  determine  the 
pressure  and  consequently  the  transmission  loss  due  to  the  point  source.  An 
example  is  shown  in  Fig.  51.  The  source  is  at  a  depth  of  50  (m)  and  is  ai  the 
center  of  the  grid,  which  is  at  a  depth  of  25  (m).  Note  the  cylindrical 
symmetry  produced  by  the  point  source,  and  the  modal  interference  pattern. 
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Fig.  50.  Schematic  Drawing  Depicting  the  Geometry  for  a  Three  Dimensional 
Solution. 
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We  allow  the  field  shown  in  Fig.  51  to  propagate  in  the  waveguide  and 
allow  it  to  insonify  the  object.  We  compute  the  field  in  three  dimensions 
between  the  ranges  of  -1300  (m)  and  5700  (m)  with  a  cross  range  of  1175  (m); 
this  represents  a  grid  containing  2856  points.  Figure  52  depicts  the  results. 
Notice  how  smooth  the  field  is.  The  lack  of  structure  is  due  to  the  small 
number  of  modes  propagating.  There  is  still  a  small  amount  of  curvature  to 
the  wavefront.  Finally  we  will  place  the  object  at  a  range  of  5000  (m)  and 
solve  for  the  scattered  field.  The  results  are  illustrated  in  Fig.  53.  The  most 
obvious  feature  is  that  the  scattered  field  is  not  azimuthally  symmetric.  The 
major  axis  of  symmetry  of  the  object  is  parallel  to  the  cross  range  axis.  Notice 
that  the  forward  scattered  field  is  slightly  stronger  than  the  backscattered 
field  and  that  the  backscattered  field  is  broader  than  the  forward  scattered 
field.  The  forward  scattered  field  also  appears  to  have  a  high  pressure  lobe  on 
each  side  of  the  main  field.  And  finally  note  that  the  field  emanating  from 
either  end  of  the  object,  (90  degrees  relative  to  either  the  forward  or  backward 
scattered  field)  is  much  lower  than  any  other  direction. 
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Fig.  53.  Three  Dimensional  Objects’  Scattered  Field  Between  a  Range  of  4300  (m)  and  5700  (m) 
from  the  Source. 
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We  now  investigate  in  more  detail  the  scattered  field  of  the  object.  If  the 


object  is  truly  incorporated  into  the  waveguide  then  we  should  see  both  free 
field  and  waveguide  characteristics  in  this  field.  Figure  54  illustrates  the  free 
field  scattered  field  (near  field)  after  the  object  has  been  insonified  by  the  7 
modes.  The  field  is  expressed  as  transmission  loss.  The  smaller  the  number, 
the  stronger  the  signal.  The  arrow  indicates  from  which  direction  the 
incident  field  is  coming. 
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Fig.  54.  The  Objects'  Free  Field  Scattered  Field  vs  Angle,  After  7  Modes 
Insonified  It. 

Comparing  Fig.  54  with  Fig.  53  we  can  make  the  following  observations. 
First,  in  both  cases,  the  forward  scattered  field  is  narrower  than  the 
backscattcred  field.  Second,  from  Fig.  54  we  see  that  the  forward  scattered  field 
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does  indeed  have  well  pronounced  side  lobes.  Third,  in  both  cases  the  field  is 
the  weakest  along  the  major  axis  of  symmetry..  We  now  look  at  the  field 
produced  by  a  point  source  at  the  objects'  depth  of  75  (m).  Figure  55  illustrates 
the  field.  We  notice  that  this  field  peaks  at  two  separate  ranges.  The  first  is  at 
approximately  100  (m)  and  the  other  is  at  approximately  400  (m).  Referring  to 
the  field  of  the  object  in  Fig.  53.  we  notice  a  very  slight  undulation  of  both  the 
forward  and  back  field.  The  peak  of  this  undulation  is  occurring  at 
approximately  400  (m)  from  the  object.  This  would  correspond  to  the  second 
peak  in  Fig.  55.  The  first  peak  is  not  observed.  It  could  be  that  the  objects'  free 
field  effects  are  dominant  close  to  the  object  since  the  field  has  not  propagated 
long  in  the  waveguide.  This  indicates  that  we  have  coupled  the  free  field 
scattered  field  into  the  waveguide. 

We  now  coherently  add  the  incident  field.  Fig.  52,  and  the  scattered  field. 
Fig.  53  to  obtain  the  total  field,  which  is  illustrated  in  Fig.  56.  There  is  no 
discemable  difference  between  the  total  field  and  the  incident  field  of  Fig.  52. 
If,  however,  we  look  at  the  difference  (in  dB  space)  between  the  total  field  and 
the  incident  field,  we  will  be  left  with  just  the  change  in  the  field  due  to  the 
object.  Figure  57  shows  this  difference:  the  effect  of  the  object  is  limited  to  the 
forward  and  backward  direction.  The  largest  effect  is  observed  in  the 
backward  direction  due  to  the  fact  that  there  is  a  minimum  in  the  incident 
field  at  this  location.  The  difference  is  less  than  2  (dB).  In  the  forward 
direction  we  see  that  the  scattered  field  does  indeed  have  an  effect,  but  it  is 
very  slight,  only  1  (dB)  or  less. 

Having  obtained  the  total  field  for  a  coherent  addition  of  the  incident 
and  scattered  fields  for  a  single-layer  waveguide,  we  proceed  to  a  multilayered 
waveguide  in  which  the  sound  speed  will  vary  with  depth. 
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Fig.  55.  Three  Dimensional  Held  Produced  by  a  Point  Source  at  the  Objects’  Depth  of  75  (m). 
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Fig.  56.  Three  Dimensional  Total  Field. 
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Fig.  57.  Three  Dimensional  Difference  Field  Between  the  Total  Field  and  the  Incident  Field. 
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CHAPTER  6 

APPLICATION  OF  HUYGENS'  PRINCIPLE  TO  A  MULTILAYERED  WAVEGUIDE 

In  this  chapter  we  develop  the  scattered  field  for  an  object  in  a 
waveguide  consisting  of  multiple  isovelocity  water  and  sediment  layers. 
Because  the  general  development  is  the  same  as  in  Chapter  S,  we  will  only 
point  out  the  difference  between  the  two  expressions  due  to  the  differences  in 
the  waveguides.  We  conclude  this  chapter  with  an  example:  we  will  explore 
what  happens  to  the  field  when  the  object  is  insonified  by  only  a  few  of  the 
propagating  modes.  This  can  happen  when  the  source  is  in  a  duct  or  sound 
channel  and  the  object  is  outside  of  this  feature..  In  Appendix  E,  we  compare 
this  method  with  the  projection  method  described  method  in  Chapter  4. 

Theory 

The  development  of  the  field  due  to  a  harmonic  point  source  in  a 
waveguide  consisting  of  many  isovelocity  layers  is  given  by  Tolstoy.  (1.2)  We 
start  with  the  expression  for  the  acoustic  field,  previously  developed 


(6-1) 


P(r,z)=C 


C 


r 


where  Zn  is  the  depth  eigenfunction,  K„  is  the  horizontal  eigenvalue,  z  is  the 
depth  of  the  field  point  and  z^  is  the  source  depth.  The  horizontal  range 
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between  the  source  and  the  field  point  (receiver)  is  given  by  r  and  (Xn  is  the 
modal  attenuation  coefficient,  (See  appendix  D  for  further  information).  The 
attenuation  term  takes  into  account  energy  losses  due  to  interaction  with  the 
boundaries  and  losses  due  to  propagation  of  the  signal  in  the  water.  The 
geometry  of  the  waveguide  with  the  object  is  the  same  as  illustrated  in  Figs.  13 
and  14.  The  derivation  of  the  incident  field  is  the  same  as  for  the  isovelocity 
case  of  Chapter  5  except  that  the  depth  eigenfunctions  of  Equation  (6-1)  take  a 
different  form.  They  are  as  follows, 

(6-2)  Z^(zp=A„sin(Y„z^)  +  B„coj(Y„z^) 

(6-3)  Zj^(z)  =  AjjSin(YnZ)  +  Bjj  cos(YnZ) 

Where  Yn  >s  the  vertical  eigenvalue. 

The  expression  equivalent  to  Eq.(5-63)  is  the  following: 


(6-4)  P(r,z)  =  C  I 
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The  incident  Held  in  terms  of  the  spherical  coordinates  is  now  given  as 


U(B=I  X  <1", 

'  n  l,m.o 


where  n  is  the  index  specifying  the  number  of  modes  and  l.m.a  are  the  indices 
for  the  spherical  harmonics.  The  partial  waves  coefficients  (Equivalent  to  Eq. 
(5-74)),  of  the  incident  field  are  given  by  the  following. 


„  A„  s  i  n(Y„z  J  +  B„  COs(y„z  J  -a  r 

a  ,  =C - , ...  . - e 

V  n  SO 


.NAd- 

ll  2i  2  f  J 


A„  s  i  n(Y„z  J  +  B„  cos(y„zJ  -a  r 

^  n  *  n  s'  n  s  n  $o 

+  C“  '  """■  —  c 


2i  2  f  J  1.0 


where  Zg  is  the  depth  of  the  source,  Zobj  is  the  depth  of  the  object  and  tjo  is  the 
range  from  the  source  to  the  object.  The  scattered  field  about  the  object  has 


the  same  form  as  Eq.  (5-75)  namely. 
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(6-7) 


f(P)=  S 

n  9  I9  in  f  o 


where. 


^n,l,  a  l,r,  o  n,r,  a 


4'™  (P)  =h,(kp)  Y™  (0,9) 

If  O  1  M  ^ 


The  development  of  the  scattered  field  far  from  the  object  is 
equivalent  to  that  developed  in  chapter  5.  The  difference  is  due  to  the  fact  that 
we  now  have  different  depth  eigenfunctions  shown  in  Eqs.  (6-2)  and  (6-3). 

The  solution  for  the  scattered  field  far  from  the  object  is  the  same  as  Eq. 
(5-76)  namely. 


I  3G(r  ,  z)  3f(r) 

U3(r3..)  =  Jjf(r)-jJ - 


where  the  surface  a  is  chosen  at  a  suitable  region  circumscribing  the  object. 
The  scattered  field  is  written  as 


(6-9) 


U^(r^,z)  =  Cl-7===e  e 


( 
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where 


B  = 
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and  Zf  is  the  depth  of  the  receiver. 

This  expression  satisfies  all  boundary  conditions  at  the  waveguide 
interfaces  and  is  continuous  throughout  all  .space. 


Examplgii 

In  this  example  we  have  placed  the  source  and  receiver  in  a  weak  sound 
channel,  and  the  object  has  been  placed  near  the  bottom,  Fig.  58  illustrates  the 
geometry.  The  waveguide  has  35  water  layers  and  3  sediment  layers  overlying 
the  semi-infinite  half  space.  The  sound  speed  profile  is  shown  in  Fig.  59.  The 
source  is  again  at  100  (Hz).  There  are  16  propagating  modes  for  this 
environment. 
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Fig.  59.  Sound  Speed  Profile. 


The  characteristics  of  the  object  are  the  same  as  for  the  example  in 
Chapter  5.  We  now  allow  the  source  to  propagate  to  the  object.  Since  the 
source  is  located  at  the  apex  of  a  sound  channel  the  field  produced  by  the  lower 
order  modes  is  concentrated  in  this  channel.  This  means  that  with  the  object 
located  outside  the  sound  channel,  those  modes  that  remain  in  the  sound 
channel  will  have  littie  or  no  effect  on  the  object.  This  can  easily  be 
illustrated  by  again  isolating  the  modes  at  the  object  and  looking  at  the 
vertical  distribution  of  the  scattered  pressure  due  to  the  individual  modes.  In 
this  example  the  first  three  modes  have  very  little  effect  upon  the  object  (See 
Figs.  61-64)  with  mode  number  4  being  the  first  mode  to  in.sonify  the  object 
with  any  appreciable  energy.  We  therefore  look  at  the  vertical  pressure  at  the 


object  due  to  the  source  insonifing  it.  Figure  60  illustrates  the  field.  Notice 
that  the  field  peaks  at  the  source  depth  of  60  (m).  The  secondary  peak  results 
from  the  modal  contributions  which  are  not  contained  in  the  sound  channel. 
This  would  correspond  to  higher  order  modes.  The  field  diminishes  very 
quickly  after  this  secondary  peak. 


0.0000  0.0005  0.0010  0.0015 

NORMALIZED  PRESSURE 

Fig.  60.  Point  Source  Pressure  vs.  Dep’h  at  the  Object,  5000  (m) 
from  the  Source. 

We  will  now  look  at  the  effect  that  each  of  the  first  four  modes  have  at 
the  object.  Figures  61-64  shows  the  pressure  vs.  depth  at  the  object  for  the 


firs!  four  modes. 
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Notice  that  for  mode  1  (Fig.  61)  the  maximum  pressure  is  contained  within  the 
sound  channel.  The  pressure  dies  off  exponentially  after  a  depth  of  150  (m).  ^ 

This  is  to  be  expected  since  this  mode  is  evanescent  after  this  depth.  This  is 

also  true  for  mode  2,  (Fig.  62)  we  see  the  exponential  decay  of  the  pressure 
field  at  a  greater  depth  than  for  mode  one.  Again  this  is  consistent  with  mode  ^ 

theory,  which  states  that  the  higher  order  modes  will  travel  at  a  steeper  angle 
and  as  such  they  will  penetrate  through  the  sound  channel  more  than  the 
lower  order  modes..  Mode  three  (Fig.  63)  is  not  excited  by  the  point  source  and  ^ 

does  not  contribute  much  at  any  depth.  Finallv  we  see  that  for  mode  4  (Fig.  64) 
the  pressure  is  significantly  larger  at  the  objects  depth  of  260  (m).  The  field 
still  decays  exponentially  starting  approximately  at  the  depth  of  the  object.  I 

Now  we  allow  this  field  (all  16  modes)  to  insonify  the  object.  Figure  65 
shows  the  resulting  field  at  a  range  of  1500  (m)  in  the  forward  direction  (270 
degree  scattered  angle).  The  field  is  approximately  3  orders  of  magnitude  less  I 

than  the  incident  field.  One  can  see  that  the  field  is  very  complex,  resulting 

from  modal  interactions.  The  field  is  not  contained  in  the  sound  channel  since 
the  object,  which  is  generating  this  field,  is  not  within  the  sound  channel.  We  ^ 

will  now  look  separately  at  the  first  four  modes  contained  in  this  scattered 
field.  Figures  66-69  depict  the  pressure  vs.  depth  for  the  first  four  modes. 
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Fig.  65.  Objects'  Pressure  vs.  Depth  at  the  Receiver,  1500  (m)  from 
the  Object.  Scattered  Angle  270  (deg). 
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Fig.  67.  Objects’  Pressure  vs.  Depth  at  the  Receiver,  1500  (m)  from 
the  Object.  Scattered  Angle  270  (deg).  Mode  2  Only. 
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Fig.  69.  Objects'  Pressure  vs.  D-pth  at  the  Receiver,  1500  (in)  from 
the  Object.  Scattered  Angle  270  (deg).  Mode  4  Only. 


Notice  that  the  excitation  of  the  first  mode  (Fig.  66)  is  very  weak  through 
the  entire  depth  of  the  waveguide.  The  magnitude  of  the  normalized  pressure 
is  on  the  order  of  10'^.  Mode  2  (Fig.  67)  is  also  very  weak.  There  is  a  slight 
variation  with  depth  of  the  field.  One  would  no»  expect  that  these  two  modes 
would  be  excited  since  they  are  both  highly  evanescent  at  the  depth  of  the 
object.  We  see  that  mode  3  (Fig.  68)  is  excited  along  with  mode  4  (Fig.  69). 

These  and  the  higher  order  modes  are  not  affected  by  the  sound  channel.  This 

kind  of  modal  behavior  is  typical  of  a  point  .source.  Since  the  object  is 
producing  this  signal,  this  is  another  indication  that  the  object  is  coupled  to 
the  waveguide. 


Wc  now  compare  the  point  source's  and  object's  transmission  loss  with 


range.  Figure  70  illustrates  the  results. 
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Fig.  70.  Comparison  of  the  Transmission  Loss  vs.  Range  of  the  Point  Source  and 
the  Objects  Scattered  Field. 


The  point  source  field  (heavy  solid  line)  decays  very  little  with  range. 
This  is  consistent  with  our  observation  that  most  of  the  energy  is  contained 
within  the  .sound  channel.  A  sharp  contrast  to  this  is  the  scattered  field  of  the 
object.  It  varies  greatly  in  amplitude  over  a  smaller  horizontal  distance.  The 
scattered  field  is  approximately  70  (dB)  smaller  than  the  incident  field.  Again 
this  is  consistent  with  what  wc  saw  when  looking  at  the  first  four  modes. 
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We  now  look  at  the  three  dimensional  field.  Figure  71  depicts  the 
incident  field  at  a  range  from  4300  (m)  to  5700  (m),  with  the  cross-range 
coordinate  varying  from  -700  to  700  (m).  There  is  more  variability  of  the  field 
than  for  the  previous  case  (Fig.  52)  again  due  to  the  fact  that  more  modes  are 
propagating.  The  wave  front  is  slightly  curved,  as  is  to  be  expected  with  a  field 
that  is  cylindrically  symmetric.  Figure  72  depicts  the  scattered  field  for  the 
same  location  in  the  waveguide,  however  the  object  is  at  the  center  of  the  grid. 
We  notice  that  this  field  looks  very  similar  to  the  scattered  field  produced  in 
the  waveguide  of  the  la.st  chapter.  This  field  appears  to  be  increasing  in 
strength  as  it  propagates  from  the  object.  We  also  see  some  of  the  same 
features  as  in  the  last  chapter,  namely  that  the  field  is  strongest  in  the 
forward  and  backward  direction,  but  still  narrower  in  the  forward  direction. 
The  side  lobes  that  were  present  in  the  last  chapter  are  not  as  prevalent.  We 
compare  this  field  to  the  free  space  scattered  field  (Fig.  73)  that  was  produced 
by  allowing  the  object  to  be  insonified  by  the  propagating  modes  and  then 
compare  it  to  the  three  dimensional  field  (Fig.  74)  produced  by  the  point 
source  at  the  objects'  depth  in  the  waveguide.  We  do  this  for  the  same  reason 
that  we  did  it  in  Chapter  5,  namely  to  verify  that  we  have  correctly  coupled  the 
objects'  field  into  the  waveguide. 
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Fig.  71.  Three  Dimensional  Point  Source  Field  Between  a  Range  of  4300  (m)  and  5700  (m). 
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Fig.  72.  Threo  Dimensional  Objects’  Scattered  Field  Between  p  Range  of  4300  (m)  and  5700  (m) 
from  the  Source. 
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Fig.  73.  The  Objects'  Free  Field  Scattered  Field  vs.  Angle,  After  16  Modes  Have 
Insonified  It. 

Notice  that  the  free  space  scattered  field  is  similar  to  the  field  ob*.,ined  in 
Chapter  5.  There  are  differences,  however.  The  field  is  stronger  than  the  field 
in  Chapter  5  and  there  was  a  peak  in  the  backscattered  field  produced  by  the 
waveguide  of  Chapter  5  which  is  absent  for  this  waveguide.  ''  •.  side  lobes  are 
present  but  they  are  not  as  pronounced  as  in  the  last  chapter  which  is 
consistent  with  what  we  saw  in  Fig.  (72).  We  also  see  that  the  forward  scattered 
field  is  narrower  than  the  backscattered  field.  There  is  a  greater  difference 
between  the  forward  and  backward  fields.  This  difference  is  difficult  to  sec  in 
Fig.  72  but  it  is  evident  in  the  tv/o  dimensional  plot  (Fig.  70),  where  it  is  seen  to 
be  about  3  (dB).  Figure  (74)  gives  a  three  dimensional  illustiaiion  of  the  field 
produced  by  a  point  source  located  at  the  same  depth  as  the  object  The 
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waveguide  effect  on  the  field  allows  the  field  to  increase  over  the  range 
increment  as  was  observed  in  Fig.  71.  It  appears  that  the  field  of  the  object  is 
being  correctly  coupled  to  the  waveguide.  We  now  determine  the  total  field, 
which  is  shown  in  Fig.  75. 

There  is  no  observable  difference  between  the  total  field  and  the 
incident  field  of  Fig.  (70).  This  was  to  be  expected  because  the  large  difference 
in  amplitude  between  the  two  fields.  On  a  difference  plot  however  (Fig.  76), 
the  difference  between  the  incident  and  total  fields  is  evident.  Notice  that  the 
largest  difference  occurs  in  the  forward  direction  where  the  incident  field  has 
a  minimum.  This  difference  is  on  the  order  of  a  half  of  a  decibel. 


Fig.  75.  Three  Dimensional  Total  Field, 


Fig.  76.  Three  Dimensional  Difference  Field  Between  the  Total  Field  and  the  Incident  Field. 
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CONCLUSIONS 


We  have  developed  a  method  to  describe  acoustical  scattering  from  an 
object  in  a  waveguide.  We  have  done  so  by  using  Normal  Mode  theory  to 
describe  the  incident  field.  Each  mode  (at  the  object)  was  decomposed  into  a 
pair  of  up  going  and  down  going  plane  waves.  We  next  used  a  Transition 
matrix  to  determine  the  near  field  scattered  field  from  these  plane  waves 
incident  upon  the  object.  By  utilizing  Huygens'  principle  we  obtained  a  far 
field  solution  which  satisfied  all  boundary  conditions  and  preserves  continuity 
of  the  solution  throughout  all  space. 

We  have  shown  through  examples  that  the  object  is  correctly  coupled  to 
the  waveguide.  We  have  shown  that  the  objects  scattered  field  acts  as  a 
secondary  source  by  showing  that  this  scattered  field  obeyed  the  same 
boundary  conditions  as  the  point  source  field. 

This  method  of  determining  the  scattered  field  from  a  three  dimensional 
object  allows  one  to  determine  not  only  the  correct  target  strength  (intensity), 
but  also  to  properly  determine  the  phase.  In  addition  this  method  allows  one  to 
investigate  the  interaction  between  the  incident  field  with  the  waveguide,  the 
incident  field  with  the  object,  and  the  object's  scattered  field  with  the 
waveguide.  In  this  manner  a  better  understanding  of  the  ongoing  physical 
processes  can  be  obtained. 


APPENDIX  A:  COMPARISON  BETW'EEN  THE  E.  B.  C  METHOD  AND  AN  ANALYTICAL 
SOLUTION  FOR  THE  SCATTERED  FIELD  PRODUCED  BY  A  PLAN.,  WAVE  INCIDENT 

UPON  A  SPHERE. 


In  this  appendix  we  will  show  a  comparison  between  the  scattered  field 
as  solved  for  by  the  Extended  Boundary  Condition  method  and  the  analytic 
solution  for  the  case  when  a  plane  wave  is  incident  upon  a  sphere.  The 
analytic  expression  is  given  in  Morse  and  Feshbach  (1)  as: 
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where  k  is  the  wave  number,  r  is  the  distance  from  the  sphere,  8^^  is  the 

phase  angle,  0  is  the  observed  angle  and  a  is  the  radius  of  the  sphere.  The 
important  parameters  are  listed  in  Table  Al. 


TABLE  A1 

PROBLEM  PARAMETERS 


Sphere  Diameter  Boundary  Conditions  Distance  from  Sphere 

(m)  (m) 

2  Dirichlet  100 


The  "Distance  from  Sphere"  ..lue  in  Table  Al  is  the  distance  from  the  sphere 
that  the  scattered  field  was  obtained.  Three  frequencies  were  examined.  They 
were,  59.68  (Hz),  596.8  (Hz),  and  5968.  (Hz).  This  resulted  in  Kl/2  values  of  .25. 
2.5  and  25.  The  three  frequencies  were  selected  because  they  represented  a 


low,  medium  and  high  frequency  case.  Figures  Al,  A2,  and  A3  shows  the 
comparison  between  the  methods  for  the  different  frequencies  (the  forward 
direction  is  at  0  deg.).  In  all  cases  the  open  circles  are  the  analytic  solution 
while  the  black  dot  is  the  E.  B.  C.  solution.  Notice  in  all  cases  the  agreement  is 
quite  good.  The  E.  B.  C.  code  has  been  tested  thoroughly  by  Dr.  Mike  Werby  of 
the  Naval  Ocean  Research  and  Development  Activity.  This  particular  test  was 
performed  as  a  check  on  the  validity  of  the  E.B.C.  code  used  in  the  present 
work. 


O  ANALYTICAL  SOLUTION 
90  •  E.  B.  C.  SOLUTION 


Fig.  Al.  Comparison  Between  the  Analytic  Solution  and  the  E.  B.  C.  Solution  for 
a  Plane  Wave  Incident  Upon  a  Sphere  for  a  Frequency  of  59.68  (Hz). 


Fig.  A2.  Comparison  Between  the  Analytic  Solution  and  the  E.  B.  C.  S 
a  Plane  Wave  Incident  Upon  a  Sphere  for  a  Frequency  of  596.8  (Hz). 
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Fig.  A3.  Comparison  Between  the  Analytic  Solution  and  the  E.  B.  C.  Solution  for 
a  Plane  Wave  Incident  Upon  a  Sphere  for  a  Frequency  of  5968.  (Hz). 
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APPENDIX  B;  DERIVATION  OF  EQ.  (3-122) 


Wc  need  to  evaluate 


(B-1) 


aV  1  9  .  ,.2  2 


-  +i_^  +  k  a  H„  (ka„r)=? 
2  T  dr  n  0  ' 


where  (k  a^r)  =  Jp(k  a„r)  +  iYp(k  a^r) . 

Since  the  Bessel  function  is  analytic  for  all  ranges,  we  need  to  show  how  Eq. 
(B-1)  behaves  when  r  ^  0  and  when  r  =  0. 

First  consider  the  case  when  r  0.  The  derivative  of  the  Neuman  function 


with  respect  to  r  is 


(B-2) 


^  ■  fr(  “nO)  =  V„(lc  a„r)  . 


We  have  the  relationship, 

Equation  (B-2)  now  takes  the  following  form. 


(B-3) 
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Taking  the  second  derivative  of  the  Yq  gives, 

(B-4>  = 

d  r  y  Y  fka  r\^ 
dr  d(ka„r)l 

2 

We  have  used  the  following  relationship  concerning  the  derivative  of  the 
Neuman  function  of  the  first  order,  (1) 


(B-5) 


iY/x)4(Y,(x)-Y^(x)), 
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Now  letting  the  operator  of  Eq.  (B-1)  operate  on  iY^(kaj^r),  we  obtain,  with 
help  from  the  above  relationships, 


(B-6) 


ar^  ^ 


•Y«{ka„r) 


=  1 


=  1 


-k  a 


^2  2  1 
k 

n 

Y  _  "  Y  +1  n 

Y^ 

1 

<  2  j 

0  1  r  J  l\  2  J 

2 

The  recursion  formula  Y  (x)  =  ~-Y  (x)  -  Y^(x)  leads  to  the  following  for 

2  0 

(1) 


=  i 


k  a 

n 


k  a. 


0  \  r  J 


2^ 

iln 

V  2  y 


v'‘^ 


)- 


=  0 


Now  when  r=0,  we  will  use  the  small  r  approximation,  (2)  for  the  Neuman 
function. 


(B-8) 


-  |(  I  n(x)  +  Y  -  i  n(2)) ,  y  =  377- 


(B-9) 


for  small  r 


and  taking  the  derivative  with  respect  to  r  of  Eq.  (B-9)  gives 


(B-10) 


Y'  H-5LfY„(ka„r)'\=:  JJ2  I  J 
0  aA  OV  n  J)  T* 


2  M  ’tr 


now  doing  the  same  to  Eq,  (B-10)  gives. 


(B-11) 

Now  solving  Eq.  (B-1)  for  the  Neuman  function  at  small  r  gives. 


(B-12) 


5,2  '  ar  « 


+  kV 


\ay  { 


liv 

2 


) 


Multiplying  Eq.  (B-12)  by  r  and  integrating  gives 


(B-13) 


fO 

2 

Ldr^  '  ”J 

iY 


(ka„r)r<lr 


ka. 


kanf 


^  2 


(“'.'/ir 

— i 

1  2  ; 

16  J 

ka 


Summarizing  the  results. 
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(B-15) 

(B-16) 

(B-17) 

We  therefore 

(B-18) 

or 

(B-19) 


dr 


2  2 

2  r  dr  n 


1  d 


2  2 


0 


have  as  the  solution  for  Eq.  (B-1)  the  following. 


2  2 
+  +  a„ 

r  dr  n 


1  d 


4i»(r) 

2tcr 
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APPENDIX  C:  VALIDATION  OF  PROPAGATION  MODELS 


This  appendix  shows  the  comparison  between  the  horizontal  eigenvalues 
generated  by  the  S.ACLANT  Center  formal  Mode  A.coustic  propagation  (SNAP) 

(1)  model  and  both  the  isovelocity  normal  mode  model  which  was  developed  for 
the  Pekeris  waveguide  (i.e.  isovelocity  water  layer  over  a  high  speed 
isovelocity  fluid  half  space),  and  for  a  multilayered  waveguide,  consisting  of 
numerous  isovelocity  water  and  sediment  layers.  Table  (Cl)  lists  the 
waveguide  parameters. 


TABLE  Cl 

PARAMETERS  FOR  THE  ISOVELOCITY  WAVEGUIDE 


Depth  (m) 

Sound  Speed  (m/s) 

Denr-ity  (g/cm**3) 

0.0 

1500.0 

1.0 

150.0 

1500.0 

1.0 

150.0+ 

1600.0 

2.0 

Table  (C2)  shows  the  comparison  between  the  Pekeris  ws ‘eguide  model  and  the 
SNAP  model  determination  of  the  horizontal  eigenvalues,  with  the  magnitude 
of  the  difference  between  the  two  models. 


TABLE a 
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CXIMPARISON  OF  HORIZONTAL  EIGENVALUES  GENERATED  BY  THE  SNAP  AND 

PEKERIS  WAVEGUIDE  MODEL 


Pekeris  Waveguide 

SNAP  Program 

Magnitude 

Program 

of  Difference 

Horizontal 

Horizontal 

Eigenvalue 

Eigenvalue 

0.4184382396D+00 

0.4184382379D+00 

.lD-08 

0.4171019040D+00 

0.4171019027D+00 

.lD-08 

0.4148351550D+00 

0.4148351369D+00 

.lD-07 

0.4115958809D+00 

0.4115958645D+00 

.lD-07 

0.4073444101D+00 

0.4073443901D+00 

.3D-07 

0.402051 1324D+00 

0.40205 10729D+00 

.5D-07 

0.3957339603D+00 

0.3957339045D+00 

.5D-07 

For  this  example  there  were 

7  modes  propagating.  The  largest 

difference 

between  the  eigenvalues  were 

on  the  order  of  1x10'^  m*I  and  the  smallest 

difference  was  1x10'^  m’l. 

Figure  (Cl)  depicts  the  transmission 

loss  vs.  range 

for  the  environment  depicted  in  Table  (Cl).  Notice  that  the  structure  between 
the  two  models  are  identical  but  that  the  transmission  loss  generated  by  the 

Pekeris  Waveguide  model  is  off  by  a  factor  of  •  It's  obvious  that  the 


difference  is  just  a  constant  since  the  level  of  variation  does  not  change  with 
range.  It  is  obvious  that  the  Pekeris  Waveguide  model  will  give  as  accurate 
results  as  SNAP  for  these  rather  simplistic  environments. 


30 
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PEKERIS  WAVEGUIDE  MODEL 
SNAP  MODEL 


I  ■  '  '  I  '  '  »  ■  r  ■-»  I  ■  I  I . . . . *■ . *■  '  '  1 

0  1000  2000  3000  4000  5000 

RANGE  (m) 


Fig.  Cl.  Comparison  of  Transmission  Loss  Between  the  SNAP  and  Pekeris 
Waveguide  Model. 


We  compare  the  horizontal  eigenvalues  generated  by  the  Pekeris  Waveguide 
model  (see  chapter  5)  and  the  MultiLayered  Waveguide  model,  the  results  are 
as  in  Table  (C3). 


TABLE a 
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COMPARISON  OF  HORIZONTAL  EIGENVALUES  GENERATED  BY  THE  MULTILAYERED 

AND  PEKERIS  WAVEGUIDE  MODEL 


Multilayered  Waveguide 

Pekeris  Waveguide 

Magnitude 

Model 

Model 

of  Difference 

Horizontal 

Horizontal 

Eigenvalue 

Eigenvalue 

0.4184382396IM)0 

0.4184382396EM)0 

0. 

0.4171019040D+00 

0.417101904QD+00 

0. 

0.4148351S50&f00 

0.414835 1550D+00 

0. 

0.4115958809D+00 

0.4115958809D+00 

0. 

0.4073444101D+00 

0.4073444101D+00 

0. 

0.402051 1324D+00 

0.40205 11324D+00 

0. 

0.3957339603D+00 

0.3957339603D+00 

0. 

Out  to  10  decimal  places  there  is  no  difference.  Figures  (C2)  and  (C3)  shows  the 
comparison  between  the  depth  eigenfunctions  for  mode  1  and  mode  2,  as 
generated  by  the  Pekeris  Waveguide  and  the  MultiLayered  Waveguide  model. 


92 


Fig.  Cl.  Comparison  of  the  First  Depth  Eigenfunction  Generated  by  the  Pekeris 
and  MultiLayered  Model. 
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Fig.  C4.  Comparison  of  the  Transmission  Loss  Generated  by  the  Pekeris  and 
MultiLayered  Waveguide  Models. 


Figure  (C5)  shows  this  comparison.  Note  the  good  agreement  between  the  two 
models. 
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Fig.  C5.  Comparison  of  the  Transmission  Loss  Generated  by  the  MuItiLaycrcd 
and  SNAP  Model. 


With  confidence  in  the  MultiLayered  model,  we  now  compare  SNAP  and  the 
MultiLayered  model  for  a  slightly  more  complex  environment,  a  three  layered 
environment.  The  environment  is  depicted  in  Fig.  (C6),  with  the  water  layer 
having  the  lowest  sound  speed,  the  sediment  layer  the  next  and  the  half  space 
the  highest. 


DEPTH  (m) 


i 

Fig.  C6.  The  Second  Waveguide  Used  to  Validate  the  MultiLayered  Waveguide 
Model. 
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We  compare  the  iran&itiission  loss  generated  by  both  models  for  this 
environment  in  Fig.  (C7). 


Fig.  Cl.  Comparison  of  the  Transmission  Loss  Generated  by  the  MultiLayered 
Waveguide  Model  and  by  SNAP. 

We  now  investigate  a  much  more  complex  environment.  This  environment  is 
shown  in  Fig.  (C8).  It  consists  of  13  water  layers,  two  sediment  layers,  and  a 
half  space.  Figures  (C9)  and  (CIO)  compares  the  depth  eigenfunction 
generated  by  SNAP  and  by  the  MultiLayered  model  for  mode  1  and  mode  2  for 
this  environment.  They  show  good  agreement  between  the  two  models.  Figure 
(Cll)  shows  the  comparison  of  the  transmission  loss  between  SNAP  and  the 
MultiLayered  model  for  this  environment. 
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Fig.  C8.  The  Third  Test  Waveguide  Used  to  Validate  the  MuItiLayered  Waveguide 
Model. 
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Fig.  Cll.  Comparison  of  Transmission  Loss  Generated  by  the  MultiLayered  and 
SNAP  Model  for  the  Third  Test  Waveguide. 
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The  final  test  of  the  MulriLayered  model  is  against  data,  which  are  taken  from 
Ingenito  and  Wolf.  (2)  The  sound  speed  profile  is  shown  in  Fig.  (Cl 2). 


SOUND  SPEED  (m/s) 

Fig.  C12.  Sound  Speed  Profile  after  Ingenito  and  Wolf. 

The  model  vs.  data  comparison  for  4(X)  (Hz)  is  shown  in  Fig.  (C13),  along  with 
the  mode!  calculation  by  Ingenito  and  Wolf.  Note  that  the  agreement  is  quite 
good.  Table  (CS)  compares  the  horizontal  eigenvalues  generated  by  each  model 
and  gives  the  absolute  difference  between  the  two.  The  largest  difference  is 
lxlO'3  m-1. 
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NORMALIZED  mCENFUNCTION 


i)  b) 

Fig.  C13.  a)  Comparison  of  the  First  Depth  Eigenfunction  for  400  (Hz) 
Generated  by  the  MultiLayered  Waveguide  Model  and  Data,  b)  Comparison  of 
Model  Results  With  Data  as  Presented  by  Ingenito  and  Wolf. 


T/BLEC4 

COMPARISON  OF  HORIZONTAL  EIGENVALUES  GENERATED  BY  THE  MULTILAYERED 
WAVEGUIDE  MODEL  AND  THE  SNAP  MODEL  FOR  A  FREQUENCY  OF 400  (Hz) 


MultiLayered  Waveguide 

SNAP  Model 

Magnitude 

Model 

of  Different 

Horizontal 

Horizontal 

Eigenvalue 

Eigenvalue 

0.1643477015D+01 

0.1643473073D+01 

.30-05 

0.1632167076D+01 

0.16321S9020D^1 

.80-05 

C.1617969309D+01 

0.16179676630-^01 

.lD-05 

0.1S98583468D+01 

0.159857706SD-tOl 

.6D-05 

0.1572341994D401 

0.15723365910401 

.5D-05 

0.1549392552D^•01 

0.1549390991D401 

.10  05 

0.1534040267D!01 

0.1534034656D401 

.5D-05 

0.I499363280D^01 

0.1499358976D401 

.40-05 

0.1470642438D-i>01 

0.1470634678D401 

.7D-05 

0.143751 1941D+01 

0J437498305D401 

.lD-04 

0.1387679442D+01 

0.1387678929D401 

.5D-06 

0.1345397981D»>01 

0.13453887650401 

.9D-05 
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Figure  (C14)  show  the  comparison  between  the  transmission  loss  generated  by 
the  two  models.  The  agreement  is  quite  good. 


RANGE  (m) 


Fig.  CI4.  Comparison  of  Transmission  Loss  Generated  at  a  Frequency  of  400  (Hz) 
by  the  MultiLayered  and  SNAP  Model  Using  the  Sound  Speed  Profile  of  Fig.  (Cl  2). 


Figure  (C15)  shows  the  model  vs.  data  comparison  at  750  (Hz).  Again  the 
agreement  is  quite  good.  Table  (CS)  compares  the  horizontal  eigenvalues  as 
determined  by  each  model.  The  largest  difference  is  again  IxlO'^  m’^ 


NORMALIZED  DEPTH 


Fig.  C15.  a)  Comparison  of  the  First  Depth  Eigenfunction  for  750  (Hz) 
Generated  by  the  MulULayered  Waveguide  Model  and  Data,  b)  Comparison  of 
Model  Results  With  Data  as  Presented  by  Ingenito  and  Wolf. 


TABLE  C5 
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COMPARISON  OF  HORIZONTAL  EIGENVALUES  GENERATED  BY  THE  MULTILAYERED 
WAVEGUIDE  MODEL  AND  THE  SNAP  MODEL  FOR  A  FREQUENCY  OF  750  (Hz) 


MultiLayered  Waveguide 

SNAP  Model 

Magnitude 

Model 

of  Difference 

Horizontal 

Horizontal 

Eigenvalue 

Eigenvalue 

0.3086479467D+01 

0.3086475129D•^01 

.4D-05 

0.3077106218D+01 

0.3077095247D-»01 

.lD-04 

0.3064396577D+01 

0.3064383452D-H01 

.lD-04 

0.3054483253D+01 

0.3054480991D-t^l 

.2D-05 

0.3039869244D+01 

0.3039859840D-^01 

.9D-05 

0.30221 17375D+01 

0.30221 11962D-K)1 

.5D-05 

0.3000974575D+01 

0.3000965802D-»01 

.80^05 

0.2976025500D+01 

0.297601 6499D-H01 

.9D-05 

0.2948016738D+01 

0.2948006350D-^01 

.lD-04 

0.2930876943D+01 

0.2930875337D-I-01 

.lD-05 

0.2914863167D+01 

0.2914853243D-f01 

.9D-05 

0.28898858530-1-01 

0.2889880707D-H01 

.5D-05 

0.28726n077D-H01 

0.2872596332D-H01 

.lD-04 

0.2838535723D-»-P'l 

0.2838528599D-t-01 

.7D-05 

0.281 1841 177D-f01 

0.2811827591D-f01 

.lD-04 
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Figure  (Cl 6)  compares  the  transmission  loss  generated  by  the  two  models; 
agreement  is  quite  good. 


Fig.  Cl 6.  Comparison  of  Transmission  Loss  Generated  at  a  Frequency  of  750  (Hz) 
by  the  MultiLayered  and  SNAP  Model  Using  the  Sound  Speed  Profile  of  Fig.  (Cl  2). 


Finally  we  test  the  model  against  data  at  a  frequency  of  1500  (Hz).  The  comparison 
is  shown  in  Fig.  (C17).  Good  agreement  is  again  obtained.  Table  (C6)  shows  the 
comparison  between  the  horizontal  eigenvalues  generated  by  the  two  models.  The 
largest  difference  is  approximately  1x10*^  m*^. 


TABLE  C6 


COMPARISON  OF  HORIZONTAL  EIGENVALUES  GENERATED  BY  THE  MULTILAYERED 
V.'AVEGUIDE  MODEL  AND  THE  SNAP  MODEL  FOR  A  FREQUENCY  OF  1500  (Hz) 


McJtiLaycrcd  Wavcgaide 

SNAP  Model 

Magnitude 

Model 

of  Difference 

Horizontal 

Horizontal 

Eigenvalue 

Eigenvalue 

0.6178208386D+01 

0.6178205701D+01 

.2D-05 

0.6170829629D+01 

0.6170825078D+01 

.4D-05 

0.6161214393D+01 

0.61tl204743D+01 

.9D-05 

0.6147548359D+01 

0.6147527882D+01 

.2D-04 

0.6134282776D+01 

0.6134262291D+01 

.2D-04 

0.6126957355D+01 

0.6126951385D+01 

.5D-05 

0.6116767898D+01 

0.6116757846D+01 

.lD-04 

0.6105159734D+01 

0.6105 157746D+01 

.lD-05 

0.6091359606D+01 

0.6091353580D+01 

.6D-05 

0.6075254880D+01 

0.6075252494D+01 

.2D-05 

0.6057488415D+01 

0.60574801 64D+01 

.8D-05 

0.6037737287D+01 

0.603772633  lD+01 

.lD-04 

0.6016358465D+01 

0.6016350572D+01 

.7D-05 

0.5992942428D+01 

0.599293 1846D+01 

.lD-04 

0.5967737123D+01 

0.5967721 190D+01 

.lD-04 
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Figure  (C18)  shows  the  comparison  between  the  two  models.  And  we  see  again 
that  there  is  good  agreement  between  the  two  models. 


RANGE  (m) 


Fig.  Cl  8.  Comparison  of  Transmission  Loss  Generated  at  a  Frequency  of 
1500  (Hz)  L'y  the  MultiLayered  and  SNAP  Model  Using  the  Sound  Speed  Profile 
of  Fig.  (C12). 


Based  on  the  preceding  comparisons  of  horizontal  eigenvalues,  transmission 
loss,  and  eigenfunctions  we  conclude  that  both  the  Pekeris  Waveguide  model 
and  the  MultiLayered  Waveguide  model  can  generate  the  acoustic  pressure  at  a 
field  point  with  sufficient  accuracy  for  the  purpose  of  this  project. 


211 


APPENDIX  D:  MODAL  ATTENUATION  COEFHCIENTS  USED  IN  THE  NORMAL  MODE 

PROGRAMS. 


This  appendix  gives  a  brief  overview  of  the  modal  attenuation 
coefficients  used  in  the  normal  mode  propagation  models.  When  a 
propagating  mode  is  excited  it  will  decay  with  range  at  a  steady  rate.  Each 
mode  however  will  not  decay  at  the  same  rate.  The  higher  order  modes  will 
decay  more  quickly  than  the  lower  order  modes.  This  can  be  visualized  by 
remembering  that  the  higher  order  modes  propagate  at  a  higher  angle  and 
thus  interact  with  the  bottom  more  than  the  lower  order  modes.  When  mode 
decay  is  present,  the  wave  number  becomes  complex,  that  is  k  -4  k  +  i  p  where 
k  is  the  wave  number  (propagation  constant)  and  P  is  the  attenuation 
coefficient  of  a  plane  wave. 

Tindle  (1)  has  developed  an  expression  for  the  modal  attenuation 
coefficients  for  a  waveguide  consisting  of  an  isovclocity  water  layer 
overlaying  a  semi-infinite  isovelocity  half  space.  The  expression  is 


(Dl) 


“n  = 


(“P'  VzPla 


in  dB/m. 


Each  parameter  is  explained  below, 
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f=frequency 

(Hz) 

tO=27Cf 

(Hz) 

P=  .056f/1000. 

(dB/m) 

C2=  Compressional  Sound  Speed 
in  Half  Space 

(m/s) 

k2=Wavc  Number  in  Half  Space 

(m*l) 

Yin= Vertical  Component  of 

Wave  Number  in  Water  Layer 

(m-1) 

Y2n=Vertical  Component  of 

Wave  Number  in  Half  Space 

(m-1) 

H=Water  Depth 

(m) 

pj=Density  of  Water 

(kg/m3) 

P2=Density  of  Half  Space 

(kg/m3) 

This  expression  is  not  valid  when  the  waveguide  consists  of  multiple 
isovelocity  water  and  sediment  layers.  In  that  case  we  have  to  resort  to 
perturbation  theory  (2)  which  produces  the  following  expression, 


(D-2) 


n 


(z)dz 


0 


Ingenito  (3)  gives  a  brief  but  through  description  of  the  development  of  Eq. 
(D-2),  where  kn  is  the  horizontal  eigenvalue,  p  is  the  attenuation  coefficient  C 

is  the  compressional  sound  speed  and  Zn(z)  is  the  depth  eigenfunction.  The 
above  equation  is  valid  throughout  the  depth  of  the  waveguide. 
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The  attenuation  coefficient  for  the  water  layers  (Po)  is  given  by  the  following 


expression,  (4) 


(D-3) 


0.  OOTf  +0.155 


1.7f 


2  2 
(1.7)  +f  • 


where  the  attenuation  coefficient  for  the  sediment  layers  and  the  haifspace  is 
an  input  into  the  model.  Physically  what  happens  when  losses  are  assumed  to 
occur  primarily  in  the  sediment  is  that  the  bottom  boundary  condition 
requires  that  for  each  mode,  power  must  flow  from  the  water  into  the  bottom. 
Therefore  for  each  unit  horizontal  distance  traveled  the  the  total  power  loss 
for  a  mode  must  equal  the  power  transmitted  into  the  bottom  by  that  mode. 

Equation  (D-1)  was  used  in  the  Pekeris  Waveguide  Model  and  Eq.  (D-2) 
was  used  in  the  MultiLayerc '  Model  to  account  for  losses  in  the  sediment,  while 
Eq.  (D-3)  was  used  for  each  model  for  the  attenuation  coefficient  in  the  water. 
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APPENDIX  E:  COMPARISON  OF  THE  SCATTERED  FIELD  OBTAINED  USING  THE 
PRESENT  METHOD  AND  A  PROJECTION  METHOD. 


The  material  in  this  appendix  was  presented  at  the  2nd  Joint  Meeting  of 
the  Acoustical  Society  of  America  and  the  Acoustical  Society  of  Japan  which 
was  held  at  Honolulu,  Hi.  (1)  The  projection  method  was  developed  by  Evans 
(2)  and  has  been  presented  by  Norton.  (3-4)  The  characteristics  of  the 
scattering  object  are  presented  in  Table  El. 


TABLE  El 

OBJECT  CHARACTERISTICS 

Length  (m)  Width  (m)  Boundary  Conditions 
25  5  Dirichlet 


The  scattering  geometry  is  shown  in  Fig.  (El). 


10  (Km) 


SOURCE 


I 


FORWARD  SCATTERED 
SIGNAL 

- ^ 


Fig.  El.  Geometry  for  the  Scattering  Problem. 

Two  waveguides  were  investigated.  The  first  had  a  isovelocity  water  layer 
overlaying  an  isovelocity  halfspace.  The  second  waveguide  consisted  of  43 


water  layers  and  4  sediment  layers  over  a  halfspace.  The  first  waveguide  had 
the  source,  object  and  receiver  all  at  a  depth  of  75(m).  The  waveguide 
characteristics  are  listed  in  Table  E2. 


TABLE  E2 

ISOVELCXTTY  WAVEGUTOE  CHARACTERISTICS 
Water  Layer 


Depth 

Compressional 

Density 

Attenuation 

(m) 

Sound  Speed  (m/s) 

(g/cm**3) 

(nepers/m) 

150 

1500 

1.0 

1.0  E-6 

Half  Space 

Depth 

Compressional 

Density 

Attenuation 

(m) 

Sound  Speed  (m/s) 

(g/cm**3) 

(dBA) 

semi-inf 

1600 

2.0 

.15 

Figure  (E2)  shows  a  comparison  of  the  transmission  loss  vs  range 
between  the  two  methods  for  a  frequency  of  100  (Hz).  The  current  method 
based  on  Huygens'  principle  is  the  solid  line  and  the  projection  method  is  the 
dotted  line  They  are  not  identical  but  the  agreement  is  very  good.  Figure  (E3) 
is  a  comparison  between  the  two  for  the  same  waveguide  at  a  frequency  of  450 
(Hz).  Again  the  agreement  is  quite  good. 


TRANSMISSION  LOSS  (dB) 
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Fig.  £2.  Comparison  Between  the  Scattered  Held  Generated  From  the  Model 
Using  Huygens'  Principle  and  the  Projection  Method  for  a  Frequency  of  100  (Hz). 
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Fig.  E3.  Comparison  Between  the  Scattered  Held  Generated  From  the  Model 
Using  Huygens'  Principle  and  the  Projection  Method  for  a  Frequency  of  450  (Hz). 


Figure  (E4)  illustrates  the  multilayered  waveguide  used.  The  source 
object  and  receiver  are  all  at  60  (m).  Note  at  the  source  depth  we  have  a 
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minimum  in  the  sound  speed,  which  results  in  a  weak  sound  channel.  Figure 
(ES)  shows  a  comparison  between  the  scattered  field  expressed  as  transmission 
loss  vs  range  for  the  two  methods.  Also  plotted  in  this  case  is  the  field 
produced  by  a  point  source.  The  range  axis  is  therefore  relative  either  to  the 
point  source  or  to  the  location  of  the  object  depending  on  which  field  is  under 
consideration.  The  fleids  were  plotted  this  way  because  the  Feld  produced  by 
the  object  should  show  the  same  structure  as  that  produced  by  a  point  source, 
albeit  the  level  should  be  much  lower.  It  is  clearly  shown  that  not  only  are  the 
two  methods  agreeing  quite  well  but  they  both  show  structure  similar  to  a 
point  source.  In  Fig.  (E6)  we  show  the  results  for  a  frequency  of  450  (Hz). 

Again  we  have  excellent  agreement  between  the  two  methods. 


219 


Fig.  E4.  Sound  Speed  Profile  Used  m  the  Multilayered  Waveguide  Comparison 
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RANGE  (m) 


Fig.  ES.  Comparison  between  the  Field  I^uced  by  a  Point  Source  and  the 
Object  Using  the  Method  Based  on  Huygens'  Principle  and  the  Projection 
Method  fc.  a  Frequency  of  100  (Hz).  There  are  a  Total  of  47  Layers  in  the 
Waveguide. 
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Fig.  E6.  Comparison  Between  the  Scattered  Field  Generated  Using  the  Method 
Based  on  Huygens'  Principle  and  the  Projection  Method  for  a  Frequency  of 
450  (Hz).  There  are  a  Total  of  47  Layers  in  the  Waveguide. 


222 


ENDNOreS 

Notes  to  Pages  4-5 

CHAPTER  TWO:  FREE  SPACE  SCATTERING 

1.  James  E.  Burke,  "Low-Frequency  Scattering  by  Soft  Spheroids,"  J.  Acoust. 
ScSLAhL.  39  (1960),  826-831. 


2.  James  E.  Burke,  "Long-Wavelength  Scattering  by  Hard  Spheroids,"  L 
Acoust.  Soc.  Am..  40  (1966),  325-330. 


3.  James  E.  Buike,  "Scattering  by  Penetrable  Spheroids,"  J.  Acoust.  Soc. 
Am*.  43  (1968),  8/1-875. 


4.  C.  Yeh,  "Scattering  of  Acoustic  Waves  by  a  Penetrable  Prolate  Spheroid. 
1.  Liquid  Prolate  Spheroid,"  J.  Acoust.  Soc.  Am..  42  (1967),  518-521. 


5.  P.  Morse  and  H.  Fcshbach,  Methods  of  Theoretical  Physics.  (New  York: 
McGraw-Hill,  1953)  494-523. 


6.  H.  Levine  and  J.  Schwinger,  "  On  the  Theory  of  Diffraction  by  an 
Aperture  in  an  Infinite  Plane  Screen.  I,"  Phvs.  Rev.  74  (1948),  958-974. 


7.  K.  K.  Mei  and  J.  G.  Van  Bladel,  "Scattering  by  Perfectly-Conducting 
Rectangular  Cylinders,"  IEEE  Trans,  on  Antenna  and  Propagation.  Vol.  11 
(March  1963),  185-192. 


8.  Mogens  G.  Andreasen,  "Scattering  from  Parallel  Metallic  Cylinders  With 
Arbitrary  Cross  Sections."  IEEE  Trans,  on  Antenna  and  Propagation.  Vol.  12 
(November  1964),  746-754. 


9.  Mogens  G.  Andreasen,  "Scattering  from  Bodies  of  Revolution,"  IEEE 
Trans,  on  Antenna  and  Propagation.  Vol.  13  (March  1965),  303-310. 


10.  Robert  P.  Banaugh  and  Werner  Goldsmith,  "Diffraction  of  Steady 
Acoustic  Waves  by  Surfaces  of  Arbitrary  Shape,"  J.  Acoust.  Soc.  Am..  35  (1963), 
1590-1601. 


( 


II.  J.  Cesar  Monzon,  ”Three*Diinensionai  Scattering  by  an  Infinite 
Homogeneous  Anisotropic  Circular  Cylinder:  A  Spectral  Approach,"  IEEE  Trans. 
on  Antenna  and  Propagation.  Vol.  AP-35  (June  1987),  670-682, 


12.  L.  H.  Chen  and  D.  G.  Schweikert,  "  Sound  Radiation  from  an  Arbitrary 
Body,"  J.  Acoust.  Soc.  Am..  35  (1963),  1620-1632. 


13.  I.  Toyoda,  M.  Matsuhara,  and  N.  Kumagai,  "Extended  Integral  Equation 
Formulation  for  Scattering  Problems  from  a  Cylindrical  Scatterer,"  IEEE  Trans. 
on  Antenna  and  Propagation.  Vol.  36  (November  1988),  1580-1586. 


14.  P.  C.  Waterman,  "  New  Foundation  of  Acoustic  Scattering,"  J.  Acoust.  Soc. 
Am..  45  (1969),  1417-1429. 


15.  P.  C.  Waterman,  "  Matrix  Formulation  of  Electromagnetic  Scattering," 
Proc.  IEEE.  53  (1965),  805-812. 


16.  W.  R.  Smyth, "Charged  Right  Circular  Cylinder,"  J.  Appl.  Phvs.  27  (1956), 
917-920. 

17.  P.  C.  Waterman,  "  Matrix  Theory  of  Elastic  Wave  Scattering  II.  A  New 
Conservative  Law,"  J.  Acoust.  Soc.  Am..  63  (1977),  1320-1325. 


18.  B.  A.  Peterson,  V.  V.  Varadan,  and  V.  K.  Varadan,  "  T-Matrix  Approach  to 
Study  the  Vibration  Frequencies  of  Elastic  Bodies  in  Fluids,"  J.  Acoust.  Soc.  Am.. 
74  (1983),  1051-1056. 


19.  A.  Bostrum,  "Scattering  of  Stationary  Acoustic  Waves  by  an  Elastic 
Obstacle  Immersed  in  a  Fluid."  J.  Acoust.  Soc.  Am..  67  (1980),  390-398. 


20.  A.  Bostrum,  "Scattering  by  a  Smooth  Elastic  Obstacle,"  Tech.  Report  79-14, 
Inst.  Theorei.  Physics,  Goteborg  University,  Goteborg,  Sweden  (1979). 


21.  W.  M.  Visscher,  "A  New  Way  to  Calculate  Scattering  of  Acoustic  and 
Elastic  Waves,"  J.  Appl.  Phvs..  7  (1980). 


22.  Y.  H.  Pao,  "The  Transition  Matrix  for  the  Scattering  of  Acoustic  Waves 
and  Elastic  Waves,"  in  MMein  Problems  in  Elastic  Wave  Propagation.  (New 
York:  Wiley  and  Interscience,  1978). 


224 

Notes  to  Pages  6-11 

23.  B.  B.  Baker  and  E.  T.  Copson,  The  Mathematical  Theory  of  Huygens' 
Principle.  (Oxford:  Clarendon  Press,  1953). 


24.  M.  F.  Werby  and  S.  Chin-Bing,  "Some  Numerical  Techniques  and  Their 
use  in  the  Extension  of  T-Matrix  and  Null-Field  Approaches  to  Scattering,"  Int. 

J.  Comp.  Math.  Appls..  11  (1985),  717-731.  ( 


25.  M.  F.  Werby  and  G.  J.  Tango,  "Application  of  the  Extended  Boundary 
Condition  Equations  to  Scattering  from  Fluid-Loaded  Bounded  Objects,"  Engr. 

Analysis.  Vol.  5,  No.  1  (1988),  12-20. 

t 


26.  P.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics.  (New  York: 
McGraw-Hill.  1953),  884-886. 


27.  F.  Tricomi,  Integral  Equations  ,(New  York:  Doyer  Press,  1967). 

28.  M.  F.  Werby,  "  A  Coupled  Higher-Order  T-Matrix,"  J.  Acoust.  Soc.  Am., 
(submitted  1987). 


i 

29.  M.  F.  Werby,  G.  Tango  and  L.  H.  Green,  "Eigenyalue  and  Similarity 
Transformation  Methods  in  the  Solution  of  Acoustical  Scattering  Problems,"  1st 
IMACS  Symposium  on  Computational  Acoustics,  New  Hayen,  CT.,  to  appear  JL 
Como.  Math.  AddIs.  (1987). 

i 

30.  S.  Basket,  V.  V.  Varadan  and  V.  K.  Varadan,  "Thin  Shell  Theories  and 
Acoustic  Waye  Scattering  by  Infinitely  Long  Cylindrical  Shells  of  Arbitrary 
Cross  Section,"  J.  Acoust.  Soc.  Am..  75  (1984),  1673-1679. 


31.  M.  F.  Werby  and  L.  H.  Green,  "An  Extended  Unitary  Approach  for 
Acoustical  Scattering  from  Elastic  Shells  Immersed  in  a  Fluids"  J.  Acoust.  Soc. 
Am..  74  (1983),  625-630. 


32.  M.  F.  Werby  and  G.  Gaunaurd,  "Transition  from  Soft  to  Hard  Behayior  in 
Scattering  from  Submerged  Thin  Elastic  Shells,"  Acoustics  Letts.  9,  (1986),  111. 


33.  G.  Gaunard  and  M.  F.  Werby,  "Proper  Background  Choice  in  Resonance 
Scattering  from  Elastic  Shells,"  Int.  J.  Solids  Structures.  22,  (1986),  1149. 


Notes  to  Pages  11-20 


225 


34.  M.  F.  Werby  and  G.  Gaunard,  "Sound  Scattering  from  Spheroidal  Shells  in 
Water,"  in  Proceedings  ASME  Pressure  Vessels  Piping  Conference.  (New  York; 
ASME  Press.  1985). 


35.  M.  F.  Werby  and  L.  H.  Green,  "A  Comparison  of  Acoustical  Scattering 
from  Fluid-Loaded  Elastic  Shells  and  Sound  Soft  Objects."  J.  Acoust.  Soc.  Am..  76 
(1984),  1227-1230. 


36.  G.  Gaunaurd  and  M.  F.  Werby,  "Resonance  Response  of  Submerged 
Acoustically-Excited  Thick  and  Thin  Shells,"  J.  Acoust.  Soc.  Am..  77  (1985),  2081- 
2093. 


37.  Y.  H.  Pao  and  V.  Varatharajulo,  "Huygens'  Principle,  Radiation 
Condition,  and  Integral  Formula  for  the  Scattering  of  Elastic  Waves,"  J.  Acoust. 
SiiSLAJlL..  59  (1976),  1361-1370. 


38.  L.  R,  Dragonettc,  "Evaluation  of  the  Relative  Importance  of  Circum- 
Ferential  or  Creeping  Waves  in  Acoustic  Scattering  from  Rigid  and  Elastic  Solid 
Cylinders  and  Cylindrical  Shells,"  NRL  Report  8216  (1978). 


39.  M.  F.  Werby  and  G.  J.  Tango,  "Numerical  Study  of  Material  Properties  of 
Submerged  Elastic  Objects  Using  Resonance  Response,"  J.  Acoust.  Soc.  Am..  79 
(1986),  1260-1268. 


40.  M.  F.  Werby,  "Reply  to  'Comment  on  'A  Comparison  of  Acoustical 
Scattering  from  Fluid-Loaded  Elastic  Shells  and  Sound-Soft  Objects,"  J.  Acoust. 
SflSLAJBL.  81  (1987),  783-787. 


41.  M.  F.  Werby  and  L.  H.  Green,  "Correspondence  Between  Acoustical 
Scattering  from  Spherical  and  End-On  Incident  Spheroidal  Sl'^lls,"  J.  Acoust. 
ScSlAdl..  81  (1983),  783-787. 


42.  G.  C.  Gaunaurd  and  H.  Uberall,  "RST  Analysis  of  Monostatic  and  Bistatic 
Acoustic  Echoes  from  Spheres,"  J.  Acoust.  Soc.  Am..  73  (1983),  1-12. 


43.  L.  Flax,  G.  C.  Gaunaurd  and  H.  Uberall,  "Theory  of  Resonance  Scattering,' 
in  Physical  Acoustics  .  XV  (Edited  by  W.  P.  Mason  and  R.  N.  Thurston),  (New 
York:  Academic  Press,  1981)  Chapter  3,  191-194. 


226 

Notes  to  Pages  23-51 

CHAPTER  THREE:  ACOUSTIC  PROPAGATION  IN  A  WAVEGUIDE 

1.  Daljit  S.  Ahluwalia  and  Joseph  B.  Keller,  "  Exact  and  Asymptotic 
Representations  of  the  Sound  Field  in  a  Stratified  Ocean."  in  Wave  Propagation 
and.  Underwater  Acoustics.  Lecture  Notes  in  Physics  70,  edited  by  Joseph  B. 
Keller  and  John  S.  Papadakis,  (Berlin:  Springer-Verlag,  1977),  14-85. 


2.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition.  1985)  57-58. 


3.  H.  Goldstein,  Classical  Mechanics.  (Reading:  Addison-Wesley,  1980)  253- 
258. 


4.  Keith  R.  Symon,  Mechanics.  (Reading:  Addison-Wesley,  3rd  Edition, 
1971)  191-201. 


5.  Lord  Rayleigh,  Theory  of  Sound  11.  (New  York:  Dover  Publications,  1945) 
69-78 


6.  W.  Heitler,  The  Quantum  Theory  of  Radiation.  (Glasgow:  Oxford 
University  Press,  1954)  40-45. 


7.  M.  A.  Biot  and  1.  Tolstoy,  "Formulation  of  Wave  Propagation  in  Infinite 
Media  by  Normal  Coordinates  with  an  Application  to  Diffraction,"  J.  Acoust.  Soc. 
AUL.  29  (1957),  381-391. 


8.  1.  Tolstoy  and  C.  S.  Clay,  Ocean  Acoustics.  Theory  and  Experiment  in 
Underwater  Sound.  (New  York:  American  Institute  of  Physics.  1987)  98-105. 


9.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  '*rd.  Edition,  1985)  97. 


10.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  604-608. 


11.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  616-621. 


Notes  to  Pages  52-72 


227 


12.  I.  Tolstoy  and  C,  S.  Clay,  Ocean  Acoustics.  Theory  and  Experiment  in 
Underwater  Sound.  (New  York:  American  Institute  of  Physics,  1987)  36. 


13.  William  Elmore  and  Mark  Heald,  Physics  of  Waves.  (New  York:  Dover 
Publications,Inc.,  1969)  69. 


14.  C.  S. 


Clay  and  H.  Medwin,  Acoustical  Oceanography:  Principles  and 
,  (New  York:  John  Wiley  and  Sons,  1977)  295-296. 


15.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  104. 


16.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  454-463. 


17.  Daljit  S.  Ahluwalia  and  Joseph  B.  Keller,  "  Exact  and  Asymptotic 
Representations  of  the  Sound  Field  in  a  Stratified  Ocean,"  in  Wave  Propagation 
and  Underwater  Acoustics.  Lecture  Notes  in  Physics  70,  edited  by  Joseph  B. 
Keller  and  John  S.  Papadakis,  (Berlin:  Springer- Verlag,  1977),  21. 


18.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  57-58. 


19.  C.  S.  Clay  and  H.  Medwin,  Acoustical  Oceanography:  Principles  and 
Applications.  (New  York:  John  Wiley  and  Sons,  1977)  32-35. 


20.  P.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics.  (New  York: 
McGraw-Hill,  1953)  467. 


21.  W.  Magnus,  F.  Obclettinger  and  R.  P.  Soni,  Formulas  and  Theorems  for 
the  Special  Functions  of  Mathematical  Physics.  (New  York:  Springer- 
Vcrlog,1966)  82. 


CHAPTER  FOUR:  METHODS  DESCRIBING  SCATTERING  FROM  OBJECTS  IN  A 
WAVEGUIDE 

1.  Richard  B.  Evans,  "Acoustic  Scattering  from  a  Rigid  Object  in  a  Halfspace 
and  a  Slab  Waveguide,"  NORDA  Tech  Note  304.  (Jan.  1985). 


228 

Notes  to  Pages  72-76 

2.  M.  F.  Werby  and  R.  B.  Evans,  "Scattering  from  Objects  in  Unbounded  and 
Bounded  Oceans,"  IEEE  Jour.  Oceanic  Engr.  OE-12  (April  1987). 


3.  R.  H.  Hackman  and  G.  S.  Sammelmann,  "Acoustic  Scattering  in  an 
Inhomogeneous  Waveguide:  Theory,"  J.  Acoust.  Soc.  Am..  80,  (1986)  1447-1458. 


4.  G.  S.  Sammelmann  and  R.  Hackman,  "Acoustic  Scattering  in  a 
Homogeneous  Waveguide,"  J.  Acoust.  Soc.  Am..  82,  (1987)  324-336. 


S.  R.  Hackman  and  G.  S.  Sammelmann,  "Long-Range  Scattering  in  a  Deep 
Occanis  Waveguide,"  J.  Acoust.  Soc.  Am..  83,  (1988)  1776-1793. 


6.  Richard  B.  Evans,  Guy  V.  Norton,  and  M.  F.  Werby,  "A  Practical  Method 
for  Computing  Acoustic  Scattering  Caused  by  a  Submerged  Object  in  a 
Waveguide,"  NORDA  Tech  Note  349.  (Aug.  1987). 


7.  M.  D.  Collins  and  M.  F.  Werby,  "Scattering  from  an  ellipsoid  Submerged 
in  the  Ocean,"  J.  Acoust.  Soc.  Am.  Suppl.  1.82.  (Fall  87),  S105. 


8,  M.  D.  Collins  and  M.  F.  Werby,  "  A  Parabolic  Equation  Model  for 
Scattering  in  the  Ocean,"  J.  Acoust.  Soc.  Am.  85,  (1989)  1895-1902. 


9.  Rahul  Sen,  "  Scattering  of  Acoustic  Waves  in  a  Waveguide,"  J.  Acoust. 
Soc.  Am.  Suppl.  1.81.  (Spring  87),  S42. 


10.  A.  A.  Kleshchev  and  L.  L.  Klyukin,  "  Spectral  Characteristics  of  Sound 
Scattering  by  a  Body  in  a  Sound  Channel,"  Sov.  Phvs.  Acoust..  20,  (Nov. -Dec. 
1974)  143-145. 


11.  A.  A.  Kleshchev,  "Sound  Scattering  by  Spheroidal  Bodies  near  an 
Interface,"  Sov.  Phvs.  Acoust..  23,  (Nov. -Dec.  1977)  400-410. 


12.  A.  A.  Kleshchev,  "Sound  Scattering  by  a  Spheroidal  Body  at  an 
Interface,"  Sov.  Phvs.  Acoust..  25,  (Jan. -Feb.  1979)  143-145. 


13.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  57-58. 


I 


Notes  to  Pages  76-89 

14.  B.  B.  Baker  and  E.  T.  Copson,  The  Mathematical  Theory  of  Huygens' 
Principle.  (Oxford:  Clarendon  Press,  1953). 


CHAPTER  HVE:  APPLICATION  OF  HUYGENS’  METHOD  TO  AN  ISOVELOCITY 
WAVEGUIDE 


1.  C.  L.  Pekeris,  "Theory  of  Propagation  and  Explosive  Sound  in  Shallow 
Water."  Geol.  Soc.  Am.  Mem.  .  27,  (1948). 


2.  C.  S.  Clay  and  H.  Medwin,  Acoustical  Oceanography:  Principles  and 
Applications.  (New  York:  John  Wiley  and  Sons,  1977)  288-317. 


3.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  57-58. 


4.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  448-451. 


5.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  97. 


6.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition.  1985)  454-463. 


7.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  616-622. 


8.  G.  Arfken.  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  454-463. 


9.  C.  S.  Clay  and  H.  Medwin,  Acoustical  Oceanography:  Principles  and 
Applications.  (New  York:  John  Wiley  and  Sons,  1977)  61-64. 


10.  C.  S.  Clay  and  H.  Medwin,  Acoustical  Oceanography:  Principles  and 
Applications.  (New  York:  John  Wiley  and  Sons,  1977)  32-35. 


11.  C.  S.  Clay  and  H.  Medwin,  Acoustical  Oceanography:  Principles  and 
Applications.  (New  York:  John  Wiley  and  Sons,  1977)  295-296. 


Notes  to  Pages  94-202 


230 


12.  Guy  V.  Norton  and  M.  F.  Werby,  "  Some  Numerical  Approaches  to 
Describe  Acoustical  Scattering  from  Objects  in  a  Waveguide,"  Mathl.  Comout. 
Modeling.  11.  (1988)  81-86. 


13.  P.  M.  Morse  and  H.  Feshbach,  Metnods  of  Theoretical  Physics  11.  (New 
York:  McGraw-Hill.  1953)  1466. 


CHAPTER  SIX:  APPLICATION  OF  HUYGENS’  METHOD  TO  A  MULTILAYERED 
WAVEGUIDE 


1.  Ivan  Tolstoy,  "Note  on  the  Propagation  of  Normal  Modes  in 
Inhomogeneous  Media,"  J.  Acoust.  Soc.  Am. .27.  (March  1955)  274-277. 


2.  Ivan  Tolstoy.  "Resonant  Frequencies  and  High  Modes  in  Layered  Wave 
Guides,"  J.  Acoust.  Soc.  Am..28.  (November  1956)  1182-1192. 


APPENDDC  A:  COMPARISON  OF  THE  SCATTERED  FIELD  USING  THE  EXTENDED 
BOUNDARY  CONDITION  METHOD  AND  THE  ANALYTICAL  SOLUTION  FOR  A  PLANE 
WAVE  INCIDENT  ON  A  SPHERE. 


1.  P.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics  II.  (New  York: 
McGraw-Hill.  1953)  1483-1486. 


APPENDIX  B:  DERIVATION  OF  EQ.  (3-122). 


1.  G.  Arfken,  Mathematical  Methods  for  Physicists.  (Orlando:  Academic 
Press,  3rd.  Edition,  1985)  602. 


2.  M.  Abramowitz  and  I.  Stegun,  Handbook  of  Mathematical  Functions. 
(New  Yoric:  Dover  Publications,  Inc.,  1970)  360. 


APPENDD(  C:  VALIDATION  OF  PROPAGATION  MODELS 


1.  Finn  B.  Jensen  and  M.  C.  Ferla, 


Acoustic  Propagation  Model."  SACLANTCEN  Mem.  SM-121.  (Jan.  1979). 


2.  Frank  Ingenito  and  Stephen  N.  Wolf,  "  Acoustic  Propagation  an  Shallow 
Water  Overlaying  a  Consolidated  Bottom,"  J.  Acoust.  Soc.  Am. .60.  (September 
1976)  611-617. 


231 

Notes  to  Pages  211-214 

APPENDIX  D;  MODAL  ATTENUATION  COEFHCIENTS  USED  IN  THE  NORMAL  MODE 
PROGRAMS. 


1.  C.  Tindle.  "Attenuation  Parameters  from  Normal  Mode  Measurements."  L. 
Acoust.  Soc.  Am. .71.  (May  1982)  1145-1148. 


2.  K.  E.  Hawker.  W.  Williams  and  T.  Foreman.  "A  Study  of  the  Acoustical 
Effects  of  Sub-Bottom  Absorption."  J.  Acoust.  Soc.  Am. .65.  (  1979)  360-367. 


3.  Finn  B.  Jensen  and  M.  C.  Ferla.  "SNAP:  The  SACLANTCEN  Normal  Mode 
Acoustic  Propagation  Model."  SAQ-ANTCEN  Mem.  SM-121.  (Jan.  1979). 


4.  Frank  Ingenito.  "Measurements  of  Mode  Attenuation  Coefficients  in 
Shallow  Water."  J.  Acoust.  Soc.  Am..53.  (  1973)  858-863. 


APPENDIX  E:  COMPARISON  OF  SCATTERED  FIELD  DEVELOPED  USING  PRESENT 
METHOD  AND  A  PROJECTION  METHOD. 


1.  Guy  Norton.  M.  F.  Werby,  and  M.  D.  Collins.  "Benchmaric  Comparison  of 
Underwater  Acoustic  Scattering  Models."J.  Acoust.  Soc.  Am..  Suppl.  1  84,  (  Fall 
1988)  S219. 


2.  Richard  B.  Evans,  Guy.  V.  Norton,  and  M.  F.  Werby,  "A  Practical  Method 
for  Computing  Acoustic  Scattering  Caused  by  a  Submerged  Object  in  a 
Waveguide,"  NORDA  Tech  Note  349.  (Aug.  1987). 


3.  Guy  V.  Norton  and  M.  F.  Werby,  "  Some  Numerical  Approaches  to 
Describe  Acoustical  Scattering  from  Objects  in  a  Waveguide,"  Mathl.  Comput. 
Modeling.  11,  (1988)  81-86. 


4.  Guy  V,  Norton,  "A  Numerical  Approach  to  Describe  Acoustical  Scattering 
from  an  Object  in  a  Waveguide  Using  a  Projection  Technique,"  Second  IMACS 
Symposium  on  Computational  Acoustics.  (March  14-17,  1989). 


232 


LIST  OF  REFERENCES 


Abramowitz,  M.,  and  I.  Stegun.  1970.  Handbook  of  Mathematical  Functions.  New 
Yo:k:  Dover  Publications,  Inc. 

Ahluwalla,  Ouljit  S.,  and  Joseph  B.  Keller.  1977.  Exact  and  Asymptotic 

Representations  of  the  Sound  Field  in  a  Stratified  Ocean.  In  Wave 
Propagation  and  Underwater  Acoustics,  lActurfc  Notes  in  Physics  70.  ed. 
Joseph  B.  Keller  and  John  S.  Papadakis,  14-85.  Berlin:  Springer-Verlag. 

Andreasen,  Mogens  G.  1964.  Scattering  from  Parallel  Metallic  Cylinders  With 
Arbitrary  Cross  Sections.  IEEE  Trans,  on  Antenna  and  Propagation  Vol. 
12:746-754. 

_ .  1965.  Scattering  from  Bodies  of  Revolution.  IEEE  Trans,  on 

Arfkcn,  G.  1985,  Tvlathematical  Methods  for  Physicists.  3rd  ed.  Orlando: 

Academic  Press. 


Baker,  B.  B.,  and  E.  T.  Copson.  1953. 

Principle.  Oxford:  Clarendon  Press. 


Banaugh,  Robert  P.,  and  Werner  Goldsmith.  1963.  Diffraction  of  Steady 

Acoustic  Waves  by  Surfaces  of  Arbitrary  Shape.  J.  Acoust.  Soc.  Am. 
35:1590-1601. 


Basker,  S.,  V.  V.  Varadan,  and  V.  K.  Varadan.  1984.  Thin  Shell  Theories  and 
Acoustic  Wave  Scattering  by  Infinitely  Long  Cylindrical  Shells  of 
Arbitrary  Cross  Section.  J.  Acoust.  Soc.  Am.  75:1673-1679. 

Biot,  M.  A.,  and  I.  Tolstoy.  1957.  Formulation  of  Wave  Propagation  in  Infinite 
Media  by  Normal  Coordinates  with  an  Application  to  Diffraction.  L. 
Acoust.  Soc.  Am.  29:381-391. 

Bostrum,  A.  1979.  Scattering  by  a  Smooth  Elastic  Obstacle.  Tech.  Report  79-14. 
Inst.  Theorct.  Physics,  Goteborg  University,  Goteborg,  Sweden. 

_ .  1980.  Scattering  of  Stationary  Acoustic  Waves  by  an  Elastic 

Obstacle  Immersed  in  a  Fluid.  J.  Acoust.  Soc.  Am.  67:390-398. 


233 

Burke.  James  E.  1960.  Low-Frequency  Scattering  by  Soft  Spheroids.  J.  Acoust. 
SjulAhl  39:826-831. 

_ .  1966.  Long-Wavelength  Scattering  by  Hard  Spheroids.  J.  Acoust. 

Soc.  Am.  40:325-330. 

_ .  1968.  Scattering  by  Penetrable  Spheroids.  J.  Acoust.  Soc.  Am. 

43:871-875. 

Chen,  L.  H.,  and  D.  G.  Schwcikert.  1963.  Sound  Radiation  from  an  Arbitrary 
Body.  J.  Acoust.  Soc.  Am.  35:1620-1632. 

Clay.  C.  S..  and  H.  Mcdwin.  1977.  Acoustical  Oceanography.  Principles  and 
Applications.  New  York:  John  Wiley  and  Sons. 

Collins.  M.  D..  and  M.  F.  Werby.  1987.  Scattering  from  an  ellipsoid  Submerged 
in  the  Ocean.  J.  Acoust.  Soc.  Am.  Suppl.  1  82:S105. 

_ .  and  M.  F.  Werby.  1989.  A  Parabolic  Equation  Model  for  Scattering 

in  the  Ocean.  J.  Acoust.  Soc.  Am.  85:1895-1902. 

Dragonette.  L.  R.  1978.  Evaluation  of  the  Relative  Importance  of  Circum- 
Ferential  or  Creeping  Waves  in  Acoustic  Scattering  from  Rigid  and 
Elastic  Solid  Cylinders  and  Cylindrical  Shells.  NRL  Rep.  8216. 

Elmore.  William,  and  Mark  Heald.  1969.  Physics  of  Waves.  New  York:  Dover 
Publications. Inc.. 

Evans.  Richard  B.  1985.  Acoustic  Scattering  from  a  Rigid  Object  in  a  Halfspace 
and  a  Slab  Waveguide.  NORDA  Tech  Note  304. 

_ .  Guy  V.  Norton,  and  M.  F.  Werby.  1987.  A  Practical  Method  for 

Computing  Acoustic  Scattering  Caused  by  a  Submerged  Object  in  a 
Waveguide.  NORDA  Tech  Note  349. 

Flax.  L..  G.  C.  Gaunaurd  and  H.  Uberall.  1981.  Theory  of  Resonance  Scattering. 

In  Physical  Acoustics,  ed.  W.  P.  Mason  and  R.  N.  Thurston,  191-194,  vol. 
XV.  New  York:  Academic  Press. 

Gaunard,  G.,  and  H.  Uberall.  1983.  RST  Analysis  of  Monostatic  and  Bistatic 
Acoustic  Echoes  from  Spheres.  J.  Acoust.  Soc.  Am.  73:1-12. 

_ ,  and  M.  F.  Werby.  1985.  Resonance  Response  of  Submerged 

Acoustically-Excited  Thick  and  Thin  Shells.  J.  Acoust.  Soc.  Am. 
77:2081-2093. 

_ ,  and  M.  F,  Werby.  1986.  Proper  Background  Choice  in  Resonance 

Scattering  from  Elastic  Shells.  Int.  J.  Solids  Structures  22- 1149. 


Gok:  jin,  H.  1980.  Classical  Mechanics.  Reading;  .Addison- Wesley. 


234 

Kacktnan,  R.  H.,  and  G.  S.  Sammeimann.  1986.  Acoustic  Scattering  in  an 

Inhomogeneous  Waveguide;  Theory.  J.  Acoust.  Soc.  Am.  80:1447-1458. 

_ .  and  G.  S.  Sammeimann.  1988.  Long-Range  Scattering  in  a  Deep 

Oceanis  Waveguide.  J.  Acoust.  Soc.  Am.  83:1776-1793. 

Hawker.  K.  E..  W  Williams,  and  T.  Foreman.  1979.  A  Study  of  the  Acoustical 
Effects  of  Sub-Bottom  Absorption.  J.  Acoust.  Soc.  Am.  65:360-367. 

Hcitler.  W.  1954.  The  Quantum  Theory  of  Radiation.  Glasgow:  Oxford  University 
Press. 


Ingenito.  Frank.  1973.  Measurements  of  Mode  Attenuation  Coefficients  in 
Shallow  Water.  J.  Acoust.  Soc.  Am.  53:858-863. 

_ .  and  Stephen  N.  Wolf.  1976.  Acoustic  Propagation  an  Shallow 

Water  Overlaying  a  Consolidated  Bottom.  J.  Acoust.  Soc.  Am.  60:611-617. 

Jensen.  Finn  B..  and  M.  C.  Ferla.  1979.  SNAP:  The  SACLANTCEN  Normal  Mode 
Acoustic  Propagation  Model  SACLANTCEN  Mem.  SM-121. 

Kleshchev.  A.  A.,  and  L.  L.  Klyukin.  1974.  Spectral  Characteristics  of  Sound 

Scattering  by  a  Body  in  a  Sound  Channel.  Sov.  Phvs.  Acoust.  20:143-145. 

_ .  1977.  Sound  Scattering  by  Spheroidal  Bodies  near  an  Interface. 

Sov.  Phvs.  Acoust.  23:400-410. 

_ .  1979.  Sound  Scattering  by  a  Spheroidal  Body  at  an  Interface.  Sov. 

Phvs.  Acoust.  25:143-145. 

Levine.  H..  and  J.  Schwinger,  1948.  On  the  Theory  of  Diffraction  by  an 
Aperture  in  an  Infinite  Plane  Screen.  1.  Phvs.  Rev.  74:958-974. 


Magnus.  W..  F.  Obelettinger  and  R.  P.  Soni.  1966.  Formulas  and  Theorems  for 

Physics.  New  York;  Springer- 


Verlog. 

Mei.  K.  K..  and  J.  G.  Van  Bladel.  1963.  Scattering  by  Perfectly-Conducting 
Rectangular  Cylinders.  IEEE  Tians-  QH  AnKnna  an!  Propagation 
11:185-192. 


Monzon.  J.  Cesar.  1987.  Three-Dimensional  Scattering  by  an  Infinite 

Homogeneous  Anisotropic  Circular  Cylinder:  A  Spectral  Approach.  IEEE 
Tcans.  lan  Antenna  and.  Propagation  Voi.  AP-35:670-682. 

Morse.  P..  and  H.  Feshbach.  1953,,  Methods  of  Theoretical  Physics.  New  York; 
McGraw-Hill. 


( 


i 


i 


i 


4 


i 


t 


235 

Norton.  Guy  V.,  and  M.  F.  Wcrby.  1988.  Some  Numerical  Approaches  to  Describe 
Acousii  ai  Scattering  from  Objects  in  a  Waveguide.  Mathl.  Compui. 
Modeling  11:81-86. 

_ ,  M.  F.  Wcrby,  and  M.  D.  Collins.  1988.  Benchmark  Comparison  of 

Underwater  Acoustic  Scattering  Models.  J.  Acoust.  Soc.  Am..  Suppl.  1 
84:S219. 

_ .  1989.  A  Numerical  Approach  to  Describe  Acoustical  Scattering 

from  an  Object  in  a  Waveguide  Using  a  Projection  Technique.  Second 

IMACS  Symposium  on  Compuiaiional  Acoustics. 

Pao,  Y.  H.,  and  V.  Varatharajulo.  1976.  Huygens'  Principle.  Radiation  Condition, 
and  Integral  Formula  for  the  Scattering  of  Elastic  Waves.  J.  Acoust.  Soc. 
Am.  59:1361-1370. 

_ .  1978.  The  Transition  Matrix  for  the  Scattering  of  Acoustic  Waves 

and  Elastic  Waves.  In  Modern  Problems  in  Elastic  Wave  Propagation. 

New  York:  Wiley  and  Interscicnce. 

Pekeris,  C.  L.  1948.  Theory  of  Propagation  and  Explosive  Sound  in  Shallow 
Water.  Gcol.  Soc.  Am.  Mem.  27. 

Peterson,  B.  A.,  V.  V.  Varadan,  and  V.  K.  Varadan.  1983.  T-Matrix  Approach  to 
Study  the  Vibration  Frequencies  of  Elastic  Bodies  in  Fluids.  J.  Acoust. 
SiliLAlIL  74:1051-1056. 

Rayleigh,  Lord.  1945.  Theory  of  Sound  II  New  York:  Dover  Publications. 

Sammelmann,  G.  S..  and  R.  Hackman.  1987.  Acoustic  Scattering  in  a 
Homogeneous  Waveguide.  J.  Acoust.  Soc.  Am.  82:324-336. 

Sen,  Rahul.  ,-?87.  Scattering  of  Acoustic  Waves  in  a  Waveguide.  J.  Acou.st.  Soc. 
Am.  Suppl.  1  81:S42. 

Smyth,  W.  R.  1956.  Charged  Right  Circular  Cylinder.  J.  Appl.  Phvs.  27:917-920. 

Symon,  Keith  R.  1971.  Mechanics.  3d  ed.  Reading:  Addison-Wesley. 

Tindle,  C.  1982.  Attenuation  Parameters  from  Normal  Mode  Measurements.  jL 
Acoust.  Soc.  Am.  71:1145-1148. 

Tolstoy,  Ivan.  1955.  Note  on  the  Propagation  of  Normal  Modes  in 
Inhomogeneous  Media.  J  Acoust.  Soc.  Am.  27:274-277. 

_ .  1956.  Resonant  Frequencies  and  High  Modes  in  Layered  Wave 

Guides.  J.  Acoust.  Soc.  Am.  28:1182-1192. 

_ .  and  C.  S.  Clay.  1987.  Ocean  Acoustics.  Thcarv  and  Experiment  in 

Underwater  Sound.  New  York:  American  Institute  of  Physics. 


Toyoda,  I..  M.  Matsuhara.  and  N.  Kumagai.  1988.  Extended  Integral  Equation 

Formulation  for  Scattering  Problems  from  a  Cylindrical  Scatterer.  IEEE 
Trans,  on  Antenna  and  Propagation  Vol.  36:1580-1586. 

Tricomi,  F.  1967.  Integral  Equations.  New  York:  Dover  Press. 

Visscher.  W.  M.  1980.  A  New  Way  to  Calculate  Scattering  of  Acoustic  and  Elastic 
Waves.  L.  Appl.  Phvs.  7. 

Waterman.  P.  C.  1965.  Matrix  Formulation  of  Electromagnetic  Scattering.  Proc. 
IEEE  53:805-812. 

_ .  1969.  New  Foundation  of  Acoustic  Scattering.  J.  Acoust.  Soc. 

Am.  45:1417-1429. 

_ .  1977.  Matrix  Theory  of  Elastic  Wave  Scattering  II.  A  New 

Conservative  Law.  J.  Acoust.  Soc.  Am.  63:1320-1325. 

Werby,  M.  F..  and  L.  H.  Green.  1983.  An  Extended  Unitary  Approach  for 

Acoustical  Scattering  from  Elastic  Shells  Immersed  in  a  Fluids.  J.  Acoust. 

74:625-630. 

_ ,  and  L  H.  Green.  1983.  Correspondence  Between  Acoustical 

Scattering  from  Spherical  and  End-On  Incident  Spheroidal  Shells.  L. 
Acoust.  Soc.  Am.  81:783-787. 

_ .  and  L.  H.  Green.  1984.  A  Comparison  of  Acoustical  Scattering  from 

Fluid-Loaded  Elastic  Shells  and  Sound  Soft  Objects.  J.  Acoust.  Soc.  Am. 
76:1227-1230. 

_ .  and  S.  Chin-Bing.  1985.  Some  Numerical  Techniques  and  Their 

use  in  the  Extension  of  T-Matrix  and  Null-Field  Approaches  to 
Scattering.  Int.  J.  Comp.  Math.  Appls.  11:717-731. 

_ ,  and  G.  Gaunard.  1985.  Sound  Scattering  from  Spheroidal  Shells  in 

Water.  In  Proceedings  ASME  Pressure  Vessels  Piping  Conference.  New 
York:  ASME  Press. 

_ .  and  G.  Gaunaurd.  1986.  Transition  from  Soft  to  Hard  Behavior  in 

Scattering  from  Submerged  Thin  Elastic  Shells.  Acoustics  Letts.  9:111. 

_ .  and  G.  J.  Tango.  1986.  Numerical  Study  of  Material  Properties  of 

Submerged  Elastic  Objects  Using  Resonance  Response.  J.  Acoust.  Soc. 

Am.  79:1260-1268. 

_ .  1987.  A  Coupled  Higher-Order  T-Matrix.  J.  Acoust.  Soc.  Am. 

submitted. 

_ ,  G.  Tango,  and  L.  H.  Green.  1987.  Eigenvalue  and  Similarity 

Transformation  Methods  in  the  Solution  of  Acoustical  Scattering 
Problems.  1st  IMACS  Symposium  on  Computational  Acoustics.  New 
Haven,  CT.,  to  appear  L Comp.  Math.  Appls.  , 


237 


_ .  1987.  Reply  to  'Comment  on  'A  Comparison  of  Acoustical 

Scattering  from  Fluid-Loaded  Elastic  Shells  and  Sound-Soft  Objects. 

Acoust.  Soc.  Am.  81:783-787. 

_ ,  and  R.  B.  Evans.  1987.  Scattering  from  Objects  in  Unbounded  and 

Bounded  Oceans.  IEEE  Jour.  Oceanic  Engr.  OE-12. 

_ ,  and  G.  J.  Tango.  1988.  Application  of  the  Extended  Boundary 

Condition  Equations  to  Scattering  from  Fluid-Loaded  Bounded  Objects. 
Engr.  Analysis.  Vol.  5.  No.  1:12-20. 

Yeh,  C.  1967.  Scattering  of  Acoustic  Waves  by  a  Penetrable  Prolate  Spheroid.  I. 
Liquid  Prolate  Spheroid.  J.  Acoust.  Soc.  Am.  42:*^  18-521. 


Distribution  List 


Assistant  Secretary  of  the  Navy 
Research,  Development  &  Acquisition 
Navy  Department 
Washington  DC  20350-1000 

Chief  of  Naval  Operations 
Navy  Department 
Washington  DC  20350-2000 
Attn  0P71 
OP-987 

Chief  of  Naval  Operations 
Oceanographer  of  the  Navy 
U.S.  Naval  Observatory 
34th  &  Massachusetts  Ave.  NW 
Washington  DC  20392-1800 
Attn:  OP-096 
OP-0961B 

David  W.  Taylor  Naval  Research  Center 
Bethesda  MD  20084-5000 
Attn:  Commander 

Fleet  Antisub  Warfare  Tng  Ctr-All 
Naval  Station 
Norfolk  VA  23511-6495 
Attn:  Commanding  Officer 

Naval  Oceanographic  Office 
Stennis  Space  Center  MS  39522-5001 
Attn*  Commanding  Officer 

Naval  Oceanography  Command 
Stennis  Space  Center  MS  39529-5000 
Attn:  Commander 


Naval  Oceanographic  &  Atmospheric 
Research  Laboratory 
Stennis  Space  Center  MS  39529-5004 
Attn;  Code  100 
Code  105 
Code  115 
Code  125L  (10) 

Code  125P 
Code  200 
Code  300 

Naval  Ocean  Systems  Center 
San  Diego  CA  92152-5000 
Attn:  Commander 

Naval  Postgraduate  School 
Monterey  CA  93943 
Attn:  Superintendent 

Naval  Research  Laboratory 
Washington  DC  20375 
Attn:  Commanding  Officer 
Library,  M.  Bradley  (3) 

Code  TD,  L.Bernard 

Naval  Underwater  Systems  Center 
Newport  Rl  02841-5047 
Attn.  Commander 

Naval  UndenA/ater  Systems  Center  Det 
New  London  Laboratory 
New  London  CT  06320 
Attn:  Officer  in  Charge 


Office  of  Naval  Research 
800  N.  Quincy  St. 

Arlington  VA  22217-5000 
Attn:  Code  10D/10P,  Dr  E.  Silva 
Code  112,  Dr  E  Hartwig 
Code  12 
Code  10 

Office  of  Naval  Research 
ONR  European  Office 
PSC  802  Box  39 
FPO  AE  09499-0700 
Attn:  Commanding  Officer 

Office  of  Naval  Technology 
800  N.  Quincy  St 
Arlington  VA  22217-5000 
Attn;  Code  20.  Dr  P.  Selwyn 
Code  228,  Dr.  M.  Briscoe 
Code  234,  Dr.  C.  Votaw 

Scripps  Institution  of  Oceanography 
University  of  California 
291  Rosecrans  St 
San  Diego  CA  92106-3505 

Space  &  Naval  Warfare  Sys  Com 
Director  of  Navy  Laboratories 
SPAWAR  005 

Washington  DC  20363-5100 
Attn:  Commander 

Woods  Hole  Oceanographic  Institution 
P.O.  Box  32 
Woods  Hole  MA  02543 
Attn*  Director 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No  0704-0188 

Public  reporting  burden  for  th>s  coiiectic  information  is  estimated  to  ifverage  t  t>our  per  response  including  ttie  ttme  tor  reviewing  .nsiruwtions  searcbmg  existing  data  sources 
gathering  and  maintaining  the  data  neewd  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  otner  aspect  of 
this  collection  of  information  including  suggestions  for  reducing  this  burden  to  Washington  HeadQjarters  Services  Directorate  for  information  Operations  and  Reports  1215  Jefferson 
Davis  Highway  Suite  120^  Arlington  VA  22202-4302  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503 

1  Agency  Use  Only  (Leave  blank)  2.  Report  Date.  3.  Report  Type  and  Dates  Covered. 

October  1991  Final 

4.  Title  and  Subtitle. 

Theory  and  Application  of  Scattering  from  an  Object  in  an  Ocean  Waveguide 

5.  Funding  Numbers. 

Job  Order  No  12211B 

Program  £/emenf  No  0601153N 

Protect  No  03202 

Task  No  350 

Accessiof}  No  DN255011 

6.  Author(s). 

Guy  V.  Norton 

7.  Performing  Organization  Name(s)  and  Address(e$). 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory 

Ocean  Science  Directorate 

Stennis  Space  Center,  Mississippi  39529*5004 

8.  Performing  Organization 

Report  Number. 

NOARL  Report  28 

9.  Sponsoring/M  on:toring  Agency  Name(s)  and  Address(es). 

10.  Sponsoring/Monitoring  Agency 
Report  Number. 

!  11.  Supplementary  Notes. 

12a.  DIstribution/Avaiiability  Statement. 

Approved  for  public  release;  distribution  is  unlimited  Naval  Oceanographic  and 
Atmospheric  Research  Laboratory,  Stennis  Space  Center,  Mississippi  39529-5004. 

12b.  Distribution  Code. 

13.  Abstract  (Maximum  200  words).  : 

A  method  was  developed  to  describe  acoustical  scattering  from  an  object  in  a  waveguide  by  using  normal  mode  theory  to  describe 
the  incident  field.  Each  mode  is  decomposed  at  the  object  into  a  pair  of  upgoing  and  downgoing  plane  waves.  A  transition  matrix  was 
used  (developed  via  the  extended  boundary  condiuon  method  of  Waterman)  to  determine  the  resulting  near-field  scattered  field.  The 
far*field  scattered  field  was  determined  by  invoking  Huygens’  principle.  This  far-field  solution  satisfies  ail  boundary  conditions  and 
preserves  continuity  of  the  solution  throughout  all  space. 

The  examples  show  that  the  object  is  correctly  coupled  to  the  waveguide.  This  was  done  by  showing  that  the  object’s  scattered 
field  acts  as  a  secondary  source  and  that  this  scattered  field  obeyed  the  same  boundary  conditions  as  the  point  source  field. 

This  method  of  determining  the  scattered  field  from  a  three-dimensional  object  allows  one  to  determine  not  only  the  correct  target 
strength  (intensity),  but  also  to  properly  determine  the  phase.  In  addition,  this  method  allows  the  investigation  of  the  interaction  between 
the  incident  field  and  the  waveguide,  the  incident  field  with  the  object,  and  the  object’s  scattered  field  with  the  waveguide.  In  this  man¬ 
ner,  a  better  understanding  of  the  ongoing  physical  processes  can  be  obtained. 


14.  Subject  Terms. 

acoustic  scattering,  shallow  water,  waveguide  propagahon 


17.  Security  Classification 

18.  Security  Classification 

IS.  Security  Classification 

20.  Limitation  of  Abstract.  1 

of  Report 

of  This  Page. 

of  Abstract. 

I 

Unclassified 

Unclassified 

Unclassified 

None 

15.  Number  of  Pages. 
_ 252 

16.  Price  Code. 


NSN  ^540-01 -280-5 >00 


Stan/Ja'd  Form  298  (Rev  2-89i 


